Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273687,3030AG56.5% 16.8 / 30.6 24intergenic (‑63/+27)corC/ybeYCNNM family magnesium/cobalt transport protein CorC/rRNA maturation RNase YbeY
Reads supporting (aligned to +/- strand):  ref base A (1/9);  new base G (13/0);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 1.22e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-05
Rejected as consensus: Frequency below/above cutoff threshold.

GTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGC  >  NZ_CP009273/687220‑687371
                                                                                   |                                                                    
gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGtgtg                                                                >  1:471633/1‑90 (MQ=255)
gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTccg                                                                >  1:141166/1‑88 (MQ=255)
gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                >  2:625572/1‑90 (MQ=255)
gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                >  1:118381/1‑90 (MQ=255)
gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                >  2:719440/1‑90 (MQ=255)
gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                >  2:37588/1‑90 (MQ=255)
gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                >  2:25665/1‑90 (MQ=255)
gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTg                                                                >  2:724074/1‑90 (MQ=255)
  gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGcc                                                              >  2:452799/1‑90 (MQ=255)
  gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGATGTCTGcc                                                              <  2:39635/90‑1 (MQ=255)
         gtcgtcGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCaaa                                                       <  2:719880/90‑1 (MQ=255)
          tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGCTGCCCGGCaaaa                                                      >  1:44330/1‑89 (MQ=255)
                        gtcgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCCGGCAAAAGCCCGGGAAaatat                                        >  1:19798/1‑90 (MQ=255)
                          cgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGTCAAATGCCTGGCAAattattt                                      >  2:132343/1‑89 (MQ=255)
                          cgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCCGTCAAAAGCCTGGCCAattattt                                      >  1:639098/1‑89 (MQ=255)
                                ggTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTtctc                                <  2:367496/90‑1 (MQ=255)
                                   tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAAAGCCTGGCCAAATATTTTTTTTCGGc                             >  1:40538/1‑90 (MQ=255)
                                   tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGc                             <  1:724074/90‑1 (MQ=255)
                                   tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGc                             <  1:112248/90‑1 (MQ=255)
                                     tcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCaa                           <  2:664009/90‑1 (MQ=255)
                                     tcATGTTAGTTAATCAAAACGCCATCGTTAATCACCAACGGCGGGGCCGTCTCCCAGTCAAAAGCCTGGCAAATTATTCTTTCTCGGCaa                           <  1:428944/90‑1 (MQ=255)
                                            agttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACgg                    <  2:667998/90‑1 (MQ=255)
                                               taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc                 <  1:36388/90‑1 (MQ=255)
                                                              cgtTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGagc  <  1:129821/90‑1 (MQ=255)
                                                                                   |                                                                    
GTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGC  >  NZ_CP009273/687220‑687371

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