Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273766,3680AT53.3% 3.7 / 15.0 15intergenic (+301/‑546)mngB/cydAmannosylglycerate hydrolase/cytochrome ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base T (8/0);  total (11/4)
Fisher's exact test for biased strand distribution p-value = 2.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCAAAAAAAGAATTAAGGTCAACCGTGCTGTTTTTGCTTCGTCTCTTTT  >  NZ_CP009273/766285‑766455
                                                                                   |                                                                                       
aaaTAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTa                                                                                   <  2:753310/90‑1 (MQ=255)
             gggTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGGGGATTTTTACCTTTTTTATTAATACTTAAActct                                                                      >  1:156922/1‑90 (MQ=255)
                 gAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGcgt                                                                  <  1:139247/90‑1 (MQ=255)
                   gTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGGGGGTTTTTTTCTTTTTTATTAATACTGAAAATCTGcgtcg                                                                >  2:47555/1‑90 (MQ=255)
                                  atTATAAGTTAACTAAATGTTAATATTGGCGGGGGGGATTTTTTCCTTTTTTTGTAATCCTGAAAATCTGCGTCGTATTAGCCAGTGAcc                                                 >  2:24121/1‑90 (MQ=255)
                                        aGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCaaaaaa                                           <  2:846422/90‑1 (MQ=255)
                                           tAACTAAATGTTAATATTGGCGGGGGGGGTTTATTCCTTTTTTTTTAATACTGAAAATTTTCGGCGGATTAGGCAGTGGCCaaaaaaaaa                                        >  2:557785/1‑88 (MQ=255)
                                            aaCTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAAACCTGAAAATCTTCGGCGTATTAGCCAGGGACCAAAAAAAGaa                                       >  1:615056/1‑90 (MQ=255)
                                            aaCTAAATGTTAATATTGGCGGGGGGGGTTTATTCCTTTTTTTTTAAAACTTAAAATTTGCGGCGTATTAGTCAGTGACCaaaaaaaaaa                                       >  2:723680/1‑90 (MQ=255)
                                            aaCTAAATGTTAATATTGGCGGGGGGGATTTTTGCCTTTTTTAGTAATACTGAAAAACTTCGTCGTATTAGCCAGTGACCAAAAAAAGaa                                       >  2:283/1‑90 (MQ=255)
                                                       aaTATTGGCGGGGGGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCAAAAAAAAAATTAAGGTCAAc                            >  1:690691/1‑90 (MQ=255)
                                                            tGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCAAAAAAAGAATTAAGGTCAACCGTGc                       <  2:845204/90‑1 (MQ=255)
                                                            tGGCGGGGGGGGTTTTTTCCTTTTTTAGTAATCCTGAAACTCTGCGTCGGATTAGCCCGTGACCAAAAAAAAAATTAAAGTCAACCGTGc                       >  2:703039/1‑90 (MQ=255)
                                                             ggCGGGGGGGGTTTATGCCTTTTTTAGTAATCCTGAAAATCTTCGTCGGATTAGCCAGTGACCAAAAAAAAAATTAAAGTCAACCGTGCt                      >  2:881569/1‑90 (MQ=255)
                                                                               cTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCAAAAAAAGAATTAAGGTCAACCGTGCTGTTTTTGCTTCGTCTCtt    >  1:60031/1‑90 (MQ=255)
                                                                                 ttattaGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCAAAAAAAGAATTa                                    >  2:361071/1‑56 (MQ=255)
                                                                                 ttattaGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCAAAAAAAGAATTa                                    <  1:361071/56‑1 (MQ=255)
                                                                                 ttattaGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCAAAAAAAGAATTAAGGTCAACCGTGCTGTTTTTGCTTCGTCTCtttt  >  1:754140/1‑90 (MQ=255)
                                                                                   |                                                                                       
AAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGACCAAAAAAAGAATTAAGGTCAACCGTGCTGTTTTTGCTTCGTCTCTTTT  >  NZ_CP009273/766285‑766455

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