Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,762,580:1 +G 23.5% coding (1885/2091 nt) PP_3325 → ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,762,5801.G23.5% 32.4 / 13.5 17coding (1885/2091 nt)PP_3325ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base . (6/7);  new base G (2/2);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

ACCAAAGGCGACGCAGCCCTGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTCCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCACTACAAGGCCAGCGAACAGGTCACCGTCGGCCTTAACCTCAACAACCTGTTCGACCGCAAGTACTA  >  NC_002947/3762433‑3762707
                                                                                                                                                    |                                                                                                                               
aCCAAAGGCGACGCAGCCCTGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTcc                                                                                                                               >  1:212206/1‑150 (MQ=255)
         gACGCAGCCCTGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTGCCGGCCTGGc                                                                                                                       >  2:73660/1‑150 (MQ=255)
         gACGCAGCCCTGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTCCGGCCAGGcc                                                                                                                      >  2:430534/1‑150 (MQ=255)
                  cTGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTCCGGCCAGGCCGGCAGCGAc                                                                                                             <  1:45263/150‑1 (MQ=255)
                  cTGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTCCGGCCAGGCCGGCAGCCAc                                                                                                             <  1:508575/150‑1 (MQ=255)
                   tGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTGCCGGCCTGGCCggagcgttc                                                                                                             <  1:73660/150‑8 (MQ=255)
                      ccAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTCCGGCCAGGCCGGCAGCGACTTCa                                                                                                         >  1:229068/1‑150 (MQ=255)
                         aaGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTCCGGCCAGGCCGGCAGCGACTTCAccc                                                                                                      >  2:316799/1‑150 (MQ=255)
                           ggCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTCCGGCCAGGCCGGCAGCGACTTCACCCTg                                                                                                    <  1:430534/150‑1 (MQ=255)
                                                 cccGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTCCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGAcc                                                                              >  1:351113/1‑150 (MQ=255)
                                                           cGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTCCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCa                                                                    <  2:190905/150‑1 (MQ=255)
                                                                     aTGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTGCCGGCCTGGCCggagc                                                                                                                 >  1:306148/1‑93 (MQ=255)
                                                                     aTGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTGCCGGCCTGGCCggagc                                                                                                                 <  2:306148/96‑4 (MQ=255)
                                                                            tgtgCGAAGTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGTGCAGCGGTGAACTACGTCGGCGAACGCTCCGGCCAGGCCGGCAGCGACTTCGCCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCACTACAAGGCCAGCGAAc                                                   <  1:455860/150‑1 (MQ=255)
                                                                            tgtACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTCCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCACTACa                                                               <  1:332475/138‑1 (MQ=255)
                                                                            tgtACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTCCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCACTACa                                                               >  2:332475/1‑138 (MQ=255)
                                                                                                  ccTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTCCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCACTACAAGGCCAGCGAACAGGTCACCGTCGGCCTTAACCt                             <  2:212206/150‑1 (MQ=255)
                                                                                                                             ggTGAACTACGTCGGCGAACGCTCCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCACTACAAGGCCAGCGAACAGGTCACCGTCGGCCTTAACCTCAACAACCTGTTCGACCGCAAGtacta  >  1:284036/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                               
ACCAAAGGCGACGCAGCCCTGCCCAAAGGCAGCGACTTGCTCGGCATCGCCCGTAACAGCGGCAGCCTGATGGGCGTGTACGAATTCCAGGACGGGGCCCTGCGCGGCTCGGATGTGGGCGCAGCGGTGAACTACGTCGGCGAACGCTCCGGCCAGGCCGGCAGCGACTTCACCCTGGACGCCTACAGCACGGTCGACCTGCTGGCGCACTACAAGGCCAGCGAACAGGTCACCGTCGGCCTTAACCTCAACAACCTGTTCGACCGCAAGTACTA  >  NC_002947/3762433‑3762707

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