Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,359,212 A→G 31.0% intergenic (‑74/+8) PP_1183 ← / ← PP_1184 siderophore biosynthesis protein/dienelactone hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,359,2120AG31.0% 6.3 / 5.7 13intergenic (‑74/+8)PP_1183/PP_1184siderophore biosynthesis protein/dienelactone hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base A (3/6);  new base G (2/2);  total (5/8)
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.35e-01

CCGGGCAATGGACGGGGCAGGGGCCAGTGGTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTTATAGCCAATGTCCGGCCCCCCATGTTCACCGTGGCTCAGCCGATCGGCATCGGGGTTGGTGAACCCGTGCTTGGCC  >  NC_002947/1359071‑1359361
                                                                                                                                             |                                                                                                                                                     
ccGGGCAATGGACGGGGCAGGGGCCAGTGGTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCttt                                                                                                                                               <  2:173510/150‑1 (MQ=255)
 cGGGCAATGGACGGGGCAGGGGCCAGTGGTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCtttt                                                                                                                                              <  1:234015/150‑1 (MQ=255)
  gggCAATGGACGGGGCAGGGGCCAGTGGTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTAGTCAAGCCTACTTTACCATAACGATTGCACAATCTTTTa                                                                                                                                             >  2:96410/1‑150 (MQ=255)
                                         ggCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCa                                                                                                      <  1:96410/150‑1 (MQ=255)
                                                                     gcggcTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTTATAg                                                                          <  2:36678/150‑1 (MQ=255)
                                                                               tgCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTTATAGCCAATGTCCg                                                                >  2:1337/1‑150 (MQ=255)
                                                                                          cAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCGCAATCTTTTATTTGAACACCTTCTTg                                                                                                                             <  1:390250/78‑1 (MQ=255)
                                                                                          cAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCGCAATCTTTTATTTGAACACCTTCTTg                                                                                                                             <  1:390343/78‑1 (MQ=255)
                                                                                          cAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCGCAATCTTTTATTTGAACACCTTCTTg                                                                                                                             >  2:390250/1‑78 (MQ=255)
                                                                                          cAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCGCAATCTTTTATTTGAACACCTTCTTg                                                                                                                             >  2:390343/1‑78 (MQ=255)
                                                                                                     ggAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTTATAGCCAATGTCCGGCCCCCCATGTTCACCGTGGCt                                          >  2:292751/1‑150 (MQ=255)
                                                                                                                                ccACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTTATAGCCAATGTCCGGCCCCCCATGTTCACCGTGGCTCAGCCGATCGGCATCGGGGTTGGTGaa               <  1:292751/150‑1 (MQ=255)
                                                                                                                                          tGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGACTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTTATAGCCAATGTCCGGCCCCCCATGTTCACCGTGGCTCAGCCGATCGGCATCGGGGTTGGTGAACCCGTGCTTg     <  1:315795/150‑1 (MQ=255)
                                                                                                                                             acaATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTTATAGCCAATGTCCGGCCCCCCATGTTCACCGTGGCTCAGCCGATCGGCATCGGGGTTGGTGAACCCGTGCTTGGcc  >  2:219662/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
CCGGGCAATGGACGGGGCAGGGGCCAGTGGTGCCGGAGCGGGGCGCAGCAGGCGGGAAGTGTGTTCATGGCGGCTATTGTGCCAGCTGTGCAGGCGCTAAGGGAGGTGGCGTGGTAACGCCTCCTTTCCCACAACGATTGCACAATCTTTTATTTGAACACCTTCTTGAAGAACGCCTGCATATCTGCCCAGGAACGCTCGTCCGCTGCCTTGTTATAGCCAATGTCCGGCCCCCCATGTTCACCGTGGCTCAGCCGATCGGCATCGGGGTTGGTGAACCCGTGCTTGGCC  >  NC_002947/1359071‑1359361

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