Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,331,2020AC52.4% 20.9 / 26.3 21G566G (GGT→GGGgyrADNA topoisomerase (ATP‑hydrolyzing) subunit A
Reads supporting (aligned to +/- strand):  ref base A (8/2);  new base C (0/11);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 2.21e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACA  >  NZ_CP009273/2331124‑2331286
                                                                              |                                                                                    
gCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGcc                                                                           >  2:413434/1‑90 (MQ=255)
     agatatGGTCGTGAGTGTTCGCCCCCAGCAGTCGGTCGAAAAATTCTTTTTTTTTAATGCGGGCGGCATATTTACCTTTCCCGCCacgac                                                                      <  1:23868/88‑1 (MQ=255)
       atatGGTCGTGGGTGTTCGCCCCCAGCAGTGGGTCAAAAAAGTCTTTTTCTTTAATCCGTGCGGCATATTTCCCTTCCCCGCCACGACGc                                                                    <  1:597115/90‑1 (MQ=255)
       atatGGTCGTGAGTTTTCCCCACCAGCAGTCGGTCGAAAAAGTCTTTTTTTTTAATACGGGCGGCAGTTTTCCCTTTCCCGCCACGACGc                                                                    <  1:625245/90‑1 (MQ=255)
         atGGTCGTGATTGTCCCCCACCAGCGGCGGGTCAAAAAATTTTTTTTTTTAAATCCGTGCGGCATTTTTCCCTTCCCCGCCACGACGCTg                                                                  <  2:460316/90‑1 (MQ=255)
               gTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTc                                                            >  1:75449/1‑90 (MQ=255)
               gTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTc                                                            >  2:367510/1‑90 (MQ=255)
               gTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTc                                                            >  2:271428/1‑90 (MQ=255)
                 gAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTcct                                                          >  2:23868/1‑88 (MQ=255)
                         ccccccaGCAGTGGGTCAAAAAAGTCTTTTTTTTTAATCCGTGCGGCAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGaa                                                  <  1:367510/86‑1 (MQ=255)
                          ccaccaGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGaaa                                                 >  1:618779/1‑90 (MQ=255)
                           ccccaGCAGTGGGTCAAAAAAGTTTTTTTTTTTAATACGGGGGGAAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAg                                                <  1:679804/88‑1 (MQ=255)
                                    tcggtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAc                                                                                     >  2:587860/1‑44 (MQ=255)
                                    tcggtcgATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTAc                                                                                     <  1:587860/44‑1 (MQ=255)
                                          gaaaaaTTCTTTTTTTTTAATGCGTGGGGCATTTTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAAc                                 <  2:75449/87‑1 (MQ=255)
                                                  cttcttcttTAATCCGTGCGGCAGATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCt                         <  1:181021/90‑1 (MQ=255)
                                                          ttAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgagag                 >  2:122314/1‑90 (MQ=255)
                                                            aaTCCGGGCGGCAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGagagag               <  1:611491/90‑1 (MQ=255)
                                                                        agttttCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACcac   <  1:653262/87‑1 (MQ=255)
                                                                         gttttCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACcaca  <  1:413434/88‑1 (MQ=255)
                                                                         gATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACcaca  <  1:122314/90‑1 (MQ=255)
                                                                              |                                                                                    
GCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACA  >  NZ_CP009273/2331124‑2331286

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