Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273687,3030AG47.1% 30.2 / 10.7 34intergenic (‑63/+27)corC/ybeYCNNM family magnesium/cobalt transport protein CorC/rRNA maturation RNase YbeY
Reads supporting (aligned to +/- strand):  ref base A (8/10);  new base G (15/1);  total (23/11)
Fisher's exact test for biased strand distribution p-value = 3.08e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAAT  >  NZ_CP009273/687219‑687380
                                                                                    |                                                                             
tgtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCt                                                                          >  2:558346/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGTCTg                                                                         >  2:432567/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                         >  1:270411/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                         >  2:242292/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg                                                                         >  2:512958/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                         >  2:287256/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg                                                                         >  1:577157/1‑90 (MQ=255)
 gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTg                                                                         >  2:84766/1‑90 (MQ=255)
  tgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGc                                                                        >  2:68686/1‑90 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGcc                                                                       >  1:224598/1‑90 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGcc                                                                       >  2:470151/1‑90 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGcc                                                                       >  2:217512/1‑90 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGcc                                                                       >  1:605967/1‑90 (MQ=255)
   gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGcc                                                                       <  2:74863/90‑1 (MQ=255)
                  tCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGc                                                        <  1:84766/90‑1 (MQ=255)
                    aTGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGTCTGCCCGTCAAAAACCTGGCaa                                                      >  1:20205/1‑90 (MQ=255)
                    aTGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGGTGCCCGGCAAAAGCCTGGcca                                                      >  1:747467/1‑88 (MQ=255)
                           cgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTAttc                                               >  2:95583/1‑90 (MQ=255)
                                    tctcATGTTATTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGGAAATTATTTTTTTTCGGc                                      >  2:33395/1‑90 (MQ=255)
                                      tcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGCTGCCAGGCAAAAGCCCGGCCAATTATTGTTTCTCCGCaa                                    >  2:859644/1‑90 (MQ=255)
                                           ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGTCAAAAGCCTGGCAAATTTTTTTTTCTCCGCCATGTAc                               >  2:48905/1‑90 (MQ=255)
                                           ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGc                                      >  2:109709/1‑83 (MQ=255)
                                           ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGc                                      <  1:109709/83‑1 (MQ=255)
                                            tagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGGCAAAAGCCTGGCAAATTATTCTTTTTCGGCAATGTACg                              >  1:280099/1‑90 (MQ=255)
                                            tagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTTTCGGCAATGTACg                              >  2:223814/1‑90 (MQ=255)
                                                taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc                          <  2:605967/90‑1 (MQ=255)
                                                taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc                          <  2:747467/90‑1 (MQ=255)
                                                taAACAAAACGCGGCGGTTAATCACCAACGACGGGGGCGTCTGCCAGTCAAATGCCTGGCAAACTATTCTTCGTCGGCAATGTACGGAtt                          <  1:660497/90‑2 (MQ=255)
                                                     aaaaCGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCAt                     <  1:89381/90‑1 (MQ=255)
                                                          gCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGccc                <  1:68686/90‑1 (MQ=255)
                                                           ccGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCa               <  2:291419/90‑1 (MQ=255)
                                                               cgtTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGagc           <  2:329581/90‑1 (MQ=255)
                                                                      cACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATa    <  1:95583/90‑1 (MQ=255)
                                                                        ccAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTTTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAAt  >  2:254579/1‑90 (MQ=255)
                                                                                    |                                                                             
TGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAAT  >  NZ_CP009273/687219‑687380

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