Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273687,3030AG52.6% 7.8 / 20.9 19intergenic (‑63/+27)corC/ybeYCNNM family magnesium/cobalt transport protein CorC/rRNA maturation RNase YbeY
Reads supporting (aligned to +/- strand):  ref base A (2/7);  new base G (10/0);  total (12/7)
Fisher's exact test for biased strand distribution p-value = 7.14e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGA  >  NZ_CP009273/687220‑687369
                                                                                   |                                                                  
gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                              >  2:776816/1‑90 (MQ=255)
gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                              >  1:611296/1‑90 (MQ=255)
gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                              >  2:51167/1‑90 (MQ=255)
                   aTGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGCTGCCCGTCAAAAGCCTGGGaa                                           >  2:867989/1‑90 (MQ=255)
                                   tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGTCAAATGCCTgg                                              >  1:494153/1‑71 (MQ=255)
                                   tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTgg                                              <  2:494153/71‑1 (MQ=255)
                                     tcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCCGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCaa                         <  2:315202/90‑1 (MQ=255)
                                         gttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCCGTCAAAAGCCTGGCAAATTATTTTTTTTCGGGAAAGTa                     >  2:856149/1‑90 (MQ=255)
                                          ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTAc                    <  1:236807/90‑1 (MQ=255)
                                          ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTAc                    <  1:856149/90‑1 (MQ=255)
                                           tagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGTCAAAAGCCTGGCAAAATATTTTTTTTCGGCAATGTACg                   >  2:236807/1‑90 (MQ=255)
                                            agttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACgg                  <  1:702227/90‑1 (MQ=255)
                                                aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCCGTCAAAAGCCTGGCAAAATATTTTTTTTCCGCCAAGTACGGATcc              >  2:175875/1‑90 (MQ=255)
                                                aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAATGGCCGGCAAATTATTTTTTTTCGGGCATGTACGGAGcc              >  1:633312/1‑90 (MQ=255)
                                                aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAAAGCCTGGCAAATTATTCTTTTTCGGCAATGTACGGATcc              >  2:702227/1‑90 (MQ=255)
                                                    aaaaCGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCAt          <  1:175875/90‑1 (MQ=255)
                                                     aaaCGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATa         >  1:440819/1‑90 (MQ=255)
                                                          ccGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGGCAAATGCCCGGCAAATTATTTTTTTTCGGCAATGTACGGGTCCTCATAGCCCa    >  1:468416/1‑90 (MQ=255)
                                                            gtcgtTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCaga  <  2:440819/90‑1 (MQ=255)
                                                                                   |                                                                  
GTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGA  >  NZ_CP009273/687220‑687369

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