Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,189,041 A→G 11.0% intergenic (+25/+17) PP_1038 → / ← tatC‑I hypothetical protein/Sec‑independent protein translocase protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,189,0410AG11.0% 111.5 / 4.6 45intergenic (+25/+17)PP_1038/tatC‑Ihypothetical protein/Sec‑independent protein translocase protein
Reads supporting (aligned to +/- strand):  ref base A (21/19);  new base G (2/3);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

TGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTC  >  NC_002947/1188900‑1189163
                                                                                                                                             |                                                                                                                           
tGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCcgcg                                                                                                                     >  2:346791/1‑150 (MQ=255)
  gAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGACACCCGTACGAGACGCCAGCCTATGAAGTATGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGaa                                                                                                                   <  2:689949/150‑1 (MQ=255)
      gTGGAGCTGGTGGTGGCTGACGAGATGATTGGGCAGTTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAggg                                                                                                              <  2:498619/150‑1 (MQ=255)
                     gCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGcc                                                                                               >  2:659656/1‑150 (MQ=255)
                       tGACGAGCTGTTTGTGCAGCTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTTTGGCGGCTCGCCGAGTACTAGAACGCATCACGGCTCCTCGCGGGCATTACCGCGAAAGGGCCCTACAATCCAGCcac                                                                                             <  1:381289/150‑1 (MQ=255)
                       tGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGaaa                                                                                                                 <  2:274308/130‑1 (MQ=255)
                       tGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGaaa                                                                                                                 >  1:274308/1‑130 (MQ=255)
                       tGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCcac                                                                                             >  2:526341/1‑150 (MQ=255)
                           gAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAAtc                                                                                         >  1:308401/1‑150 (MQ=255)
                               tGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCt                                                                                     >  1:136567/1‑150 (MQ=255)
                                                gCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGc                                                                    <  1:21634/150‑1 (MQ=255)
                                                  tGCGGTCAAGCAGAGCCACGCGCACGAGACGCCAGCCTATGAAGTCTGGCGGTTCGCCGAGTTCTAGAACGCATCACTGCTATTCGGGGGCACGCCCGTGAAAGGGCCCTACATTCCAGCCACAATCTCCTTGGCTGGCTCCGGCCGCtt                                                                  <  1:564831/150‑1 (MQ=255)
                                                  tGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCAAGCAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCtt                                                                  >  1:250647/1‑150 (MQ=255)
                                                   gcgcTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTc                                                                 <  2:75619/150‑1 (MQ=255)
                                                      cTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGa                                                                                                                    <  1:712786/97‑1 (MQ=255)
                                                      cTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGa                                                                                                                    >  2:712786/1‑97 (MQ=255)
                                                      cTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAAc                                                              >  2:131610/1‑150 (MQ=255)
                                                       tCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACa                                                             >  1:714409/1‑150 (MQ=255)
                                                          aGCAGAGCCACCCGTACGAGACGCCAGCATATGAAGTCGGGCTGCTCGCCGAGTTCTAGAACGCATCACGCCAGTTTGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTc                                                          <  1:52277/150‑1 (MQ=255)
                                                               agCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACaa                                                     <  2:136567/150‑1 (MQ=255)
                                                                   aCCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCGCGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCa                                                 <  1:131610/150‑1 (MQ=255)
                                                                    cccGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAg                                                >  1:437412/1‑150 (MQ=255)
                                                                           gagaCGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAACACACGTCTGAACTCCAGTACAATTGCTGTTCCTGGCTCCGGCCGCTTCAACATTACACAAGCCAGCAAcccc                                         >  1:93598/1‑150 (MQ=255)
                                                                                gCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCACtctc                                    >  1:134401/1‑150 (MQ=255)
                                                                                   aGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCCCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGaa                                 <  2:308401/150‑1 (MQ=255)
                                                                                      cTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACac                              >  2:208385/1‑149 (MQ=255)
                                                                                        aTGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCa                            <  2:714409/150‑1 (MQ=255)
                                                                                             gTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAAc                                            <  2:116969/129‑1 (MQ=255)
                                                                                             gTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAAc                                            >  1:116969/1‑129 (MQ=255)
                                                                                                  gcggcTCGCCGAGTTCTAGAACGCATCACGGGTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCgg                  <  1:100627/150‑1 (MQ=255)
                                                                                                  gcggcTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCgg                  >  1:684250/1‑150 (MQ=255)
                                                                                                  gcggcTCGCCGAGTTCTAGAACGCATCACGGCTATTCACGCGCACGCCCGAGAAAAGACACAACAATCCAACCACAATCTCCTTTCCTGTCTCCGGCCTCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCgg                  >  2:290981/1‑150 (MQ=255)
                                                                                                  gcggcTCGCCGAGTCCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGGCCTAGAATCCGGCCACAATCTGCTTGCCTGGGTCCGGCCGCTGCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCgg                  >  2:291160/1‑150 (MQ=255)
                                                                                                        cGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGc            <  1:346791/150‑1 (MQ=255)
                                                                                                         gCCGAGTTCGAGACCGCATCACGGTTTTTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCTCTTCAACTGCTGACCAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGcc           <  1:15522/150‑1 (MQ=255)
                                                                                                               ttCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACt              <  2:742475/141‑1 (MQ=255)
                                                                                                               ttCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACt              >  1:742475/1‑141 (MQ=255)
                                                                                                               ttCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATc     <  1:659656/150‑1 (MQ=255)
                                                                                                               ttCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATc     >  2:509818/1‑150 (MQ=255)
                                                                                                                 cTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGt   <  1:526341/150‑1 (MQ=255)
                                                                                                                  tAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTc  >  1:555929/1‑150 (MQ=255)
                                                                                                                             aCGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCAccc                         <  2:783743/116‑1 (MQ=255)
                                                                                                                             aCGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCAccc                         >  1:783743/1‑116 (MQ=255)
                                                                                                                               gccctttCGCGGGCGTGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGaa                                 <  2:104154/101‑1 (MQ=255)
                                                                                                                               gccctttCGCGGGCGTGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGaa                                 >  1:104154/6‑106 (MQ=255)
                                                                                                                                             |                                                                                                                           
TGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTC  >  NC_002947/1188900‑1189163

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