Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,581 C→T 61.5% 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,5810CT61.5% 37.4 / 58.7 52V65V (GTC→GTTaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (6/14);  new base T (24/8);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 3.35e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGG  >  W3110S.gb/70511‑70648
                                                                      |                                                                   
tGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTc                                                                     <  1:1439036/71‑1 (MQ=255)
       ttATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGAc                                                              <  1:1646301/71‑1 (MQ=255)
       ttATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGAc                                                              <  1:690993/71‑1 (MQ=255)
              ctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGc                                                                  >  1:1227793/1‑60 (MQ=255)
              ctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGc                                                                  >  1:675103/1‑60 (MQ=255)
              ctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTg                                                        <  1:305137/70‑1 (MQ=255)
                     tcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCggg                                                          >  1:1598033/1‑61 (MQ=255)
                       aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAttt                                                  <  1:1813671/67‑1 (MQ=255)
                       aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAttt                                                  <  1:622253/67‑1 (MQ=255)
                             cgcgGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTccc                                        >  1:648787/1‑71 (MQ=255)
                             cgcgGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTccc                                        >  1:758752/1‑71 (MQ=255)
                                   cAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                   >  1:1702042/1‑70 (MQ=255)
                                    aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                   >  1:1288297/1‑69 (MQ=255)
                                    aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                   >  1:1236723/1‑69 (MQ=255)
                                    aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                   >  1:1466535/1‑69 (MQ=255)
                                    aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                   >  1:1522131/1‑69 (MQ=255)
                                    aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                   >  1:250541/1‑69 (MQ=255)
                                    aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTccagccag                                   >  1:1401824/1‑69 (MQ=255)
                                      ggggTGGTGAAAAATCAGGGACGAGAATTTGTTTGcc                                                                 >  1:1015614/1‑37 (MQ=255)
                                      ggggTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGga                                 >  1:1646732/1‑69 (MQ=255)
                                      ggggTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgaga                               <  1:1273212/71‑1 (MQ=255)
                                        ggtggtGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGCTGATATTTTGCTGTTCCCGCCAGga                                 >  1:1627924/1‑67 (MQ=255)
                                          tggtgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCTCCAGGAGAGAtt                           <  1:1066180/71‑1 (MQ=255)
                                            gtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAttt                                                  <  1:378868/46‑1 (MQ=255)
                                            gtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTg                                                <  1:1263713/48‑1 (MQ=255)
                                             tgAAATATCAGGGACGAGTATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                        <  1:340069/71‑1 (MQ=255)
                                             tgAAAAATCAGGGACGAGTATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                        <  1:1657994/71‑1 (MQ=255)
                                             tgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTg                                                <  1:368577/47‑1 (MQ=255)
                                              gAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                        >  1:945439/1‑70 (MQ=255)
                                                aaaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCAc                     <  1:217416/71‑1 (MQ=255)
                                                aaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCAc                     <  1:916219/71‑1 (MQ=255)
                                                aaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCAc                     <  1:1811360/71‑1 (MQ=255)
                                                    tCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACg                 <  1:1764504/71‑1 (MQ=255)
                                                    tCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACg                 >  1:778899/1‑71 (MQ=255)
                                                    tCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACg                 <  1:159358/71‑1 (MQ=255)
                                                    tCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACg                 >  1:1155269/1‑71 (MQ=255)
                                                    tCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACg                 >  1:1234267/1‑71 (MQ=255)
                                                      aGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                      >  1:1781475/1‑64 (MQ=255)
                                                      aGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                      >  1:38228/1‑64 (MQ=255)
                                                               aaTTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgat                                <  1:101459/45‑2 (MQ=255)
                                                               aaTTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgag                                <  1:520911/45‑1 (MQ=255)
                                                               aaTTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgag                                <  1:124463/45‑1 (MQ=255)
                                                                aTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCgg     <  1:1430563/71‑1 (MQ=255)
                                                                 tttGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                     <  1:470023/38‑1 (MQ=255)
                                                                   tGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAgg  >  1:49957/1‑71 (MQ=255)
                                                                   tGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAgg  >  1:1226812/1‑71 (MQ=255)
                                                                   tGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAgg  >  1:1804712/1‑71 (MQ=255)
                                                                   tGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAgg  >  1:203486/1‑71 (MQ=255)
                                                                   tGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAg   >  1:432455/1‑70 (MQ=255)
                                                                   tGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAg   >  1:836417/1‑70 (MQ=255)
                                                                   tGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                   >  1:727905/1‑38 (MQ=255)
                                                                   tGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                   >  1:1694051/1‑38 (MQ=255)
                                                                   tGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                      >  1:936590/1‑51 (MQ=255)
                                                                      |                                                                   
TGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGG  >  W3110S.gb/70511‑70648

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