Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,269,870 C→G 16.7% intergenic (‑526/‑35) PP_3740 ← / → mrdA‑I MFS transporter/transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,269,8700CG16.7% 55.3 / 6.5 24intergenic (‑526/‑35)PP_3740/mrdA‑IMFS transporter/transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (12/8);  new base G (2/2);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

ATGACTCCGCAATGAACAAATCTTCATATTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGCCTGTGTGCCGTCTTGGTGGG  >  NC_002947/4269732‑4270014
                                                                                                                                          |                                                                                                                                                 
atGACTCCGCAATGAACAAATCTTCATATTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATTGc                                                                                                                                       >  1:226423/1‑150 (MQ=255)
      ccGAAATGAACAAATCTGCATATTCCCCCCCTCCGCCCTGATCATGGGGCACACAGCAGATGCCAGGCAGCGGGCAGCCGCGGCTCGTCTTGACGGTGCGGTACGGTTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATGGCCAcccc                                                                                                                                 <  2:129416/150‑1 (MQ=255)
                  aaTCTTCATATTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATAGCCACCCCTCgggtgcgact                                                                                                                     <  2:596642/150‑8 (MQ=255)
                          tatTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATg                                                                                                             >  1:209795/1‑150 (MQ=255)
                           atTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTcaca                                                                                                                 >  1:378344/1‑145 (MQ=255)
                           atTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGAGCGCGCGCGTATTGCCACCCCTCAGGATCGATTcaca                                                                                                                 <  2:378344/145‑1 (MQ=255)
                                                      gCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACg                                                                                 >  1:140576/1‑150 (MQ=255)
                                                                     cagccagcGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCAcc                                                                  <  2:356406/150‑1 (MQ=255)
                                                                      agccagcGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAAc                                                                                                      <  1:226269/113‑1 (MQ=255)
                                                                      agccagcGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAAc                                                                                                      >  2:226269/1‑113 (MQ=255)
                                                                      agccagcGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCt                                                                 >  2:367965/1‑150 (MQ=255)
                                                                       gccagcGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTg                                                                <  2:209795/150‑1 (MQ=255)
                                                                                                        cGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCCCGAAAAGGAAATGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATg                               <  2:226423/150‑1 (MQ=255)
                                                                                                          atcctgaggggtggcaataCGCGCGCGATCGGGCCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTaaa                                                                                         >  1:65847‑M1/20‑90 (MQ=255)
                                                                                                          atcctgaggggtggcaataCGCGCGCGATCGGGCCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTaaa                                                                                         <  2:65847‑M1/71‑1 (MQ=255)
                                                                                                              gggTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGAACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTc                         >  2:69031/1‑150 (MQ=255)
                                                                                                                 ttGCGTAGCATCCGCGGCCCGATCGCGCGCGTCTTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGc                      >  2:227907/1‑150 (MQ=255)
                                                                                                                 ttGCGTAGCATCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGc                      >  2:405089/1‑150 (MQ=255)
                                                                                                                     gTAGCATCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGCCtgt                  >  1:212214/1‑150 (MQ=255)
                                                                                                                         aataCGCGCGCGATCGGGCCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACgaaaa                                                                             >  1:83640‑M1/5‑87 (MQ=255)
                                                                                                                         aataCGCGCGCGATCGGGCCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACgaaaa                                                                             <  2:83640‑M1/83‑1 (MQ=255)
                                                                                                                          aTCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGCCTGTGTGcc             >  2:279154/1‑150 (MQ=255)
                                                                                                                                GGCCCGATCGCGCGCGTATTGCCACCCCTGAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGCCTGTGTGCCGTCTTg       >  2:91736‑M1/1‑150 (MQ=255)
                                                                                                                                     gATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGCCTGTGTGCCGTCTTGGTggg  <  1:405089/150‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
ATGACTCCGCAATGAACAAATCTTCATATTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGCCTGTGTGCCGTCTTGGTGGG  >  NC_002947/4269732‑4270014

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