Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,124,7990TC53.8% 9.7 / 17.7 13pseudogene (969/1065 nt)yghQlipopolysaccharide biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base T (6/0);  new base C (0/7);  total (6/7)
Fisher's exact test for biased strand distribution p-value = 5.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCA  >  NZ_CP009273/3124666‑3124896
                                                                                                                                      |                                                                                                 
ccGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTccgcc                                                                                              >  2:90814/1‑139 (MQ=255)
                     ttccccACCTCCCAAAAACCAGGCCATGCGCTGCGAAATTCTAACCGCAAAACACCAGAAAATGGGGGGCCTCGCCGCCTACCACCCCCGCGCGCCCCACCAAAACACCGTTCCCCCCCGCCAGTAAACCGGACTTca                                                                          <  2:561978/133‑1 (MQ=255)
                                 ccaaaaaCAGGTCATCCGCTTCGAGATACTTAACGCCAAACCCCAGCGAAAGGGGCGTCTTGCCGACGATCAACCCCGCGAGCCCCACCAGAATCCCGATCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCa                                                             <  1:41752/136‑1 (MQ=255)
                                 ccaaaaCCAGGCCATACGCTCCGAGATACTAAACGCCAACCCCAAGCGAAATGGGCGGCTTGCCGACAACCAACACCGCGGGCGCCACCAGAACACCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCa                                                             <  1:588534/136‑1 (MQ=255)
                                      ataaGGTCACCCCCTTGGAGATACTAACCGCCAAACCCCAGCAAAAGGCGGGGCCCGCCGACGACCCCCCCCGCGAGCCCCACCAAAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTc                                                        <  2:357548/136‑1 (MQ=255)
                                           gTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCtggtg                                                   >  1:561893/1‑139 (MQ=255)
                                           gTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCtggtg                                                   >  1:576211/1‑139 (MQ=255)
                                                     tCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGAt                                         >  1:245631/1‑139 (MQ=255)
                                                          atacaaaaccccaaccccCAGCGAAGTGCGGGGCTGCCCGCCGACCAACCCGGCGGCCGCCCCCAGAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTaaa                                    <  2:561893/133‑1 (MQ=255)
                                                          aTACTAACCCCCAAACACAAGAAAAATGGGCGGCTGGCCGACGACCAACCCCGCGACCCCCACAAAAATCCCGTTCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTaaa                                    <  2:245631/139‑1 (MQ=255)
                                                                       aaCACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTg                       >  1:69661/1‑139 (MQ=255)
                                                                             aGCAAAATGGGCGGCTGCCCGCCTACAACCCCCGCGAGCCCCCCCAGAATCCCGATCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTaaa                 <  1:458147/139‑1 (MQ=255)
                                                                                              tGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCa  >  1:467162/1‑138 (MQ=255)
                                                                                                                                      |                                                                                                 
CCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCA  >  NZ_CP009273/3124666‑3124896

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: