Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,898,571 A→T 12.2% L300I (TTA→ATA)  PP_4307 ← DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,898,5710AT12.2% 121.3 / 6.3 49L300I (TTA→ATA) PP_4307DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base A (24/19);  new base T (2/4);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAATTGCCGCAATC  >  NC_002947/4898427‑4898719
                                                                                                                                                |                                                                                                                                                    
tccGTGCAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGGAGGTGTTGGTCAAGGGTGTCGATGGGCGATGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAAccc                                                                                                                                                  <  7:360441/148‑1 (MQ=255)
   gTGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCcgcg                                                                                                                                               >  6:203458/1‑149 (MQ=255)
     ggAACTGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTTTACCcgcgaa                                                                                                                                             <  4:226447/149‑3 (MQ=255)
      gAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGCGCGGCAAGCAGGTGTTGCTCAAGAGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTTTACCcgcgaag                                                                                                                                            <  5:365276/149‑4 (MQ=255)
       aaCAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGAGCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTTTACCcgcgaaga                                                                                                                                           <  3:231464/149‑5 (MQ=255)
            gccgcggGGTTCCTCGTGGGCAGGGTCGAGGGCGGCAAGCCGCGGGTGCTCCCGGGTGTCGATGGGCGAGGCCCTGGCGGATCGCAGGCTGGCAAACAGCGGCCTTCTCCCTGGCGGGCCGCTTCGGGGGTAAACCCGCGCCTACAGGa                                                                                                                                      <  4:237396/146‑1 (MQ=255)
                    tttCCTGTGGGGCAGGGTCGAGTGCGGCAAGCAGGTGTTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCgcgc                                                                                                                              <  6:222389/149‑1 (MQ=255)
                           tgGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTGTCGATTGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCTCGCCTACAGGATGCCGCGCCTAGc                                                                                                                         <  6:35007/147‑1 (MQ=255)
                           tgGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTGTCGATGGGCGATGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGc                                                                                                                         <  6:371712/147‑1 (MQ=255)
                           tgGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGc                                                                                                                         <  7:277535/147‑1 (MQ=255)
                            ggggCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCAGCGCCTACAGGATGCCGCGCCGAGCgg                                                                                                                       >  7:202244/1‑148 (MQ=255)
                                cAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAacca                                                                                                                  <  7:109121/149‑1 (MQ=255)
                                   ggTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCTGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCAGAGCGGAacaaccc                                                                                                               >  5:321612/1‑148 (MQ=255)
                                           gcggcAAGCAGGTGTTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCctgct                                                                                                       >  4:38663/1‑149 (MQ=255)
                                            cggcAAGCAGGTGTTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCctgctg                                                                                                      >  2:148662/1‑149 (MQ=255)
                                               cAAGCAGGTGTTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCAt                                                                                                   <  5:203458/149‑1 (MQ=255)
                                                   cAGGTGTTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGccc                                                                                               <  2:210682/149‑1 (MQ=255)
                                                        gtTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTgc                                                                                          >  7:267580/1‑149 (MQ=255)
                                                         ttGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTgcc                                                                                         >  8:197052/1‑149 (MQ=255)
                                                             tCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTgccgc                                                                                       >  1:187880/1‑147 (MQ=255)
                                                                         attggCGAGGTGATGGCTGATCGAAGTCTGGCAATAAGCGGCTTTCTTCCTTGCGTGCCTCTTCGCTGGTAAACCCGCGCCTACTGGATGCCGCGCCGAGCGGAACCACAACCCCTGCTGGATGCCCGCTGCCGCGATCAGGATCAacc                                                                         <  8:202244/146‑1 (MQ=255)
                                                                             gCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCGCCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCGaccccgc                                                                     >  8:128609/1‑147 (MQ=255)
                                                                                         cGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAaccccac                                                                     <  8:318503/137‑1 (MQ=255)
                                                                                         cGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAaccccac                                                                     >  7:318503/1‑137 (MQ=255)
                                                                                                cAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAgcggc                                                  <  1:148662/149‑1 (MQ=255)
                                                                                                 aGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCCCCACCCCTGCTGCATGCACGCTGCCGCCAGCAGCATCAACCCCACCCCCAGCCACACTAgcggcg                                                 >  1:337425/1‑149 (MQ=255)
                                                                                                   gCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACt                                                          >  8:116696/1‑138 (MQ=255)
                                                                                                   gCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACt                                                          <  7:116696/138‑1 (MQ=255)
                                                                                                     tGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTATCCCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCAGGCCCGCTGCCGCCGGCAGGGTCGGCCCGGCCCCCAGCCACTGTAGCGGCGCCAg                                             >  6:13938/1‑149 (MQ=255)
                                                                                                     tGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAg                                             >  5:358856/1‑149 (MQ=255)
                                                                                                      ggCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAgg                                            >  7:99702/1‑149 (MQ=255)
                                                                                                         aaaaaGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCCCCCCCAGCCACTGTAGCGGCGCCAGGCgg                                         >  7:210479/1‑149 (MQ=255)
                                                                                                           aaaGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCcg                                                                                                                                 >  4:356369/1‑59 (MQ=255)
                                                                                                           aaaGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCcg                                                                                                                                 <  3:356369/59‑1 (MQ=255)
                                                                                                                ggCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGa                                  >  1:214540/1‑149 (MQ=255)
                                                                                                                 gccttCTCCCTTGCGGGCGTCTTCGCGGGTAACCCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCGGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCTGCGCCAGGCGGTGCCCGaa                                 <  6:321612/146‑1 (MQ=255)
                                                                                                                  cATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAAc                                <  6:358856/149‑1 (MQ=255)
                                                                                                                       tACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCAc                           >  4:96745/1‑149 (MQ=255)
                                                                                                                         ccTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCAcc                          <  2:187880/148‑1 (MQ=255)
                                                                                                                          cTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCAcc                          >  2:316550/1‑147 (MQ=255)
                                                                                                                             gCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGt                      <  3:38663/148‑1 (MQ=255)
                                                                                                                                ggCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTccac                  >  1:352889/1‑149 (MQ=255)
                                                                                                                                   ctctTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTccac                  >  4:130276/1‑146 (MQ=255)
                                                                                                                                       gcgcgGGTTTACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAAtt           >  4:155697/2‑149 (MQ=255)
                                                                                                                                         gcgGGTTTACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAATTGc         <  3:155697/149‑1 (MQ=255)
                                                                                                                                         gcgGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAATTGc         <  8:267580/149‑1 (MQ=255)
                                                                                                                                          cgGGTAAACCCGCGCCTATAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAATTGcc        >  4:62896/1‑149 (MQ=255)
                                                                                                                                            ggTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAATTGCCg       >  6:185035/1‑148 (MQ=255)
                                                                                                                                            ggTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAATTGCCg       >  3:265431/1‑148 (MQ=255)
                                                                                                                                                aaCCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGCTCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAATTGCCGCAATc  >  3:298212/1‑149 (MQ=255)
                                                                                                                                                aaCCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAATTGCCGCAATc  >  3:293105/1‑149 (MQ=255)
                                                                                                                                                aaCCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAATTGCCGCAATc  <  3:130276/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ACCGTGGAACAGGCGGCGGGTTTCCTGTGGGGCAGGGTCGAGGGCGGCAAGCAGGTGTTGCTCAAGGGTGTCGATGGGCGAGGTCATGGCGGATCGCAGGCTGGCAAAAAGCGGCATTCTACCTTGCGGGCCTCTTCGCGGGTAAACCCGCGCCTACAGGATGCCGCGCCGAGCGGAACCACCACCCCTGCTGCATGCCCGCTGCCGCCAGCAGGATCAACCCCACCCCCAGCCACTGTAGCGGCGCCAGGCGGTGCCCGAACGCCACCCAGTCCACCAGAATTGCCGCAATC  >  NC_002947/4898427‑4898719

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