Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,458,921 T→C 42.6% K6K (AAA→AAG nusG ← transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,458,9210TC42.6% 30.0 / 85.2 68K6K (AAA→AAGnusGtranscription termination factor
Reads supporting (aligned to +/- strand):  ref base T (27/12);  new base C (14/15);  total (41/27)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTCACGCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGAATACCAT  >  W3110S.gb/3458853‑3458989
                                                                    |                                                                    
tGCTCACGCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGcttctt                                                                    <  1:4332719/71‑1 (MQ=255)
tGCTCACGCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCtttt                                                                     <  1:464384/70‑1 (MQ=255)
  cTCACGCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCttttt                                                                    >  1:3877885/1‑69 (MQ=255)
  cTCACGCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCttttt                                                                    >  1:1562635/1‑69 (MQ=255)
     aCGCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAg                                                                  >  1:3137916/1‑68 (MQ=255)
     aCGCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAg                                                                  >  1:3131710/1‑68 (MQ=255)
       gCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAg                                                                  >  1:3393248/1‑66 (MQ=255)
       gCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAg                                                                  >  1:4373115/1‑66 (MQ=255)
        cAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCtt                                                            <  1:2785390/71‑1 (MQ=255)
            gACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGa                                                        >  1:835629/1‑71 (MQ=255)
             aCGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGAc                                                       <  1:2920771/71‑1 (MQ=255)
                  gCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTa                                                                   >  1:3742858/1‑54 (MQ=255)
                  gCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTa                                                                   >  1:5073908/1‑54 (MQ=255)
                  gCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAg                                                                  >  1:4147067/1‑55 (MQ=255)
                  gCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAg                                                                  >  1:2611097/1‑55 (MQ=255)
                  gCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACAtct                                                   >  1:3919600/1‑70 (MQ=255)
                  gCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACAtct                                                   >  1:3483874/1‑70 (MQ=255)
                  gCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACAtct                                                   >  1:1539442/1‑70 (MQ=255)
                  gCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACAtct                                                   >  1:4348803/1‑70 (MQ=255)
                      cgcgGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGaa                                              >  1:1840327/1‑71 (MQ=255)
                      cgcgGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGaa                                              <  1:1027038/71‑1 (MQ=255)
                      cgcgGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGa                                               <  1:3468602/70‑1 (MQ=255)
                      cgcgGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAg                                                <  1:4108184/69‑1 (MQ=255)
                      cgcgGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGaa                                              >  1:3326013/1‑71 (MQ=255)
                      cgcgGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGa                                               <  1:3003086/70‑1 (MQ=255)
                           ccTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCt                                           >  1:3352055/1‑69 (MQ=255)
                            cTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAg                                                         <  1:4348064/54‑1 (MQ=255)
                             ttCAAAACCGGAAAACGCCTGAACGACGTACCAGCGcttctt                                                                    <  1:4327800/42‑1 (MQ=255)
                             ttCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCa                                                          <  1:1950566/52‑1 (MQ=255)
                              tCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGAc                                                       >  1:1791910/1‑54 (MQ=255)
                              tCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGc                                      >  1:4325674/1‑71 (MQ=255)
                              tCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGc                                      >  1:4746088/1‑71 (MQ=255)
                              tCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGc                                      >  1:1183883/1‑71 (MQ=255)
                              tCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGc                                      >  1:2728983/1‑71 (MQ=255)
                               cAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAg                                                                  >  1:4917658/1‑42 (MQ=255)
                               cAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGc                                      >  1:4711135/1‑70 (MQ=255)
                               cAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAgg                                                                 >  1:634743/1‑43 (MQ=255)
                               cAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAg                                                               >  1:4186903/1‑45 (MQ=255)
                               cAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAg                                                               >  1:4635951/1‑45 (MQ=255)
                               cAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGa                                                        >  1:3791125/1‑52 (MQ=255)
                                aaaaCCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTa                                                                   >  1:2378759/1‑40 (MQ=255)
                                aaaaCCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGcc                                     <  1:3615988/70‑1 (MQ=255)
                                aaaaCCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGcc                                     <  1:3645875/70‑1 (MQ=255)
                                aaaaCCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGc                                      >  1:3942340/1‑69 (MQ=255)
                                   aCCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTg                                 <  1:2063311/71‑1 (MQ=255)
                                   aCCGGAAAACGCCTGAACGACGTACCAGCGCTTCTTAgg                                                                 <  1:4271497/39‑1 (MQ=255)
                                      ggAAAACGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATa                              <  1:1915742/71‑1 (MQ=255)
                                         aaaCGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGAt                               >  1:1592595/1‑67 (MQ=255)
                                         aaaCGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAg                           <  1:3340241/71‑1 (MQ=255)
                                         aaaCGCCTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAg                           >  1:2668778/1‑71 (MQ=255)
                                          aaCGCCAGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAg                           >  1:4584/1‑70 (MQ=255)
                                              ccTGAACGACGTACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAg                           <  1:4832532/66‑1 (MQ=255)
                                                tGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATAcc                     <  1:2634571/70‑1 (MQ=255)
                                                      acgTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCt                                           >  1:3106490/1‑42 (MQ=255)
                                                      acgTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGaa              >  1:3248535/1‑71 (MQ=255)
                                                      acgTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGaa              >  1:1809168/1‑71 (MQ=255)
                                                         tACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCg                <  1:4700710/66‑1 (MQ=255)
                                                         tACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCa           <  1:62103/71‑1 (MQ=255)
                                                         tACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCa           <  1:1095103/71‑1 (MQ=255)
                                                         tACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCa           <  1:4006862/71‑1 (MQ=255)
                                                         tACCAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCa           <  1:104278/71‑1 (MQ=255)
                                                           ccAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATa                       >  1:179340/1‑57 (MQ=255)
                                                           ccAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGc                 >  1:3382160/1‑63 (MQ=255)
                                                           ccAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAg          <  1:5144942/70‑1 (MQ=255)
                                                            ccgcgcTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAg          <  1:3371355/67‑1 (MQ=255)
                                                            cAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGa                         >  1:4246580/1‑54 (MQ=255)
                                                            cAGCGCTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAg          <  1:5139877/69‑1 (MQ=255)
                                                             aGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGaa        >  1:4047698/1‑70 (MQ=255)
                                                               cgcTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGAATa      >  1:4842895/1‑70 (MQ=255)
                                                                gcTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGAATAcc    <  1:4622559/71‑1 (MQ=255)
                                                                gcTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGAATAcc    <  1:2901116/71‑1 (MQ=255)
                                                                gcTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGAATAcc    <  1:3608590/71‑1 (MQ=255)
                                                                gcTTCTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGAATAc     >  1:754848/1‑70 (MQ=255)
                                                                    tttAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGAATACCAt  <  1:103448/69‑1 (MQ=255)
                                                                    |                                                                    
TGCTCACGCAGCGACGTTGCTACGCGGCCTTCAAAACCGGAAAACGCCTGAACGACGTACCAGCGCTTTTTAGGAGCTTCAGACATCTCAGAACCTCAGGCCAGTGATAAAGGATACCAGGCGAACCAGAATACCAT  >  W3110S.gb/3458853‑3458989

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: