Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,899,298 T→G 7.6% I72I (ATA→ATC)  PP_4307 ← DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,899,2980TG 7.6% 150.5 / 2.6 56I72I (ATA→ATC) PP_4307DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base T (22/29);  new base G (2/3);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCTTTTCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAACCAGC  >  NC_002947/4899154‑4899446
                                                                                                                                                |                                                                                                                                                    
gATCTTTTCGCCCAGAAACAGCGCCGCCAGCCCCCCCAGCATGAACGGATGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGc                                                                                                                                                 <  1:1036194/150‑1 (MQ=255)
  tCTTTTCGCCCAGAAACAGCGCCGCCAGCCCAAACCGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCa                                                                                                                                               <  2:621991/150‑1 (MQ=255)
     tttCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATTAACGGCTGCACGTTGTACACCGCCGTGCCGATGGAAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCGATAGCCAGcg                                                                                                                                            <  1:1119101/150‑2 (MQ=255)
       tCGCCCAGAAACAGCGCCGCCAGCCACACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCGATAGCCAGcggg                                                                                                                                          <  1:62267/150‑4 (MQ=255)
          cccAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAAt                                                                                                                                       >  1:1240159/1‑150 (MQ=255)
            cAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGc                                                                                                                                     >  1:1083532/1‑150 (MQ=255)
             aGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTTGCTATCGCCAGCAAAAATGcc                                                                                                                                    <  2:386114/150‑1 (MQ=255)
             aGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGcc                                                                                                                                    <  2:866664/150‑1 (MQ=255)
                 aCAGCGCCGCCAGCCGCACCAGCATGAACGGCTGCACGTTGTACACACCCGTTCCGCTGGGAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGc                                                                                                                                >  2:586162/1‑150 (MQ=255)
                 aCAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTCCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGc                                                                                                                                >  1:1125069/1‑150 (MQ=255)
                   aGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCgg                                                                                                                              >  2:1224698/1‑150 (MQ=255)
                    gcgcCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGc                                                                                                                             >  1:740910/1‑150 (MQ=255)
                       ccgccAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGa                                                                                                                          <  2:552641/150‑1 (MQ=255)
                             gCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCAccc                                                                                                                    <  1:78376/150‑1 (MQ=255)
                              ccccACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCAcccc                                                                                                                   <  2:806451/150‑1 (MQ=255)
                              ccccACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCAcccc                                                                                                                   >  2:149287/1‑150 (MQ=255)
                               cccACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCGATAGCCAGcgg                                                                                                                                           >  1:817953/1‑123 (MQ=255)
                               cccACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCGATAGCCAGcgg                                                                                                                                           <  2:817953/125‑3 (MQ=255)
                                 caccaGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCgg                                                                                                                <  1:1224698/150‑1 (MQ=255)
                                 caccaGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCgg                                                                                                                >  2:760810/1‑150 (MQ=255)
                                      gCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCa                                                                                                           <  1:148684/150‑1 (MQ=255)
                                          gAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGaaa                                                                                                       >  2:491767/1‑150 (MQ=255)
                                          gAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGaaa                                                                                                       <  1:406255/150‑1 (MQ=255)
                                                  gCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGc                                                                                                     >  1:1203739/1‑144 (MQ=255)
                                                  gCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGc                                                                                                     <  2:1203739/144‑1 (MQ=255)
                                                             acCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGa                                                                                                                          >  2:1006650/1‑112 (MQ=255)
                                                             acCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGa                                                                                                                          <  1:1006650/112‑1 (MQ=255)
                                                                      ccGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCaaa                                                                                      <  1:1016615/139‑1 (MQ=255)
                                                                      ccGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCaaa                                                                                      >  2:1016615/1‑139 (MQ=255)
                                                                      ccGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAg                                                                           >  1:117314/1‑150 (MQ=255)
                                                                      ccGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAg                                                                           <  2:1083532/150‑1 (MQ=255)
                                                                       cGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAgg                                                                          >  2:550871/1‑150 (MQ=255)
                                                                           ggCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCAc                                                                      <  2:345581/150‑1 (MQ=255)
                                                                                tCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCAc                                                                 <  2:117314/150‑1 (MQ=255)
                                                                                    cGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGaa                                                             <  1:149287/150‑1 (MQ=255)
                                                                                              tATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGc                                                                   >  1:443404/1‑134 (MQ=255)
                                                                                              tATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGc                                                                   <  2:443404/134‑1 (MQ=255)
                                                                                                tGAAGCGAACCACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCTTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCa                                                 <  1:859391/150‑1 (MQ=255)
                                                                                                tGAAGCGAACAACAGCAGCCAGTTGCCGACCATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCa                                                 <  1:586162/150‑1 (MQ=255)
                                                                                                    gCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCAc                                                                 <  1:100339/130‑1 (MQ=255)
                                                                                                    gCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCAc                                                                 >  2:100339/1‑130 (MQ=255)
                                                                                                    gCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCCCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAac                                             <  2:740910/150‑1 (MQ=255)
                                                                                                      gAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAacac                                           >  1:638224/1‑150 (MQ=255)
                                                                                                               cagcCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAg                                  <  2:539718/150‑1 (MQ=255)
                                                                                                                          cGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGc                       >  2:580898/1‑150 (MQ=255)
                                                                                                                           gACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACgg                                                                    >  2:282591/1‑104 (MQ=255)
                                                                                                                           gACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACgg                                                                    <  1:282591/104‑1 (MQ=255)
                                                                                                                           gACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGcc                      >  2:1054860/1‑150 (MQ=255)
                                                                                                                            aCAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTAACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCg                     <  1:550871/150‑1 (MQ=255)
                                                                                                                                 ggCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGacac                <  2:1125069/150‑1 (MQ=255)
                                                                                                                                     aCCCCGCTGGCGATCGCCAGCAAAAATGCCGTGCGGATGATCATTTCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGACTACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGc            >  1:1125362/1‑150 (MQ=255)
                                                                                                                                      ccccGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCa           >  2:18632/1‑150 (MQ=255)
                                                                                                                                          gCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATTCGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGaa       <  1:580898/150‑1 (MQ=255)
                                                                                                                                              gCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAACCa    >  2:1213935/1‑149 (MQ=255)
                                                                                                                                              gCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAACCa    <  1:1213935/149‑1 (MQ=255)
                                                                                                                                               cTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAACCAGc  <  1:760810/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GATCTTTTCGCCCAGAAACAGCGCCGCCAGCCCCACCAGCATGAACGGCTGCACGTTGTACACCGCCGTGCCGATGGCAATCGACGCCCGGGCGTATGAAGCGAACAACAGCAGCCAGTTGCCGACAATGGCCACCCCGCTGGCTATCGCCAGCAAAAATGCCGTGCGGCTGATCACCCCCGGGCGCAGAAAGCCAAACGCGGCGCAAATCGCCAGCAAGGTCACGGCACCGAATACGCAGCGCCAGAACACCACATCCAGCACCGGCTGGCCGGACACCAGCACGAACCAGC  >  NC_002947/4899154‑4899446

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