Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273982,9500TG62.5% 6.0 / 16.3 17intergenic (‑512/+91)ompF/asnSporin OmpF/asparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/6);  new base G (10/0);  total (11/6)
Fisher's exact test for biased strand distribution p-value = 1.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGG  >  NZ_CP009273/982823‑983071
                                                                                                                               |                                                                                                                         
cAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGaaataaat                                                                                                                >  2:786185/1‑139 (MQ=255)
                    gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTGCTTTACTTaacaaaaaa                                                                                             >  2:111014/1‑135 (MQ=255)
                               tAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTGTTTGAATTTACACAAACATAAAtttata                                                                                  >  2:320241/1‑134 (MQ=255)
                                            aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGCTTCTCATAGATTTCTTTACTTAACCAAAAGAAAAATTTAGAGATCAGATaaaaa                                                                     >  1:5709/1‑138 (MQ=255)
                                                            cGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTTTTTTACTTTAAAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAAACAGCGGGGGCt                                                    >  1:856965/1‑139 (MQ=255)
                                                            cGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGCTAATAGATATATTTGTTTAACAAAAAAAAAACATGAAGGAAAAGATAAAAAAAAAACCCGCGGGGGCt                                                    >  2:372270/1‑139 (MQ=255)
                                                                       ttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGaa                                          <  1:797288/138‑1 (MQ=255)
                                                                          cgcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACg                                      <  1:111014/139‑1 (MQ=255)
                                                                                       cAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGc                         <  1:786185/139‑1 (MQ=255)
                                                                                                 ttCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcg               <  1:795199/139‑1 (MQ=255)
                                                                                                    cGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgcg             <  1:348728/138‑1 (MQ=255)
                                                                                                    cGTAACAGGATGATCGTTATGGCGGGGGTTCAAATAAATTTCTTTACTTTAAAAAAAAATAACTTTAAAGAAAAGAAAAAAAAAAAGCCCGCGGGGGCCGGGGCTTGAAGACGCCAAAATAGAAACCGGCGGTAcgcg             >  2:41370/1‑138 (MQ=255)
                                                                                                      tAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATTCTTGAATTAAAAAAAAAATAAATTTACAGAAAAAATAAAAAAAAAAACAGCGGGGGCTGGCTCTTTAAGACGCCAAATTAAAAAATGGCGGTTCGCgggg          >  2:348728/1‑137 (MQ=255)
                                                                                                      tAACAGGATGATCGTTATGGCGGGGGTGCAAATAGATTGCTTGACTTTACAAAAAGAAAACTTGACAGAAAAGATAAAAAAAAAGCCCGCGGGGGCTGGCCCTTGTAGAAGCCAAAATAGAAACTGGCGGTACGCGGAg          >  1:215860/1‑139 (MQ=255)
                                                                                                      tAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTTCTTGTCTTTAAAAAAAGAAAAAATGACAGAAAAAAAAAAAAAAAAAACAACGGGCGGCGGGTCTTTTAAAAGGTAAATTAGAAACTGGCGGTAAGCgggg          >  2:795199/1‑137 (MQ=255)
                                                                                                            gatgatCGTTATGGCGGGGGTGGAAATAAATTTTTTTACTTTAAAAACAGAAAAATTGAAAGAAAAGATAAAAAAAAAACCCACGGGCGCCGGCTCTTTGAGACGCTAAATTAAAAGCTGGGGTTACGCCGAGTACggg    >  2:486872/1‑137 (MQ=255)
                                                                                                              tgatCGTTATTGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTggg  <  1:372270/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                         
CAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGG  >  NZ_CP009273/982823‑983071

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