Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,581 C→T 32.4% 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,5810CT32.4% 46.2 / 29.9 37V65V (GTC→GTTaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (11/14);  new base T (8/4);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

AGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGG  >  W3110S.gb/70515‑70648
                                                                  |                                                                   
aGGTTATATTCTCAATCTCACCATTCGCGGGCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGcc                                                                 <  1:761217/71‑1 (MQ=255)
  gTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGa                                                               >  1:1034158/1‑71 (MQ=255)
  gTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGa                                                               >  1:1022208/1‑71 (MQ=255)
    tataTTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGAcc                                                             >  1:485180/1‑71 (MQ=255)
     atatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCg                                                            >  1:319593/1‑71 (MQ=255)
     atatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCg                                                            >  1:1456487/1‑71 (MQ=255)
      tatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCgg                                                           <  1:913129/71‑1 (MQ=255)
          ctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGa                                                       >  1:810810/1‑71 (MQ=255)
             aaTCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGatat                                                    >  1:959592/1‑71 (MQ=255)
              aTCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAtt                                                   >  1:384374/1‑71 (MQ=255)
                 tcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTg                                                >  1:951628/1‑71 (MQ=255)
                   aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGc                                               <  1:1325420/70‑1 (MQ=255)
                    ccATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGc                                               <  1:1641351/69‑1 (MQ=255)
                    ccATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTg                                             <  1:1342734/71‑1 (MQ=255)
                         cgcgGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGcc                                                                 >  1:1319953/1‑46 (MQ=255)
                            ggTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGc                                      <  1:511280/70‑1 (MQ=255)
                                aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGga                                 <  1:291527/71‑1 (MQ=255)
                                aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                   >  1:1406825/1‑69 (MQ=255)
                                  ggggTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgaga                               <  1:555879/71‑1 (MQ=255)
                                  ggggTGGTGAAAAATCAGGGACGAGAATTTCTCTGCGGAcc                                                             <  1:958117/41‑1 (MQ=38)
                                    ggtggtGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagaga                             <  1:1741872/71‑1 (MQ=255)
                                      tggtgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGAtt                           >  1:1035484/1‑71 (MQ=255)
                                       ggtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTccc                                        >  1:576938/1‑57 (MQ=255)
                                         tgAAATATCAGGTACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGAt                            <  1:1715750/67‑1 (MQ=255)
                                           aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagaga                             >  1:164175/1‑64 (MQ=255)
                                           aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                        <  1:239527/69‑1 (MQ=255)
                                            aaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGc                                               <  1:1038253/45‑1 (MQ=255)
                                             aaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                      >  1:1771482/1‑69 (MQ=255)
                                                  aGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTg                                             <  1:389325/41‑1 (MQ=255)
                                                  aGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgt               >  1:1121802/1‑71 (MQ=255)
                                                     gACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtcg             >  1:1397119/1‑70 (MQ=255)
                                                      aCGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagag                              <  1:158951/52‑1 (MQ=255)
                                                       cGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                        <  1:1535171/57‑1 (MQ=255)
                                                          gaATTTGTCTGCCGCCCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATcc       >  1:39680/1‑71 (MQ=255)
                                                             tttgtttgCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGa    <  1:970273/71‑1 (MQ=255)
                                                             tttGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCgg     <  1:461208/70‑1 (MQ=255)
                                                             tttGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCgg     <  1:1753072/70‑1 (MQ=255)
                                                               tGTTTGACGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAgg  >  1:2167/1‑71 (MQ=255)
                                                                  |                                                                   
AGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGG  >  W3110S.gb/70515‑70648

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: