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

CGCAGGCCGGCCAGCTTCTGGCGGTCGGCCTTGGGCTTGCCGACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGGATGTCACAGAGCCCGAAGACGGCGG  >  ACB143/2043446‑2043681
                                                                                                                                    |                                                                                                        
cgcAGGCCGGCCAGCTTCTGGCGGTCGGCCTTGGGCTTGCCGACTTCGGCACGGGCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAggcg                                                                                                      <  3:143689/137‑1 (MQ=255)
        ggccAGCTTCTGGCGGTCGGCCTTGGGCTTGCCGACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAggcgg                                                                                                     >  7:253589/1‑130 (MQ=255)
        ggccAGCTTCTGGCGGTCGGCCTTGGGCTTGCCGACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAggcgg                                                                                                     <  8:253589/130‑1 (MQ=255)
            tgCTTCGGGGGGGCGGCGTTGGGCTGGCCGACTTCGGCAGGGGGGTGGGGGATCGGTTTTACGCCGCGAACTTCGGCCTTCAACAGGTAAAAATCGTCGTCTTGCATGATGCCTCCACGTCGGCGGCGTAGTTAAcg                                                                                         <  4:315822/136‑1 (MQ=255)
            aggttCTGGCGGTCGGCCTTGGGCTTGCCGACTTCGCCACGGTGGGGGCGTTTCGGATTAACCCCGCGCACTTCGGCCTTGACCTGGTAAACAGCGCGTTCTTGCATGATGCCTCCACTTCGGCGGCGTAGTTAAcg                                                                                         <  6:160076/134‑1 (MQ=255)
                   ggCGGTCGGCCTTGGGCTTGCCGACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTTACCAGGTAAAAATCGTCGTCTTCCATGATGCCTCCACGTAGGCGGCGTAGTTAACGCGACtgt                                                                                  <  8:180811/137‑1 (MQ=255)
                     ggggcGGCCTTGGGCTTGCCGACTTCGGCACGGGCGGGGCGGATCGGCTTTACGCCGCGCACTGCGGCCTTGAACCGGGAAAAATCGTCGTCTTGCATGAGGCCCCCACGTCGGCGGCGTAGTTTACGCGACTGTgg                                                                                <  8:125557/133‑1 (MQ=255)
                           gCCTTGGGCTTGCCGACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGgt                                                                          <  7:308261/137‑1 (MQ=255)
                           gCCTTGGGCTTGCCGACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGgt                                                                          >  2:62425/1‑137 (MQ=255)
                               ggggCTTGCGGACTTCGGCTGGGTCGTGGGGTTTCTGTTTTACGCCGCGATCTTGGGCTTTTAACTGGAAAATTTGGTGGTCTTGATTGAGGCCTCCAGGTTGGGGGGTTAGTAACCGCGACTGTGGGGCGGTGTGc                                                                      <  4:4798/136‑1 (MQ=255)
                                     tGCCGACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTCAACCGGCAAAAAGCGTCGTCTTGCATGATGCCGCCACGTCGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCtg                                                                <  6:80468/137‑1 (MQ=255)
                                         gACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCg                                                            >  3:145244/1‑137 (MQ=255)
                                         gACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCg                                                            >  1:274351/1‑137 (MQ=255)
                                               ggCACGGGCGTGGCGGAGCGGCTTAACGCCGCGCACTTCGGCCTTCAACAGGGAAAAATCGTCGTCTTCCATGAGGCCTCCACGTCGGCGGCGTAGTTTCCGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCctc                                                      <  7:1052/137‑1 (MQ=255)
                                                gCAGGGTGGGGGCGGATCGGCTTAACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTCCATGATGCCTCCAGGTCGGCGGCGTAGTTTCCGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCctct                                                     <  2:72648/137‑1 (MQ=255)
                                                            cGGATCGGCTTAACGCCGCGCAGTTCGGCCTTGAACAGGCAAAAATCTTGGTCTTGCCTGAGGCCTCCACGTAGGGGGCGTATTTACCGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAAc                                         <  7:65528/137‑1 (MQ=255)
                                                             gtatgggCTTTACGCCGCGCACTGCGGCCTTTAACAGGCAAAAATCGGCGTCTTGCATGATCCCTCCCCGTCGGCGGCGTAGTTAACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAAcc                                        <  4:255093/132‑1 (MQ=255)
                                                                ggggCTATACGCGGCGAAGTTGGTCGTTGAACTGGTAAAAATCGTCGTATTTCATGATGCCACAACGTCGGCGGGTTATTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGc                                     <  4:190020/135‑1 (MQ=255)
                                                                   gCTTTCCGCCGCGCACTTGGGCGTTGAACAGGTAAATCTCGTCGTCTTGCATGAGGCCTCCCCGTCGGGGGCGTTGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTcc                                  <  1:289924/137‑1 (MQ=37)
                                                                     ttcaCGCCGCGCACTTCGGCCTGCAACCGGGAAAACTCGTCGTCTTCAATGATGCCTCCACGGCGGGGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTCCTa                                <  6:72840/134‑1 (MQ=37)
                                                                        ccgccgcGCACTTCGGCCTTGAACCGGGAAAAAGCGTCGGCTTGCAGGAGGCCTCCACGTCGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTCCTACgg                             <  6:332458/136‑1 (MQ=38)
                                                                              cgcACTGCGGCATTAAACCTTCAAAAATCGGCGTCTTCCTTGAGGCCGCCCCTTGGGCGGCTTTTTTCACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGGATGt                       <  5:243655/137‑1 (MQ=16)
                                                                                                   aaaaaaGCGGGGTCTTCGATGATGCCCCCTCTTTGGGGGGTTTTTTACGGCTACTGTGGGGCGGTGTGCAGCCGGTGCGGGCCTCTTCGCGGTTGAACCGGCTCCTACGGGGATGTCACAGACCCCGATGAcggcgg  <  4:124901/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                        
CGCAGGCCGGCCAGCTTCTGGCGGTCGGCCTTGGGCTTGCCGACTTCGGCACGGTCGTGGCGGATCGGCTTTACGCCGCGCACTTCGGCCTTGAACAGGGAAAAATCGTCGTCTTGCATGATGCCTCCACGTAGGCGGCGTAGTTTACGCGACTGTGGGGCGGTGTGCAGCCTGTGCGGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGGATGTCACAGAGCCCGAAGACGGCGG  >  ACB143/2043446‑2043681

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: