Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273653,3560AC57.9% 17.2 / 21.4 19G13G (GGT→GGGcrcBfluoride efflux transporter CrcB
Reads supporting (aligned to +/- strand):  ref base A (8/0);  new base C (0/11);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 1.32e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-05
Rejected as consensus: Frequency below/above cutoff threshold.

ATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAAC  >  NZ_CP009273/653267‑653437
                                                                                         |                                                                                 
aTCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTa                                                                                   >  1:449268/1‑90 (MQ=255)
                 cgcccccAACGGAATCGCCTGATGCAGGGGTTTAAATCGCATACTTAACAGCCATCTCGCCCCCCTCCCCGTCCCACCGCCAATAAAAAc                                                                  <  2:310463/87‑1 (MQ=255)
                    ccccAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGc                                                               >  1:615255/1‑90 (MQ=255)
                    ccccAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGGACCACCGCCAATAAAAACTGc                                                               >  1:580253/1‑90 (MQ=255)
                            gAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAaagaa                                                       >  2:68571/1‑90 (MQ=255)
                            gAATCGCCTGATGCAGTGGGTAAAACCGCATATTTAACAGCCATCTCGCCACGCTCCCCGTCCCACCGCCAATAAAAACTGCTAAaagaa                                                       <  2:583152/90‑1 (MQ=255)
                            aaaTCCCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCCCCCCGCTTCCCGTCCCACCGCCAATAAAAACTGCTAAaagaa                                                       <  2:214250/89‑1 (MQ=255)
                                   cTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTaa                                                >  2:561007/1‑90 (MQ=255)
                                           gTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCaa                                        >  2:335115/1‑90 (MQ=255)
                                                       ccATACTAAACAGCCATCTCCCCACGCTTCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTatt                            <  2:449268/89‑1 (MQ=255)
                                                                 cacccATCTCGCCCCGCTTCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTAt                  <  2:573549/87‑1 (MQ=255)
                                                                 cacccACCTCCCCCCGCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTAt                  <  1:68571/87‑1 (MQ=255)
                                                                   gCCATCCCGCCCCTCTTCCCGTCCCCCCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGa                <  1:284921/90‑1 (MQ=255)
                                                                      aTCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTg             >  1:427218/1‑90 (MQ=255)
                                                                      aTCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTg             >  2:598273/1‑90 (MQ=255)
                                                                      aTCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTg             >  2:760661/1‑90 (MQ=255)
                                                                         cccccACCCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATa          <  1:129362/87‑1 (MQ=255)
                                                                         ccGCCACGCTTCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATa          <  2:745585/89‑1 (MQ=255)
                                                                         ccGCCACGCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATa          <  2:21646/89‑1 (MQ=255)
                                                                                 ctccccGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAAc  <  2:13975/87‑1 (MQ=255)
                                                                                         |                                                                                 
ATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAAC  >  NZ_CP009273/653267‑653437

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