Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,802,847 T→C 6.9% M944T (ATG→ACG)  dnaEA → DNA polymerase III subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,802,8470TC 6.9% 167.5 / 2.6 58M944T (ATG→ACG) dnaEADNA polymerase III subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base C (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TTCCATGACGAGATCAAGGCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGACACCCTTGGCCTGTACCTCACCGGCCACCCGATCGACGAGTACGAGACCGAGAT  >  NC_002947/1802705‑1802988
                                                                                                                                              |                                                                                                                                             
ttCCATGACGAGATCAAGGCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGa                                                                                                                                         >  4:18581/1‑149 (MQ=255)
  ccATGACGAGATCAAGGCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGt                                                                                                                                             >  3:156084/1‑143 (MQ=255)
  ccATGACGAGATCAAGGCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGt                                                                                                                                             <  4:156084/143‑1 (MQ=255)
   cATGACGAGATCAAGGCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGGAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGAcgc                                                                                                                                      >  8:199663/1‑149 (MQ=255)
   cATGACGAGATCAAGGCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGAcg                                                                                                                                       <  5:216738/148‑1 (MQ=255)
       aCGAGATCAAGGCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGc                                                                                                                                   <  1:189492/148‑1 (MQ=255)
             tCAAGGCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTc                                                                                                                            >  7:23971/1‑149 (MQ=255)
             tCAAGGCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTc                                                                                                                            <  4:360954/149‑1 (MQ=255)
              cAAGGCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCACGTTCGACGCGGCAGATGTc                                                                                                                            >  1:272861/1‑148 (MQ=255)
              cAAGGCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCACGTTCGACGCGGCAGATGTCg                                                                                                                           >  3:256473/1‑149 (MQ=255)
                 ggCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACg                                                                                                                        >  4:152377/1‑149 (MQ=255)
                 ggCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCACTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACg                                                                                                                        >  1:88869/1‑149 (MQ=255)
                  gCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGt                                                                                                                       <  3:18581/149‑1 (MQ=255)
                    cTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTc                                                                                                                      >  5:78320/1‑148 (MQ=255)
                     tACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGTAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGAt                                                                                                                               >  8:346361/1‑138 (MQ=255)
                     tACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGTAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGAt                                                                                                                               <  7:346361/138‑1 (MQ=255)
                         aCGCCAACATCGACATCAACCGAGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGc                                                                                                                 >  8:56862/1‑148 (MQ=255)
                                     aCATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAg                                                                                                    <  5:150824/149‑1 (MQ=255)
                                       aTCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCACGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGCCCTGTTCGGCAGCATGTTCTACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGt                                                                                                  >  4:43779/1‑149 (MQ=255)
                                       aTCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGt                                                                                                  >  8:111499/1‑149 (MQ=255)
                                         cAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTgc                                                                                                >  5:211603/1‑149 (MQ=255)
                                           aCCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTgcgc                                                                                              <  8:23971/149‑1 (MQ=255)
                                                    ccTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCAcc                                                                                     >  7:142784/1‑149 (MQ=255)
                                                          tGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTc                                                                                  >  5:339811/1‑146 (MQ=255)
                                                            tCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAg                                                                               <  8:320071/147‑1 (MQ=255)
                                                            tCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAg                                                                               >  7:132557/1‑147 (MQ=255)
                                                              ggCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGc                                                                           <  3:152377/149‑1 (MQ=255)
                                                              ggCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGc                                                                           >  8:301934/1‑149 (MQ=255)
                                                               gcgcTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCg                                                                          >  6:348416/1‑149 (MQ=255)
                                                                  cTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCt                                                                       >  5:98699/1‑149 (MQ=255)
                                                                      gCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAg                                                                   <  6:78320/149‑1 (MQ=255)
                                                                       cGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAgg                                                                  >  5:114668/1‑149 (MQ=255)
                                                                        gAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAggg                                                                 <  4:59077/149‑1 (MQ=255)
                                                                        gAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAggg                                                                 <  5:348416/149‑1 (MQ=255)
                                                                         aGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAgggg                                                                <  7:56862/149‑1 (MQ=255)
                                                                             cATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGgaga                                                            <  7:111499/149‑1 (MQ=255)
                                                                              aTCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGgag                                                             >  7:148650/1‑147 (MQ=255)
                                                                              aTCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGaa                                                           <  2:88869/149‑1 (MQ=255)
                                                                                 aaGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAgg                                                         <  6:211603/148‑1 (MQ=255)
                                                                                 aaGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGa                                                        <  7:301934/149‑1 (MQ=255)
                                                                                  aGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGac                                                       >  6:146687/1‑149 (MQ=255)
                                                                                   ggCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCAATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGaca                                                      >  2:305585/1‑149 (MQ=255)
                                                                                        ccGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGACACCCtt                                                 <  8:142784/149‑1 (MQ=255)
                                                                                                    cgcATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCACGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGACACCCTTGGCCTGTACCt                                      <  4:256473/148‑1 (MQ=255)
                                                                                                    cgcATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCACGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGACACCCTTGGCCTGTACCt                                      <  2:272861/148‑1 (MQ=255)
                                                                                                         aCCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATg                                                                                                                              <  5:384574/55‑1 (MQ=255)
                                                                                                         aCCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATg                                                                                                                              >  6:384574/1‑55 (MQ=255)
                                                                                                          ccGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGACACCCTTGGCCTGTACCTCACCGGc                               <  8:235661/149‑1 (MQ=255)
                                                                                                                aCAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACACTCAAGGAGCGTCTGAAGGGGGAGAAGGACACCCTTGGCCTGTACCTCACCGGCCAccc                          >  4:147245/1‑148 (MQ=255)
                                                                                                                   gTGGCAATGTCGACCTGTTCGGCAGCATGGTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCCAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGACACCCTTGGCCTGTACCTCACCGGCCACCCGATc                      <  3:43779/149‑1 (MQ=255)
                                                                                                                    tGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCt                                                                                   <  8:199077/87‑1 (MQ=255)
                                                                                                                    tGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCt                                                                                   >  7:199077/1‑87 (MQ=255)
                                                                                                                      gCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGACACCCTTGGCCTGTACCTCACCGGCCACCCGATcgac                   <  8:132557/149‑1 (MQ=255)
                                                                                                                          tGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGACACCCTTGGCCTGTACCTCACCGGCCACCCGATCGACGAg                >  5:7084/1‑148 (MQ=255)
                                                                                                                          tGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGACACCCTTGGCCTGTACCTCACCGGCCACCCGATCGACGAGt               <  8:148650/149‑1 (MQ=255)
                                                                                                                            tCGACCTGTTCGGCAGCATGTTCGACGCTGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGACACCCTTGGCCTGTACCTCACCGGCCACCCGATCGACGAGTAc             <  4:35583/149‑1 (MQ=255)
                                                                                                                            tCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGACACCCTTGGCCTGTACCTCACCGGCCACCCGATCGACGAGTAc             <  1:305585/149‑1 (MQ=255)
                                                                                                                                       ggCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGACACCCTTGGCCTGTACCTCACCGGCCACCCGATCGACGAGTACGAGACCGAGAt  <  6:7084/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
TTCCATGACGAGATCAAGGCCTACCACGCCAACATCGACATCAACCGGGCAACCTTGCTGTCGGCGCTGGGCGAGGCCATCAAGGCAGCCGAGCAGGCTGCGCATACCGCTGACAGTGGCCATGTCGACCTGTTCGGCAGCATGTTCGACGCGGCAGATGTCGACGTCTATGCCAACCATCGCAAGGTGCGCGAGCTCACCCTCAAGGAGCGTCTGAAGGGGGAGAAGGACACCCTTGGCCTGTACCTCACCGGCCACCCGATCGACGAGTACGAGACCGAGAT  >  NC_002947/1802705‑1802988

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