Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,388,9230AC54.2% 24.1 / 28.5 24G151G (GGT→GGGnuoJNADH‑quinone oxidoreductase subunit J
Reads supporting (aligned to +/- strand):  ref base A (7/4);  new base C (0/13);  total (7/17)
Fisher's exact test for biased strand distribution p-value = 9.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-04
Rejected as consensus: Frequency below/above cutoff threshold.

TTCGCGCTGTCGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGC  >  NZ_CP009273/2388846‑2389010
                                                                              |                                                                                       
ttcGCGCTGTCGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGagcagc                                                                             <  2:540595/90‑1 (MQ=255)
                  ccaTTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGGGGAAGGCCACCACCACACCTGCGAGCAGCAGCATAGAAGCCAGGTcc                                                           >  1:377143/3‑90 (MQ=255)
                  cGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTcc                                                           >  2:173775/1‑90 (MQ=255)
                  cGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTcc                                                           >  2:339473/1‑90 (MQ=255)
                          cAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGagcagc                                                                             <  2:26129/64‑1 (MQ=255)
                          cAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGagcagc                                                                             >  1:26129/1‑64 (MQ=255)
                                tcccccccCACGCTCTTCACGCCCGACGTGGAAGGCCACAACCAGCCCTGCGAGCAGCAGCATAGAAGCCAGTCCCACCGCCAGTACGTa                                             <  2:522647/86‑1 (MQ=255)
                                     ccgcccgcTCTCCACGCCCGACGGGGAAGCCCACAACCAGCCCGGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGcc                                        <  1:173775/90‑1 (MQ=255)
                                          ccctctCCCCGCCCGACGGGAAAGCCCAAAACCAGCCCGGCGAGCAGCAGAAAAAAAGCCATTTCCACCGCCAGTACGTAAGGCCCGAAc                                   <  1:262995/88‑1 (MQ=255)
                                          cGCTTTTCACGCCCGCCGGGGAAGCCCACACCCAGCCCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAAc                                   <  2:454685/90‑1 (MQ=255)
                                          cGCTCTTCACGCCCGACGTGGAAGGCCACAACCAGCCCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAAc                                   <  2:282072/90‑1 (MQ=255)
                                          cGCTCTTCACGCCCGAGGGGAAGGCCCACAACCAGCCCGGGGAGAAGAAGAATAAAAGCCATTTCCACCGCCAGTACGTAAGGCCCGAAc                                   <  1:393117/90‑1 (MQ=255)
                                          cGCTCTTCACGACCGACGGGGAAGGCCACACCCAGCCCGGCGAGCAGCAGCATAAAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAAc                                   <  1:38859/90‑1 (MQ=255)
                                                 cACGCCCGACGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTaa                            <  2:138470/90‑1 (MQ=255)
                                                 cACGCCCGACGGGGAAGGCCACAACCAGCCCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTaa                            <  2:523032/90‑1 (MQ=255)
                                                 aaCGCCCGCCGGGGAAGCCCACACCCAGCCCGGCGAGCAGCAGAAAAAAACCCATTCCCACCGCCAGTACGTAAGGCCCGAACAGCGTaa                            <  1:272205/89‑1 (MQ=255)
                                                      ccgacGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATAccc                       >  2:32168/1‑90 (MQ=255)
                                                        gacGTGGAAGGCCACACCCAGCCCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCAc                     <  2:247457/90‑1 (MQ=255)
                                                        gacGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGcc                                                     >  1:443653/1‑58 (MQ=255)
                                                        gacGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGcc                                                     <  2:443653/58‑1 (MQ=255)
                                                         ccGTGGAAGCCCACAACCAGCCCGGCGAGCAGCAACATAAAAGCCATTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACt                    <  1:2075/89‑1 (MQ=255)
                                                         acGTGGAGGCCCACAACCAGCCCGGCGAGCAGCAGCAAAAAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACt                    <  1:141278/90‑1 (MQ=255)
                                                         acGGGGAAGGCCAAAACCAGCCCTGCGAGCAGCAGCATAAAAGCCAGTCCCCCCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACt                    <  2:470623/90‑1 (MQ=255)
                                                                            aGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGc  >  1:280055/1‑90 (MQ=255)
                                                                              |                                                                                       
TTCGCGCTGTCGTCTTTACGATTGCTCAGCACTTCACCCGCACGCTCTTCACGACCGACGTGGAAGGCCACAACCAGACCTGCGAGCAGCAGCATAGAAGCCAGTTCCACCGCCAGTACGTAAGGCCCGAACAGCGTAATACCCACTGCTTTAGCACTGATTGGC  >  NZ_CP009273/2388846‑2389010

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