Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,288,0410TG52.4% 8.2 / 17.6 21intergenic (+4/+66)slyB/nhaA_1Outer membrane lipoprotein SlyB/Na( )/H( ) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base T (5/5);  new base G (11/0);  total (16/5)
Fisher's exact test for biased strand distribution p-value = 1.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACG  >  ACB122/1287936‑1288174
                                                                                                         |                                                                                                                                     
aTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGGTTCGCGCGGGCCGGTGTGGAGTAATAAAAAAACCCGCCAGGgc                                                                                                          >  6:206097/1‑135 (MQ=255)
  gcgcgcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGGACGGGCCGGAGTGGGGTAATTAAAAAAACCGACAAGgctg                                                                                                        >  1:382232/1‑135 (MQ=255)
    gcgcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGTCCGGACTGGAGTAAAAAAAAAAACCCACAAGgctgct                                                                                                      >  7:90052/1‑135 (MQ=255)
              gCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGGTTAGCACTCCCCGGACTGGAGTAAAAAAAAAACCCCACAAGGCGGGGGCCACTCggg                                                                                            >  3:381003/1‑135 (MQ=255)
                        cAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTGGCCCTGACCGGAGTGCAGTAATAAAAAACCCCCCCAAGGCTGCTGCCAGTCGGGGGGTTTTGAt                                                                                  >  3:92254/1‑135 (MQ=255)
                         aGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGGTTCGCACTGACCGGAGTGCAGTAATAAAAAAACCCGCCAAGGCTGCGGCCAGTCGGGGGTTTTTGAtt                                                                                 >  5:158867/1‑134 (MQ=255)
                         aGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGGTTCGCACTGACCGGAGTGCAGTAATAAAAAAACCCCCCAAGGCGGCTGGCCCTCGGGGGTTTTTTAtt                                                                                 >  6:291221/1‑134 (MQ=255)
                          ggAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCCCTGCCCGGAGTGCAGGAATAAAAAACCCCCACATGGCTGCTGCCAGGCGGGGGTTTTTAATgg                                                                                >  6:83622/1‑134 (MQ=255)
                                 gAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTTTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGcc                                                                         <  2:115159/135‑1 (MQ=255)
                                 gAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCAATGAACGTAGTGGAGTAAAAAAAAAACCCGCCATTGCCGCTGCCACTAGGGGTTTTTTAAATCCTCAga                                                                          >  7:72301/1‑133 (MQ=255)
                                        tCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGGAGTGGAGTAATAAAAAAACCCGACAATGCTGCTGCCAGTCGGGGTTTTTTCATTCCGCAGCCAAGTCAg                                                                  >  4:252746/1‑135 (MQ=255)
                                        tCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGGTTGGCACTGGCCGGAGTGGAGTAATAAAAAAACCCCCCCAGGCTGCGGCCAGGGGGGGGGTTTTTATTGCTCAGCCaacatag                                                                  >  4:203343/1‑130 (MQ=255)
                                                                      cTGTCGACGGTACCCGGCGCGTGTCGCACTGACCGTAGTGCATTAAAAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCg                                    <  1:349076/135‑1 (MQ=255)
                                                                            aCGGTACCAGTCGCGGTTAGCACTGCCCGTAGTGCAGTAATAAAAAAACCCGACAAGGCTGCTGCCAGTCGGGGTTTTTTGATTCCTCACCCAAGTCAGGCAATGCTGGGTTTGCGGGTGGCCATTGCCGTTAcc                              >  8:209221/1‑134 (MQ=255)
                                                                              ggTACCAGTCGCGTTTCGCACTGACCGGAGTGGAGTAATAAAAAACCCCGCCATGGCTGCTGCCAGTCGGGGGTTTTTCATGCCTCACCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCCTTGCGGTggcggc                            >  7:319245/1‑135 (MQ=255)
                                                                                        gcgTTTAGCACGGACCGTAGTGCAGTAATAAAAAACCCCGACAAGGCTGCTGGCAGTCGGGGGTTTTTAATGCCTCAGCCAAGTCAGGCAATTCTGGTTTTTCGGCTGGCCCATTCCGTGACCGCGGTGCCGCTc                  >  2:10577/1‑135 (MQ=255)
                                                                                         cgTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCGGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCa                 <  8:65020/135‑1 (MQ=255)
                                                                                                      cggTAGGGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAga     <  8:72301/132‑1 (MQ=255)
                                                                                                      ccGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAga     <  6:307158/134‑1 (MQ=255)
                                                                                                        gTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaac   >  8:38397/1‑134 (MQ=255)
                                                                                                        gTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTACTGGCTCAGCCAAGTCAGGTAATGCTGGTTTTGCGGCTGGCCATTGCGGTGCCGGGGGTGCCGGTCGCCGCAGGCAAgcacg  >  3:64799/1‑135 (MQ=255)
                                                                                                         |                                                                                                                                     
ATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCAGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACG  >  ACB122/1287936‑1288174

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