Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,800,993 G→A 10.2% R326H (CGC→CAC)  dnaEA → DNA polymerase III subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,800,9930GA10.2% 101.8 / 2.3 39R326H (CGC→CAC) dnaEADNA polymerase III subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (17/18);  new base A (2/2);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

ACTTCCCCACGCCCAACGGCATGGGCATCGACGACTACCTGCGGCATGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAACGACGTACCGGTGGGCCCGGGTCGTGGTTCGGGTGCCGGCTCGCTGGTGGCCTATGT  >  NC_002947/1800852‑1801140
                                                                                                                                             |                                                                                                                                                   
aCTTCCCCACGCCCAACGGCATGGGCATCGACGACTACCTGCGGCATGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGCCCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAg                                                                                                                                              >  5:272567/1‑149 (MQ=255)
    ccccACGCCCAACGGCATGGGCATCGACGACTACCTGCGGCATGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTc                                                                                                                                           <  5:501/148‑1 (MQ=255)
    ccccACGCCCAACGGCATGGGCATCAACGACTACCTGCGGCATGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTc                                                                                                                                           <  5:126841/148‑1 (MQ=255)
      ccACGCCCAACGGCATGGGCATCGACGACTACCTGCGGCATGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAg                                                                                                                                        >  8:184081/1‑149 (MQ=255)
                   cATGGGCATCGACAACTACCTGCGGCATGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGAcatcatca                                                                                                                           >  8:292208/1‑149 (MQ=255)
                    aTGGGCATCGACGACTACCTGCGGCATGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGAcatcatcat                                                                                                                          >  1:272291/1‑149 (MQ=255)
                        gCATCGACGACTACCTGCGGCATGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACagcatcatc                                                                                                                         >  8:328923/1‑146 (MQ=255)
                         cATCGACGACTACCTGCGGCATGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGa                                                                                                                     <  8:284074/149‑1 (MQ=255)
                                       tGCGGCATGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGAcatc                                                                                                                               >  7:268770/1‑125 (MQ=255)
                                       tGCGGCATGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGAcatc                                                                                                                               <  8:268770/125‑1 (MQ=255)
                                        gCGGCATGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGt                                                                                                       >  7:17886/1‑148 (MQ=255)
                                        gCGGCATGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGt                                                                                                       >  6:88176/1‑148 (MQ=255)
                                             aTGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTc                                                                                                 <  2:272291/149‑1 (MQ=255)
                                              tgtgGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTcc                                                                                                <  8:17886/149‑1 (MQ=255)
                                                tgGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTg                                                                                              <  5:88176/149‑1 (MQ=255)
                                                 ggCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGGGGCCGAAAGACACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCGGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGa                                                                                             <  6:272567/149‑1 (MQ=255)
                                                  gCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGa                                                                                             <  7:292208/148‑1 (MQ=255)
                                                  gCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGAt                                                                                            <  5:347225/149‑1 (MQ=255)
                                                         aaGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCTGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATg                                                                                     <  1:7251/149‑1 (MQ=255)
                                                              cTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACt                                                                                 >  7:208678/1‑148 (MQ=255)
                                                                    gagCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTtcatca                                                                           >  3:86162/1‑148 (MQ=255)
                                                                     agCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGTCATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCaaa                                                                         <  1:284028/149‑1 (MQ=255)
                                                                      gcgcCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTtc                                                                               <  1:168838/142‑1 (MQ=255)
                                                                      gcgcCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATCGGCTTCCCCGGTTACTTCCTCATTGGTATGGACTtc                                                                               >  2:168838/1‑142 (MQ=255)
                                                                                      tgGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTGCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAacg                                                         >  5:309138/1‑148 (MQ=255)
                                                                                      tgGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAacg                                                         <  6:309138/148‑1 (MQ=255)
                                                                                          cGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAACGACGTa                                                    >  7:74399/1‑149 (MQ=255)
                                                                                             aaGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTTGACATCATCATCCAGATTGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAACGACGTAcc                                                  <  7:184081/148‑1 (MQ=255)
                                                                                                agaCCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCACCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAacg                                                         <  2:320173/138‑1 (MQ=255)
                                                                                                agaCCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCACCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAacg                                                         <  2:315124/138‑1 (MQ=255)
                                                                                                agaCCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCACCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAacg                                                         >  1:320173/1‑138 (MQ=255)
                                                                                                agaCCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCACCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAacg                                                         >  1:315124/1‑138 (MQ=255)
                                                                                                  aCCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGACAAGAACAACGACGTACCGGTggg                                            >  3:249088/1‑149 (MQ=255)
                                                                                                  aCCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAAAGGGCCAAGAACAACGACGTACCGGTggg                                            <  1:200696/149‑1 (MQ=255)
                                                                                                             ttATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAACGACGTACCGGTGGGCCCGGGTCGTg                                 <  4:86162/149‑1 (MQ=255)
                                                                                                                       aaGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAACGACGTACCGGTGGGCCCGGGTCGTGGTTCGGGTGc                       >  3:61590/1‑149 (MQ=255)
                                                                                                                               ggTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAACGACGTACCGGTGGGCCCGGGTCGTGGTTCGGGTGCCGGCTCGc               >  3:331115/1‑149 (MQ=255)
                                                                                                                               ggTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAACGACGTACCGGTGGGCCCGGGTCGTGGTTCGGGTGCCGGCTCGc               >  6:16301/1‑149 (MQ=255)
                                                                                                                                            cGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAACGACGTACCGGTGGGCCCGGGTCGTGGTTCGGGTGCCGGCTCGCTGGTGGCCTATGt  <  8:208678/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
ACTTCCCCACGCCCAACGGCATGGGCATCGACGACTACCTGCGGCATGTGGCCCACGAAGGCCTGGAAGAGCGCCTGGCCGTGCTGTGGCCGAAAGAGACCACGCCCAATTATGAAGAAAAGCGTCAGGTTTACCTGGACCGCCTGAAGTTCGAGCTGGACATCATCATCCAGATGGGCTTCCCCGGTTACTTCCTGATTGTTATGGACTTCATCAAATGGGCCAAGAACAACGACGTACCGGTGGGCCCGGGTCGTGGTTCGGGTGCCGGCTCGCTGGTGGCCTATGT  >  NC_002947/1800852‑1801140

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