Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,018,3730CG52.9% 14.2 / 24.0 17A326P (GCA→CCA) yccSYccS family putative transporter
Reads supporting (aligned to +/- strand):  ref base C (1/7);  new base G (9/0);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 4.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAG  >  NZ_CP009273/1018286‑1018450
                                                                                       |                                                                             
gTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTgcgc                                                                             <  1:1387629/90‑1 (MQ=255)
gTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTgcgc                                                                             <  1:680112/90‑1 (MQ=255)
gTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTgcgc                                                                             >  1:574031/1‑90 (MQ=255)
     gATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGttt                                                                        >  1:96398/1‑89 (MQ=255)
     gATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGGGGGCCGttt                                                                        >  1:732582/1‑89 (MQ=255)
           aTTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGGGGGCCGTTTTTgcgc                                                                  >  1:132561/1‑90 (MQ=255)
               ttgCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATcc                                                              <  2:461041/90‑1 (MQ=255)
                tgCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGGGGGGCGTTATTGCGGATTCg                                                             >  2:995753/1‑90 (MQ=255)
                tgCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGGGGGCCGTTATTGCGCATccc                                                             >  1:411666/1‑89 (MQ=255)
                tgCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGGGGGCCGGTTTTGCGCCTccc                                                             >  1:906408/1‑89 (MQ=255)
                       tGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAgc                                                      <  2:150493/90‑1 (MQ=255)
                         ggCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAgcgc                                                    <  1:1222855/90‑1 (MQ=255)
                               aaTGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCC‑‑CGCAGCATCaa                                            <  2:1502620/90‑1 (MQ=255)
                                   gcgcGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAata                                          <  1:431908/90‑1 (MQ=255)
                                   gcgcGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAata                                          <  2:1311162/90‑1 (MQ=255)
                                     gcgTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGGGGCGGTATCCCCCAACCCCTCCCCCGCCGCATCAatata                                        >  2:409158/1‑89 (MQ=255)
                                                   gCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGGGGGCCGTTTTTGCGCCTCCCCCCCACCGCCGCATCAATAATCGTAAAAacgcgcc                          >  2:579356/1‑89 (MQ=255)
                                                                   tttGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATg          <  2:1428528/90‑1 (MQ=255)
                                                                           aaTCGGCGGGGGGGGCGGTTTTGCGCATCCGCTCCCGCCCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAg  >  2:1193402/1‑90 (MQ=255)
                                                                           aaTCGGCGGGGGGGCCGGTTTCGCGCCTCCGCTCCCACGCCGCAACAAAATGCGTTAAAGGGCGCCCAAAATGCGGATCATGGTGATAAg  >  1:681915/1‑90 (MQ=255)
                                                                                       |                                                                             
GTTCTGATTCAATTGTTGCCAGTTGGGCATCAATGGCGCGTAAATTGTTCAGCAAAAATCCCAGTGTTTTGAGTAAATCGGCGGGTGCGCCGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAG  >  NZ_CP009273/1018286‑1018450

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