Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273687,3030AG54.5% 20.6 / 43.0 33intergenic (‑63/+27)corC/ybeYCNNM family magnesium/cobalt transport protein CorC/rRNA maturation RNase YbeY
Reads supporting (aligned to +/- strand):  ref base A (3/12);  new base G (18/0);  total (21/12)
Fisher's exact test for biased strand distribution p-value = 1.28e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-06
Rejected as consensus: Frequency below/above cutoff threshold.

ACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGA  >  NZ_CP009273/687217‑687391
                                                                                      |                                                                                        
actGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCgg                                                                                       >  1:12543/1‑89 (MQ=255)
 ctGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCggc                                                                                      >  1:163738/1‑88 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGTCTg                                                                                    >  2:58002/1‑90 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                                    >  2:357573/1‑90 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                                    >  1:1112746/1‑90 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                                    >  2:1088521/1‑90 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                                    >  1:529560/1‑90 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                                    >  1:210202/1‑90 (MQ=255)
   gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTg                                                                                    >  1:159169/1‑90 (MQ=255)
     gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGcc                                                                                  <  1:357573/90‑1 (MQ=255)
        aaTTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGCTGcccgg                                                                               >  2:1000756/1‑87 (MQ=255)
        aaTTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTTccgcg                                                                               >  2:875780/1‑87 (MQ=255)
             tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGCTGCCAGTCaaaa                                                                          >  2:27918/1‑89 (MQ=255)
             tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCCGTCAAAt                                                                          >  1:1119482/1‑90 (MQ=255)
             tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGGCAAAt                                                                          >  2:452009/1‑90 (MQ=255)
                       tGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGACAGTCAAATGCCTGGCaaa                                                                <  1:114887/90‑1 (MQ=255)
                                 tGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCCCTCAAAAGCCTGGCCAATTAttttttt                                                      >  2:951132/1‑89 (MQ=255)
                                      tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAAAGCCTGGGAAAATATTTTTTTTTGGc                                                 >  2:103877/1‑90 (MQ=255)
                                      tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGc                                                 <  1:791765/90‑1 (MQ=255)
                                      tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGGCTGCCCGTCAAAAGCCTGGCAAAATATTCTTTCTCGGc                                                 >  1:1021520/1‑90 (MQ=255)
                                            gttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCCGTCAAAAGCCTGGCAAAATATTTTTTTTCCGCAATGTa                                           >  2:453400/1‑90 (MQ=255)
                                             ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTAc                                          <  2:159169/90‑1 (MQ=255)
                                               agttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGTCAAATGCCTGGCAAATTATTTTTTTTCGGCAATGTACgg                                        >  1:635130/1‑90 (MQ=255)
                                               agttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGTCAAAAGCCTGGCAAAATATTTTTTCTCGGCAATGTACgg                                        >  1:890683/1‑90 (MQ=255)
                                                  taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc                                     <  2:434735/90‑1 (MQ=255)
                                                  taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc                                     <  2:43880/90‑1 (MQ=255)
                                                  taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc                                     <  1:504710/90‑1 (MQ=255)
                                                      cAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCa                                 <  2:1021520/90‑1 (MQ=255)
                                                            gCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGccc                           <  1:1000756/90‑1 (MQ=255)
                                                               gtcgtTAATCACCAACGGCGGGGACGTCTGCCCGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCAATGTACGGATCCTCATAGCCCaga                        >  1:769315/1‑90 (MQ=255)
                                                                 cgtTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGagc                      <  1:103877/90‑1 (MQ=255)
                                                                            aaCGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAAtct           <  1:1039561/90‑1 (MQ=255)
                                                                                     gACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGa  <  1:843863/90‑1 (MQ=255)
                                                                                      |                                                                                        
ACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGA  >  NZ_CP009273/687217‑687391

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