Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,053,9480AC54.5% 3.8 / 10.7 11T137P (ACC→CCC) torCpentaheme c‑type cytochrome TorC
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base C (0/6);  total (4/7)
Fisher's exact test for biased strand distribution p-value = 1.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA  >  NZ_CP009273/1053872‑1054010
                                                                            |                                                              
acCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTg                                                   >  1:173752/1‑90 (MQ=255)
aaccgAAAAATTCGAAGCCAAACTCGCGAACCTTGCCGCGGATAAATGGGCGCGCATGGAAGAAAACAACTCGACAACCTGCCGCTCCTg                                                   >  2:688431/5‑90 (MQ=255)
  cggAAAAATTCAAACCCAACCGCGCGGAACTTGCCGGGCGTGAATGGGCGCAAATGAAAAAAAACAACTCGGCACCCTGCCGCTCCTGcc                                                 <  2:189373/88‑1 (MQ=255)
              aaGCCAACCGCGCGGACCTTGCCGAGCGTGAAGGGGCGCGAATAAAAAAAAACAACTCGGCACCCTGCCGCTCCTGCCATAACTAcgatg                                     <  1:461010/90‑1 (MQ=255)
                 ccAACCGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACACCTCGGCACCCTGCCGCTCCTGCCATAACTAcgatgcga                                  <  1:246330/90‑1 (MQ=255)
                 ccAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTAcgatgcga                                  <  1:38221/90‑1 (MQ=255)
                 caaaaCGCGCGAAACTTGCGGGGCGTAAAGGGGCGCAAATAAAAAAAAACAACCCGGCACCCTGCCGCCCCTGCCATAACTAcgatgcga                                  <  2:84786/88‑1 (MQ=255)
                                   ccGGGCGTAAAGGGGCGCAAAGAAAAAAAAACACCCCGCCACCCTGCCGCCCCTGCCATAACTACGATGCGATGGATCATGCGAagcagc                <  2:195769/90‑1 (MQ=255)
                                   ccGGGCGTAAAGGGGCGCAAAAAAAAAAAAACAACCCGGCACCCTGCCGCCCCTGCCATAACTACGATGCGATGGATCATGCGAagcagc                <  2:195194/90‑1 (MQ=255)
                                      aGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATc             >  1:708251/1‑90 (MQ=255)
                                                 gcgcGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAgcagca  >  2:410765/1‑90 (MQ=255)
                                                                            |                                                              
ACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA  >  NZ_CP009273/1053872‑1054010

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: ‑