Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,156,112 T→C 100% A117A (GCT→GCC holB → DNA polymerase III, delta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,156,1120TC89.5% 156.1 / 8.6 57A117A (GCT→GCCholBDNA polymerase III, delta prime subunit
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base C (27/24);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

AAATACGCTGGGCGTTGATGCGGTACGTGAGGTCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCTGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCC  >  NC_000913/1156016‑1156204
                                                                                                |                                                                                            
aaaTACGCTGGGCGTTGATGCGGTACGTGAGGTCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATgccgcc                                                                                           >  1:1326090/1‑100 (MQ=255)
     cGCTGGGCGTTGATGCGGTACGTGAGGTCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACt                                                                                      >  1:2345824/1‑100 (MQ=255)
       cTGGGCGTTGATGCGGTACGTGAGGTCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTaa                                                                                    >  2:1527883/1‑100 (MQ=255)
         gggCGTTGATGCGGTACGTGAGGTCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAAcc                                                                                  <  1:2545889/100‑1 (MQ=255)
            cGTTGATGCGGTACGTGAGGTCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGAc                                                                               <  1:237833/100‑1 (MQ=255)
            cGTTGATGCGGTACGTGAGGTCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGAc                                                                               <  2:1316872/100‑1 (MQ=255)
                  tGCGGTACGTGAGGTCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCTGCCTTACTAACCGACGCcgcg                                                                         <  1:796563/100‑1 (MQ=255)
                  tGCGGTACGTGAGGTCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCcgcg                                                                         >  1:2057241/1‑100 (MQ=255)
                        aCGTGAGGTCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAAc                                                                   <  1:2234367/100‑1 (MQ=255)
                         cGTGGGGTCACCGAAAAGCTGAATGAGCACGCACGTTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGATGTCGCGGCTAAc                                                                   <  2:637929/99‑1 (MQ=255)
                           ttaGGTCACCGAAATGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACg                                                                  <  2:206497/96‑1 (MQ=255)
                              ggTCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTg                                                             <  2:1673750/100‑1 (MQ=255)
                              ggTCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTg                                                             >  2:1120956/1‑100 (MQ=255)
                                tCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCt                                                           <  2:1998598/100‑1 (MQ=255)
                                tCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCt                                                           >  1:946975/1‑100 (MQ=255)
                                 cACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGc                                                            >  2:370331/1‑98 (MQ=255)
                                   ccGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGaa                                                        >  1:1974251/1‑100 (MQ=255)
                                   ccGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGaa                                                        >  2:262834/1‑100 (MQ=255)
                                    cGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGaaa                                                       <  1:1337816/100‑1 (MQ=255)
                                    cGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGaaa                                                       <  1:1248235/100‑1 (MQ=255)
                                      aaaaGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAAc                                                     >  2:2490488/1‑100 (MQ=255)
                                         aGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCTGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCt                                                  <  1:1909134/100‑1 (MQ=255)
                                          gCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCtt                                                 >  2:231730/1‑100 (MQ=255)
                                           cTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTg                                                <  2:1080907/100‑1 (MQ=255)
                                               aTGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAagn                                            <  2:1857967/100‑2 (MQ=255)
                                               aTGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAaga                                            <  1:1017193/100‑1 (MQ=255)
                                                tGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCTGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAagag                                           <  1:1291939/100‑1 (MQ=255)
                                                tGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAagag                                           >  2:2579940/1‑100 (MQ=255)
                                                tGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAagag                                           <  2:2345811/100‑1 (MQ=255)
                                                   gcacgcacGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGcc                                         <  2:1502693/99‑1 (MQ=255)
                                                     acgcacgcTAAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAaccacc                                      >  2:861100/1‑100 (MQ=255)
                                                        cacgcTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACcagc                                   >  1:698706/1‑100 (MQ=255)
                                                           gcTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGa                                >  1:1723966/1‑100 (MQ=255)
                                                              tAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAAc                             <  1:1464076/100‑1 (MQ=255)
                                                                ggtggtGCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACtt                           >  2:816026/1‑100 (MQ=255)
                                                                 gtggtgCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTg                          <  2:2067274/100‑1 (MQ=255)
                                                                   ggtgCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGt                        >  2:1946136/1‑100 (MQ=255)
                                                                    gtgCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGtt                       >  1:234225/1‑100 (MQ=255)
                                                                     tgCGAAAGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGttt                      <  2:1433836/100‑1 (MQ=255)
                                                                       cGAAAGTCGTTTGGGTAACCGATGCTGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGttttt                    >  1:1811759/1‑100 (MQ=255)
                                                                         aaaGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCAGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTcc                  <  1:2505670/100‑1 (MQ=255)
                                                                           aGTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGTATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTg                <  2:770363/100‑1 (MQ=255)
                                                                            gTCGTTTGGGTAACCGATGCTGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTgg               <  2:2141511/100‑1 (MQ=255)
                                                                            gTCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCGGAAACTTGGTTTTTCCTgg               <  2:1178153/100‑1 (MQ=255)
                                                                             tCGTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGc              <  1:330792/100‑1 (MQ=255)
                                                                               gTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTa            >  1:446507/1‑100 (MQ=255)
                                                                               gTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTa            >  1:311093/1‑100 (MQ=255)
                                                                               gTTTGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTa            >  1:311044/1‑100 (MQ=255)
                                                                                  tGGGTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTAccc         >  1:2042050/1‑100 (MQ=255)
                                                                                   gggTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCcg        >  1:2276254/1‑100 (MQ=255)
                                                                                   gggTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCcg        <  1:837638/100‑1 (MQ=255)
                                                                                    ggTAACCGATGCTGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCcgc       >  1:2273666/1‑100 (MQ=255)
                                                                                    ggTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCcgc       >  2:127908/1‑100 (MQ=255)
                                                                                    ggTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCcgc       >  2:1333327/1‑100 (MQ=255)
                                                                                     gTAACCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCcgcg      >  1:1727031/1‑100 (MQ=255)
                                                                                        aCCGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGc   <  1:1122090/100‑1 (MQ=255)
                                                                                         ccGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGcc  >  1:1611385/1‑100 (MQ=255)
                                                                                         ccGATGCCGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGcc  >  1:1006763/1‑100 (MQ=255)
                                                                                                |                                                                                            
AAATACGCTGGGCGTTGATGCGGTACGTGAGGTCACCGAAAAGCTGAATGAGCACGCACGCTTAGGTGGTGCGAAAGTCGTTTGGGTAACCGATGCTGCCTTACTAACCGACGCCGCGGCTAACGCATTGCTGAAAACGCTTGAAGAGCCACCAGCAGAAACTTGGTTTTTCCTGGCTACCCGCGAGCC  >  NC_000913/1156016‑1156204

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