Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,109 T→C 36.8% T241T (ACT→ACC araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb71,1090TC36.8% 30.7 / 38.6 38T241T (ACT→ACCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (13/11);  new base C (6/8);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTT  >  W3110S.gb/71043‑71168
                                                                  |                                                           
aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAccc                                                         >  1:723223/1‑71 (MQ=255)
aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAcc                                                          >  1:1811410/1‑70 (MQ=255)
 gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                        <  1:1586818/71‑1 (MQ=255)
 gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                        <  1:1212758/71‑1 (MQ=255)
 gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                        <  1:480224/71‑1 (MQ=255)
 gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                        <  1:504736/71‑1 (MQ=255)
        gCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                        <  1:2195743/64‑1 (MQ=255)
           tCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                        <  1:588106/61‑1 (MQ=255)
             ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGa                                                      >  1:1373399/1‑61 (MQ=255)
             ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTa                                                            >  1:752692/1‑55 (MQ=255)
             ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                        <  1:536925/59‑1 (MQ=255)
               aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGa                                                      >  1:195739/1‑59 (MQ=255)
               aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGc                                            <  1:1084591/69‑1 (MQ=255)
               aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                                          >  1:504204/1‑71 (MQ=255)
                     gcgcGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtc                                     <  1:1362532/70‑1 (MQ=255)
                     gcgcGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgt                                      <  1:1075670/69‑1 (MQ=255)
                      cgcgAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCgtcgg                                   <  1:119649/71‑1 (MQ=255)
                        cgAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtcggtc                                 >  1:934848/1‑71 (MQ=255)
                         gAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtcggtc                                 >  1:1689097/1‑70 (MQ=255)
                           ggcccAACGCATTAGTCAGGCGAAGCTGCTTTTGAGAACTACCCGGATGCCTATCGCCACCGTCGGTCGc                               >  1:2128127/4‑70 (MQ=255)
                            gACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCaa                             >  1:926012/1‑71 (MQ=255)
                              ccAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTa                                                <  1:1206202/50‑1 (MQ=255)
                                 aCGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTAt                                               >  1:1860587/1‑48 (MQ=255)
                                 aCGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTg                        >  1:1104435/1‑71 (MQ=255)
                                 aCGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTg                        >  1:706743/1‑71 (MQ=255)
                                   gCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGt                      >  1:1249649/1‑71 (MQ=255)
                                        agccAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTg                  <  1:371229/67‑1 (MQ=255)
                                        aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGt                          <  1:462567/62‑1 (MQ=255)
                                         gTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGAc                <  1:1719171/71‑1 (MQ=255)
                                          tCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTg                        >  1:1622363/1‑62 (MQ=255)
                                           cAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGa              <  1:245708/71‑1 (MQ=255)
                                            aGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGa              <  1:1356082/70‑1 (MQ=255)
                                             ggCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGAt             <  1:1617277/70‑1 (MQ=255)
                                                    ctgctTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTa     <  1:188874/71‑1 (MQ=255)
                                                      gctTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGttt                    >  1:2024278/1‑54 (MQ=255)
                                                       ctTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAttt  >  1:1317140/1‑71 (MQ=255)
                                                         tttGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATc            >  1:1958573/1‑59 (MQ=255)
                                                          ttGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGAc                >  1:207903/1‑54 (MQ=255)
                                                                  tACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt                   <  1:1792988/43‑1 (MQ=255)
                                                                  |                                                           
AGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTT  >  W3110S.gb/71043‑71168

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: