Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 236,541 C→T R159* (CGA→TGA)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913236,5410CT94.6% 172.8 / 0.6 56R159* (CGA→TGA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (32/21);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGA  >  NC_000913/236445‑236635
                                                                                                |                                                                                              
gTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaa                                                                                             >  1:917061/1‑100 (MQ=255)
  cACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaacg                                                                                           >  2:212273/1‑100 (MQ=255)
   aCCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGc                                                                                          <  2:1197803/100‑1 (MQ=255)
     cGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTg                                                                                        >  2:1917479/1‑100 (MQ=255)
           ccTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACggg                                                                                  >  1:1756304/1‑100 (MQ=255)
                cGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCAtt                                                                             >  2:1103968/1‑100 (MQ=255)
                cGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCAtt                                                                             <  1:1596521/100‑1 (MQ=255)
                   tGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGTGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACt                                                                          <  1:1360999/100‑1 (MQ=255)
                     gCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCg                                                                        >  2:829324/1‑100 (MQ=255)
                      cGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGGTGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGa                                                                       >  2:674575/1‑100 (MQ=255)
                         ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGc                                                                    >  1:68305/1‑100 (MQ=255)
                           aaaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGccc                                                                  <  2:1137079/100‑1 (MQ=255)
                            aaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGCTAACAGTAAATGAACGCTGCACGGGGAATTACTCGATGcac                                                                  >  1:1504084/1‑97 (MQ=255)
                            aaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCa                                                                 >  2:951451/1‑100 (MQ=255)
                            aaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCa                                                                 >  1:138443/1‑100 (MQ=255)
                             aaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAg                                                                >  1:346176/1‑100 (MQ=255)
                                 tctCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATcc                                                            <  2:194435/98‑1 (MQ=255)
                                  ttCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCt                                                           >  1:393954/1‑100 (MQ=255)
                                    cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATccaga                                                         >  1:194434/1‑97 (MQ=255)
                                    cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTg                                                         <  1:989350/100‑1 (MQ=255)
                                    cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTg                                                         >  1:1674921/1‑100 (MQ=255)
                                     ccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGc                                                        <  1:1148776/100‑1 (MQ=255)
                                      cGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCg                                                       <  1:91349/100‑1 (MQ=255)
                                          aaGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAg                                                   <  1:181060/100‑1 (MQ=255)
                                          aaGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAg                                                   >  1:1092839/1‑100 (MQ=255)
                                           aGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGACATTACTCGATGCCCAGATCCTTGCGGAAGt                                                  >  2:500695/1‑100 (MQ=255)
                                             cgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGttt                                                >  2:815337/1‑100 (MQ=255)
                                             cgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGttt                                                <  1:1359280/100‑1 (MQ=255)
                                              gcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTa                                               <  2:1848535/100‑1 (MQ=255)
                                                aaCAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGACCGCTGCACGGGGCATTACTCGATGCCCACATCCTTGCGCAAGTTTATc                                             >  1:565901/1‑100 (MQ=255)
                                                 aCAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCt                                            >  1:558225/1‑100 (MQ=255)
                                                   aGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTgg                                          >  1:1729506/1‑100 (MQ=255)
                                                    gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                         <  2:114778/100‑1 (MQ=255)
                                                       tCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGAt                                      >  2:355447/1‑100 (MQ=255)
                                                       tCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGAt                                      <  2:789089/100‑1 (MQ=255)
                                                               ttATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtg                              <  1:1334444/100‑1 (MQ=255)
                                                                tATGTGCTCGCGACGAAATAGATAACAGTAAACGAACGGTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTAGCTGGCGATGACCggtgg                             <  1:1845582/100‑1 (MQ=255)
                                                                 aTGTGCTCGCTACTAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtggt                            <  2:484211/100‑1 (MQ=255)
                                                                     gCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCaaa                        >  2:1032838/1‑100 (MQ=255)
                                                                      cTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAAc                       >  1:653164/1‑100 (MQ=255)
                                                                       tCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACg                      >  2:950809/1‑100 (MQ=255)
                                                                       tCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACg                      >  2:276732/1‑100 (MQ=255)
                                                                          cTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCg                   >  2:1700527/1‑100 (MQ=255)
                                                                           tACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGa                  >  2:1172949/1‑100 (MQ=255)
                                                                            aCGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGa                  >  1:1253729/1‑99 (MQ=255)
                                                                              gAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATgg               <  1:881616/100‑1 (MQ=255)
                                                                                aaTAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCAATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCt             >  2:1044952/1‑100 (MQ=255)
                                                                                     aTAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGc        <  1:537600/100‑1 (MQ=255)
                                                                                     aTAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGc        <  2:36157/100‑1 (MQ=255)
                                                                                     aTAACAGTAAATGAACGCGGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTg         <  1:437762/99‑1 (MQ=255)
                                                                                      tAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCg       >  1:1394212/1‑100 (MQ=255)
                                                                                       aaCAGTAAATGAACGATGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGa      <  2:1912037/100‑1 (MQ=255)
                                                                                        aCAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCg       >  1:641643/1‑98 (MQ=255)
                                                                                        aCAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGAt     <  1:674574/100‑1 (MQ=255)
                                                                                         cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg    >  2:1408697/1‑100 (MQ=255)
                                                                                          aGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATgg   >  1:1140465/1‑100 (MQ=255)
                                                                                          aGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATgg   >  2:606279/1‑100 (MQ=255)
                                                                                           gTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGa  >  2:1690584/1‑100 (MQ=255)
                                                                                                |                                                                                              
GTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGA  >  NC_000913/236445‑236635

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: