Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,043,5780AC54.5% 6.4 / 18.7 23intergenic (‑14/+161)smrA/CJ019_01939putative DNA endonuclease SmrA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (5/5);  new base C (0/12);  total (5/18)
Fisher's exact test for biased strand distribution p-value = 9.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGCAGGCCGGCCAGCTTCTGGCGGTCGGCCTTGGGCTTGCCGACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGGATGTCACAGAGCCCGAAGACGGCGGGGTGCCT  >  KT2440/2043446‑2043688
                                                                                                                                    |                                                                                                              
cgcAGGCCGGCCAGCTTCTGGCGGTCGGCCTTGGGCTTGCCGACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAggcg                                                                                                            <  3:332768/137‑1 (MQ=255)
cgcAGGCCGGCCAGCTTCTGGCGGTCGGCCTTGGGCTTGCCGACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAggcg                                                                                                            <  6:65943/137‑1 (MQ=255)
                     ggggcGGCGTTGGGCTTGCCGACTTCGGCACGGGCGGGGCGTAGCGGCTTTACGCCGCGCACTTCGGCCTTCAACAGGCAAAAATCGTCGTCTTGCATGATGCCCCCACGTCGGCGGCGTAGTTTCCGCGACTGTgg                                                                                       <  2:228824/133‑1 (MQ=255)
                     ggggcGGCCTTGGGCTTGCCGACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGCAAAAATCGTCGTCTTGCATGATGCCTCCACGTCGGCGGCGTAGTTTACGCGACTGTgg                                                                                       <  7:44793/133‑1 (MQ=255)
                     gggTCGGCGTGGGGCTTGCCGACTTCGGCACTGGGGGGGGGGCTCGGTTTAACCCCGCGCAGTTCGGCCTTCACCAGGAAAAAAGCGGCGTCTTGCATGATCCCTCACCGGCGGCGGCGTATTCCCCGCGACTGTgg                                                                                       <  2:325381/136‑1 (MQ=255)
                     gggTCGGCCTTGGGCTTGCCGACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTgg                                                                                       <  6:199273/136‑1 (MQ=255)
                     gggTCGGCCTTGGGCTTGCCGACTGCGCCACGGTCGGGGCGCATCGGCTTTCCGCCGCGCACTTCGGCCTTGAACCGGTAAAAATCGGCGTCTTGCAGGATGCCTCCACGTCGGCGGCGTAGTTTCCGCGACTGTgg                                                                                       <  6:70510/136‑1 (MQ=255)
                     cGGGCGGCGTTGGGCTTGCCGACTTCGGCAGGGTCGGGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGCAAAAATCGTCGTCTTGCATGAGGCCTCCACGTCGGCGGCTTTGTTACCGCGACTGTgg                                                                                       <  8:13758/137‑1 (MQ=255)
                                         gACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGccctgtct                                                                   >  5:220508/1‑131 (MQ=255)
                                         gACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCg                                                                   >  4:53755/1‑137 (MQ=255)
                                         gACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCATGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCg                                                                   >  6:193000/1‑137 (MQ=255)
                                         gACTTCGGCACGGTATTTGAATATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCg                                                                   >  5:161463/1‑137 (MQ=255)
                                                     gTCGGGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTCAACCGGCAAAAAGCGTCGTCTGCCATGATGCCCCCCCGTCGGGGGCGTAGTTTCCGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTcgcg                                                       <  5:370369/137‑1 (MQ=255)
                                                            gggttCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGGCGTCTTGCATGATGCCTCCACGTCGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAAc                                                <  2:256842/133‑1 (MQ=255)
                                                            gggATCGGCTTTACGCCGCGCACTTCGGCCTTCAACAGGCAAAAATCGTCGTCTTGAATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAAc                                                <  6:273496/136‑1 (MQ=255)
                                                            cgtatgggCTTTACGCCGCGCACTTCGGCCTTGAACAGGTAAAAATCGGCGTCTTGCATGATGCCGCCACGTCGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAAc                                                <  8:39857/131‑1 (MQ=255)
                                                                   gCTTTACGCCGCGCACTTGGGCCTTGAACTGGCAAAAATCGTCGTGTTGCAGGATGCAGCCACGTCGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTcc                                         <  6:161463/137‑1 (MQ=38)
                                                                     tttACGCCGCGCAGTTCGGCCTTCACCAGGCAAAATGGGGCGTGTTGCATGATCCCTCCACGTCGGCGGCGTATTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTCCTa                                       <  2:35298/137‑1 (MQ=37)
                                                                     ttaccgccgcGCACTTCGGCCTTCAACTGGTAAAACTCGGCGTCTTGCATCATCCCTCCCCGTCGGCGGCGTAGTTTCCGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTCCTa                                       <  3:53755/133‑1 (MQ=37)
                                                                         cgccgcGCACTTCTGCCTTGAACAGGTAAAAATCGTCGTCTTGCATGTTGCCACCACGTTGGCGGGGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTCCTACggg                                   <  3:342797/137‑1 (MQ=38)
                                                                                 aCTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGgt                                                                                 <  7:211089/83‑1 (MQ=255)
                                                                                 aCTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGgt                                                                                 >  8:211089/1‑83 (MQ=255)
                                                                                                          cggcgacTTGAATCAGGCCTCCAGTTCGGCGGGGTAGTTACCGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTCAAACGGGGATGTCACAGAGCCCGAAGACGGCGGGGTGCCt  <  5:186854/131‑1 (MQ=37)
                                                                                                                                    |                                                                                                              
CGCAGGCCGGCCAGCTTCTGGCGGTCGGCCTTGGGCTTGCCGACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGGATGTCACAGAGCCCGAAGACGGCGGGGTGCCT  >  KT2440/2043446‑2043688

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