Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,581 C→T 34.1% V65V (GTC→GTT araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb70,5810CT34.1% 51.6 / 40.4 44V65V (GTC→GTTaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (18/11);  new base T (8/7);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

AAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCT  >  W3110S.gb/70513‑70650
                                                                    |                                                                     
aaaGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAAtttgtttg                                                                     <  1:2287412/71‑1 (MQ=255)
aaaGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTg                                                                     >  1:1529699/1‑71 (MQ=255)
  aGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGcc                                                                   >  1:743003/1‑71 (MQ=255)
   ggTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGc                                                                    >  1:2927564/1‑69 (MQ=255)
      tataTTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGAcc                                                               >  1:3077608/1‑71 (MQ=255)
          ttCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCTACCggg                                                            <  1:5097/70‑1 (MQ=255)
           tctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTc                                                                       >  1:184137/1‑58 (MQ=255)
              cAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGAc                                                                >  1:1043936/1‑62 (MQ=255)
                     aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAttt                                                    <  1:2383166/67‑1 (MQ=255)
                      ccATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTg                                               >  1:2981211/1‑71 (MQ=255)
                        aTTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGtt                                             <  1:2121798/71‑1 (MQ=255)
                        aTTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGtt                                             >  1:1318711/1‑71 (MQ=255)
                         ttCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGa                                                                 >  1:1391573/1‑50 (MQ=255)
                               ttCAGGGGGGGGTGAAAAATGAGGGACGAGAATTTGTCTGCCg                                                                  <  1:2253153/42‑1 (MQ=38)
                               gTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGc                                                                    >  1:305338/1‑41 (MQ=255)
                               gTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGc                                        >  1:526167/1‑69 (MQ=255)
                                  aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGga                                   <  1:1864104/71‑1 (MQ=255)
                                  aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGt                                                           >  1:670648/1‑47 (MQ=255)
                                    ggggTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgaga                                 <  1:393141/71‑1 (MQ=255)
                                    ggggTGATGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGt                                              <  1:1190433/58‑1 (MQ=255)
                                     gggTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagag                                <  1:2904783/71‑1 (MQ=255)
                                       gttgtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                       <  1:421738/59‑1 (MQ=255)
                                           tgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                          <  1:1665628/71‑1 (MQ=255)
                                             aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAgg                                    >  1:84768/1‑59 (MQ=255)
                                             aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                        >  1:2600539/1‑71 (MQ=255)
                                             aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                          >  1:121308/1‑69 (MQ=255)
                                              aaaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                        >  1:2157045/1‑70 (MQ=255)
                                               aaaTCAGTGACGAGAATTTGTTTGCCGACCGGGTGat                                                        >  1:3044561/1‑37 (MQ=37)
                                                aaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgag                                  <  1:2185370/58‑1 (MQ=255)
                                                aaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTa                     >  1:74110/1‑71 (MQ=255)
                                                  tCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACg                   >  1:250574/1‑71 (MQ=255)
                                                  tCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACg                   <  1:1901980/71‑1 (MQ=255)
                                                    aGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATCCATCACTACGgt                 <  1:2690652/71‑1 (MQ=255)
                                                      ggACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgaga                                 <  1:1051011/53‑1 (MQ=255)
                                                        aCGAGAGTTTGTCTGCCGACCGGTTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtcgt              >  1:408439/1‑70 (MQ=255)
                                                        aCGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGAtt                             <  1:1618316/55‑1 (MQ=255)
                                                           agaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATc          <  1:892684/71‑1 (MQ=255)
                                                            gaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATcc         >  1:524831/1‑71 (MQ=255)
                                                            gaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATcc         >  1:1724926/1‑71 (MQ=255)
                                                             aaTTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCg        <  1:3112474/71‑1 (MQ=255)
                                                              aTTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCgg       <  1:2594879/71‑1 (MQ=255)
                                                               tttGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtc                >  1:228690/1‑61 (MQ=255)
                                                                 tGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                     >  1:1617214/1‑38 (MQ=255)
                                                                 tGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAgg    >  1:2261971/1‑71 (MQ=255)
                                                                  gTCTGCCGACCGGGTGATATTTTGATGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGc   >  1:1430669/1‑71 (MQ=255)
                                                                   tCTGGCGCCCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATAACTACGGTCGTCATCCGGAGGCt  >  1:2964732/1‑71 (MQ=255)
                                                                    |                                                                     
AAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCT  >  W3110S.gb/70513‑70650

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