Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,128,8120TG61.3% ‑1.2 / 11.3 32intergenic (‑47/+18)CJ019_04823/aer_3hypothetical protein/Aerotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (9/3);  new base G (19/0);  total (29/3)
Fisher's exact test for biased strand distribution p-value = 4.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CTTTGCGGG‑CCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAAT  >  KT2440/5128681‑5128946
                                                                                                                                    |                                                                                                                                       
gggTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CGTTggg                                                                                                                                    >  5:63470/1‑135 (MQ=255)
     cGGCTGTCGAGCGAGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGACCG‑CGTTTCGGG‑gcc                                                                                                                              >  2:377044/1‑137 (MQ=255)
      ggCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CGTTTGGGG‑cccc                                                                                                                             >  7:136055/1‑137 (MQ=255)
      ggCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CGTTGCGGG‑cccc                                                                                                                             >  7:86553/1‑137 (MQ=255)
              aGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CGTTGGGGG‑GCCCGCGAgttt                                                                                                                     >  4:17524/1‑136 (MQ=255)
                     ggATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CTTTGCGGG‑CCCCGCGATTTCAGCGGtt                                                                                                              >  8:10105/1‑137 (MQ=255)
                     ggATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CGTTGCGGG‑CCCCGCGATTTCAGCgggg                                                                                                              >  2:168926/1‑135 (MQ=255)
                        tCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CTTTGCGGG‑CCCCGCGCTTTCAGCGGTTGaa                                                                                                           >  5:291221/1‑137 (MQ=255)
                                  tAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CGTTGCGGG‑CCCCGCGCGTTCAGCGGGTGAAACGTTCCGcc                                                                                                 >  3:238798/1‑137 (MQ=255)
                                    ggCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CGTTGCGGG‑CCCCGCGATTTAAGCGGTTGAACCTTTCCGCCac                                                                                               >  4:92176/1‑136 (MQ=255)
                                      cAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CGTTGCGGG‑CCCCGCGCTTTGAGCGGTTGTACCTTTCCGCCCGGc                                                                                             >  3:106349/1‑137 (MQ=255)
                                           gTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑GAGTGCGGG‑CCCCGCGCTTTCAGCGGTTGTAGCTTTACGCCCggggttgc                                                                                        >  7:234023/1‑131 (MQ=255)
                                           gTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CTTTGCGGG‑CCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATgg                                                                                        >  7:41384/1‑136 (MQ=255)
                                            tCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑GGTTGGGGG‑CCCCGCGGTTTCAGCGGGTGGAACGTTACGGCAGGCGATggc                                                                                       >  7:47308/1‑135 (MQ=255)
                                            tCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CTTTGCGGG‑CCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGGTAttcc                                                                                       >  7:265141/1‑134 (MQ=255)
                                            tCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CGTTGCGGG‑CCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCa                                                                                       >  7:185402/1‑137 (MQ=255)
                                             cGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CGTTGCGGG‑CCCCGCGATTTCACCGGGTGAACCGTTACGCCCGGCTATggca                                                                                      >  5:46614/1‑134 (MQ=255)
                                                        tGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CGTTGCGGG‑CCCCGCGCTTTCCGCGGGTGACCCTTTCCGCCAGGCTATGGAGGTACTACGccc                                                                           >  7:79390/1‑136 (MQ=255)
                                                         ggCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCC‑CGTTGGGGG‑CCCCGGGGTTTGAGCGGGTGGACCGTTACGCCCGGCTATGGAGGTTCAGCGCCAt                                                                          >  7:73120/1‑137 (MQ=255)
                                                                                   aTGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CGTTGCGGG‑CCCCCCGCTTTCAGCGGGTGCACCGTTACGCCCGGGTATGCCGGTACTCCGCCCTGGGCGCGAGGTCCTGGCTGCTTTCcg                                                >  7:211072/1‑137 (MQ=255)
                                                                                        ggTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CTTTGCGGG‑CCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGGTATGCAGGTACTCCGCCATGCGCTCGAGGTACTGGCTGATTTCCGCGCCt                                           >  8:352455/1‑137 (MQ=255)
                                                                                        ggTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CTTTGCGGG‑CCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTACGCGCCt                                           >  7:332409/1‑137 (MQ=255)
                                                                                        ggTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CGTTTCGGG‑CCCCGCGATTTGAGCGGTTGAACCATTACGCCCCGCGATGCACGGAATCCGCCAATGGCTCGACGTGCAGGGGGAATTCCGCGcca                                           >  4:343999/1‑136 (MQ=255)
                                                                                         gTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CTTTTCGGG‑GCACGCGATTTGATCGGTTTAACCTTTACGCCAGGCTATTCAGGTTCTCCCCCATGCGCGCGAGGTCCTGGGTGATTTACGCGCCtt                                          >  8:158238/1‑137 (MQ=255)
                                                                                         gTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CGTTGCGGG‑CCCCGCGATTTCAGCGGTTGAACCTTTACGCCAGGCTATGCAGGTACTCCGCCATGCGCGCGAGGTCCTGGCTGATTTCCGCGCCtt                                          >  6:158214/1‑137 (MQ=255)
                                                                                                         cAGTTTGCCACATAACCGGGGCCCG‑CGTTGCGGG‑CCCCGCGCTTTGAGCGGTTGAACCTTTACGCCAGGCTATGCAGGGACTCCGCCAAGCGCGCGAGGTCCTGGCGGATTTCCGCGCCATGCCTGGCATGCCCGGc                          >  4:71467/1‑137 (MQ=255)
                                                                                                           gTTTGCCACATAACCGGGGCCCG‑CGTTTCGGG‑CCCCGCGATTTCAGCGGTTTAACCTTTACGCCCGGCTATGCAGGTACTCCGCCATTCGCTCGAGGTCCTGGCTGGTTTGCGCGCCTTGCCTGGCTTGCCCGGccc                        >  7:116527/1‑135 (MQ=255)
                                                                                                                   aaaacccGGGGCCCGTTTTTGGGGG‑CCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTc                 <  6:291221/132‑1 (MQ=255)
                                                                                                                      acccGGGGCCCG‑GTTGGGGGGCCCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGca              <  4:106349/135‑1 (MQ=255)
                                                                                                                         cGGGGCCCG‑CTTTGCGGG‑CCCCGCGATTTCAGCGGTTGAACCTTTACGCCACGCTTTGCAGGTACTCCGCCATGCGCTCGGGGTGCTGGGGGATTTGCGCGCCGTGCCTCGCTTGCCCGGCACTCGGGTCGGACAGc          >  6:278796/1‑137 (MQ=255)
                                                                                                                          gggcCCCG‑TTTGGCGGG‑CCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCt         <  6:316413/137‑1 (MQ=255)
                                                                                                                                 g‑CTTTGCGGG‑CCCCGCGATTTAAGCGGTTGAACCTTTACGCCCGGCGATTGAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAAt  >  2:263571/1‑137 (MQ=255)
                                                                                                                                    |                                                                                                                                       
GGGTGCGGCTGTCGAGCGCGGGGATCAGTGGGCATAGGCAATGGTCGAGCGGGCGCTGGCAGCGCTCGCAGCGGGGTCTGGGCATGGCGGTCTCCTTGGCGGGGACAGTTTGCCACATAACCGGGGCCCG‑CTTTGCGGG‑CCCCGCGATTTCAGCGGTTGAACCTTTCCGCCAGGCTATGCAGGTACTCCGCCATGCGCTCGAGGTCCTGGCTGATTTCCGCGCCTTGCCTGGCTTGCCCGGCACTCTGGTCGCACAGCTGGGCAAT  >  KT2440/5128681‑5128946

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