Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,581 C→T 57.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,5810CT57.4% 19.5 / 59.2 47V65V (GTC→GTTaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (5/15);  new base T (16/11);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 3.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

TGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGG  >  W3110S.gb/70511‑70645
                                                                      |                                                                
tGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTc                                                                  >  1:573178/1‑71 (MQ=255)
tGAAAGGTTATATTCTCAATCTCAACATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTc                                                                  >  1:2153179/1‑71 (MQ=255)
  aaaGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAAtttgtttg                                                                >  1:1215911/1‑71 (MQ=255)
        tataTTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGAcc                                                          >  1:2597581/1‑71 (MQ=255)
        tataTTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGAcc                                                          >  1:2486532/1‑71 (MQ=255)
        tataTTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGAcc                                                          >  1:1248170/1‑71 (MQ=255)
         atatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCg                                                         >  1:831568/1‑71 (MQ=255)
         atatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCg                                                         >  1:828559/1‑71 (MQ=255)
              ctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGa                                                            >  1:2555035/1‑63 (MQ=255)
              ctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATAAGGGACGAGAATTTGTTTGCAGACCGGGTGa                                                    >  1:2571103/1‑71 (MQ=255)
                    ctcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATAtttt                                              >  1:1789722/1‑71 (MQ=255)
                    ctcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATAttt                                               >  1:926318/1‑70 (MQ=255)
                     tcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTg                                             >  1:1244956/1‑71 (MQ=255)
                     tcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTg                                             >  1:1326997/1‑71 (MQ=255)
                     tcACCATTCACGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTg                                             >  1:2070425/1‑71 (MQ=255)
                       aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAAtttgtttg                                                                >  1:1159350/1‑50 (MQ=255)
                       aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAAtttgtttg                                                                >  1:2160239/1‑50 (MQ=255)
                       aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGc                                            <  1:2119456/70‑1 (MQ=255)
                       aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGc                                            <  1:184221/70‑1 (MQ=255)
                        ccATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTg                                          <  1:1510445/71‑1 (MQ=255)
                         cATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTg                                          <  1:2504002/70‑1 (MQ=255)
                           ttCGCGGTCAGTGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGata                                                  <  1:628512/60‑1 (MQ=255)
                           ttCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGata                                                  <  1:1923065/60‑1 (MQ=255)
                           ttCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGata                                                  <  1:707133/60‑1 (MQ=255)
                            tCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGata                                                  <  1:1359364/59‑1 (MQ=255)
                                ggTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                  <  1:1275462/71‑1 (MQ=255)
                                       gggTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCt                                           >  1:1066270/1‑55 (MQ=255)
                                       gggTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCt                                           >  1:1020189/1‑55 (MQ=255)
                                        ggtggtGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagaga                          <  1:423344/71‑1 (MQ=255)
                                        ggtggtGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagaga                          <  1:2243473/71‑1 (MQ=255)
                                         gtggtgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGtt                                        <  1:271408/56‑1 (MQ=255)
                                         gtggtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagaga                          <  1:1312698/70‑1 (MQ=255)
                                          tggtgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                  <  1:1281210/61‑1 (MQ=255)
                                          tggtgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGGGATATTTTGCTGTTCCCGcc                                  <  1:811983/61‑1 (MQ=255)
                                          tggtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGAtt                        <  1:1182936/71‑1 (MQ=255)
                                             tgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                     <  1:1010996/71‑1 (MQ=255)
                                             tgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtca                      <  1:1122053/70‑1 (MQ=255)
                                             tgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAtttt                                              <  1:2402927/46‑1 (MQ=255)
                                             tgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAtttt                                              <  1:2117810/46‑1 (MQ=255)
                                              gAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                     >  1:313783/1‑70 (MQ=255)
                                                aaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCAc                  <  1:1438647/71‑1 (MQ=255)
                                                 aaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTcc                                      >  1:72712/1‑50 (MQ=255)
                                                 aaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtca                      >  1:127742/1‑66 (MQ=255)
                                                 aaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtca                      >  1:863691/1‑66 (MQ=255)
                                                      aGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                   >  1:861519/1‑64 (MQ=255)
                                                      aGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACgg             >  1:2654870/1‑70 (MQ=255)
                                                            gagaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTcc                                      >  1:1006711/1‑39 (MQ=255)
                                                               aaTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGc                                   <  1:2086455/39‑1 (MQ=255)
                                                               aaTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCg   <  1:141195/71‑1 (MQ=255)
                                                                aTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCgg  <  1:29830/71‑1 (MQ=255)
                                                                aTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCgg  <  1:133064/71‑1 (MQ=255)
                                                                aTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCgg  <  1:952293/71‑1 (MQ=255)
                                                                      |                                                                
TGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGG  >  W3110S.gb/70511‑70645

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: