Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,028,4280CA25.6% 76.1 / 11.9 41intergenic (‑97/‑150)CJ019_03802/mcpU_3hypothetical protein/Methyl‑accepting chemotaxis protein McpU
Reads supporting (aligned to +/- strand):  ref base C (26/3);  new base A (0/10);  total (26/15)
Fisher's exact test for biased strand distribution p-value = 4.50e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCAGCATGATGCCGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAA  >  ACB143/4028295‑4028564
                                                                                                                                     |                                                                                                                                         
tCAGCATGATGCCGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACaaa                                                                                                                                        >  8:30757/1‑137 (MQ=255)
     atgatgCCGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTaa                                                                                                                                   >  7:237468/1‑137 (MQ=255)
     atgatgCCGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTaa                                                                                                                                   >  7:237470/1‑137 (MQ=255)
     atgatgCCGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGAGGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTaa                                                                                                                                   >  8:190509/1‑137 (MQ=255)
           ccGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑aa                                                                                                                            >  2:260889/1‑137 (MQ=255)
           ccGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑aa                                                                                                                            >  3:151798/1‑137 (MQ=255)
                 gcggcgATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAACGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑ATTGTTTg                                                                                                                      >  7:89185/1‑137 (MQ=255)
                 gcggcgATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCTCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTg                                                                                                                      >  7:154233/1‑137 (MQ=255)
                 gcggcgATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTg                                                                                                                      >  5:59331/1‑137 (MQ=255)
                                           cTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGa                                                                                            >  1:270632/1‑137 (MQ=255)
                                              aCCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGatat                                                                                         >  5:313341/1‑137 (MQ=255)
                                                 tCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATtgt                                                                                      >  4:30721/1‑137 (MQ=255)
                                                        gACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGc                                                                               >  1:105405/1‑137 (MQ=255)
                                                            gAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCg                                                                           >  7:248117/1‑137 (MQ=255)
                                                             aaaaaaGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGa                                                                          >  2:326483/1‑137 (MQ=255)
                                                             aaaaaaGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGACGATATTGTGGGAGGCAGCgc                                                                          >  8:71729/1‑136 (MQ=255)
                                                                aaaGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGt                                                                       >  2:64158/1‑137 (MQ=255)
                                                                aaaGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGAGGCGCTTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGTGAGCCAGCGAGGt                                                                       >  5:278996/1‑137 (MQ=255)
                                                                                 gccgggggCGGCATTTGTCAAGCGGGCCTGGTTGCGGCGGGTAACCGGCGGATAAATTAAAATTG‑ATTGTTTAAATATCGAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCtt                                                      <  8:270702/132‑1 (MQ=255)
                                                                                 gCAGTGGGCGGCATATGTCAAACGCGCCTTGAGGGGGCGCTTTACCGGCGCAAAAATTTAAATTG‑AATGTTTTAAAAAAAAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCtt                                                      <  4:249608/137‑1 (MQ=255)
                                                                                        gcggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGTCAGCGAGTTGCCTGGCCGGAGCGCTTTGGCAcg                                               >  1:188055/1‑136 (MQ=255)
                                                                                        gcggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCAcc                                               >  1:41127/1‑137 (MQ=255)
                                                                                        gcggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCAcc                                               >  7:274005/1‑137 (MQ=255)
                                                                                          ggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACctc                                             >  5:100245/1‑137 (MQ=255)
                                                                                                   aaaacGCGCCTTGATGCGGCGCTTTAACGGCGCAAAAATTTAAATTT‑AATGTTTAAAAAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCa                                    <  4:151798/136‑1 (MQ=255)
                                                                                                    aaaaGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCAc                                   >  6:36990/1‑137 (MQ=255)
                                                                                                          gcCTTTATGCGGCTGTTTAACGCCGTAAAAATTTAAATTT‑AATGTTTAAAAAAACAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGa                             <  5:249419/137‑1 (MQ=255)
                                                                                                          gcCTTGATGGGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTTAATAATGAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGa                             <  8:237470/137‑1 (MQ=255)
                                                                                                          gcCTTGATGCGGCGCGTTAACGGCGGAAAAATTTAAATTT‑AATGTTTGAATAAAGAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGa                             <  8:237468/137‑1 (MQ=255)
                                                                                                          cccTTGTGGCGGCGCTTTAACGGCGTAAAATTTAAAATTT‑ATTTTTTAAATAAAAAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGa                             <  5:174938/136‑1 (MQ=255)
                                                                                                          cccTTCATGCGGCGCTTAACCGCCGCAAAAATAAAAATTG‑AATGTTTAAAAAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGa                             <  1:79823/136‑1 (MQ=255)
                                                                                                               gATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAg                        >  1:86384/1‑137 (MQ=255)
                                                                                                                 tGCGGCGGTTAAAGGGCGTAAAAATTTAAATTT‑ATTTTTTTAATAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGcc                      <  1:25606/137‑1 (MQ=255)
                                                                                                                 tGCGGCGCTTTAACGGGGGAAAAATTTAATTTG‑AATGTTTGAAAAAATAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGcc                      <  5:263532/137‑1 (MQ=255)
                                                                                                                  gcggggCGTTACCGGCGCA‑AAATTTAAAATGCATTGTTTCAAAACAAAAACCACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCa                     <  4:123163/137‑1 (MQ=255)
                                                                                                                        ggTTAACGGCTCACAAATTAAAATTT‑AATGTTTAAATAAATAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAAccc               <  7:303738/136‑1 (MQ=255)
                                                                                                                        cgTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAAccc               >  5:132551/1‑137 (MQ=255)
                                                                                                                         tttAACGGCGGAAAAATTTAAATTG‑AATGTTTTAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCt              <  8:274005/136‑1 (MQ=255)
                                                                                                                         gTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGACATAACCCt              >  3:64873/1‑137 (MQ=255)
                                                                                                                             aCGGCGTAAAAATTTAAATTT‑ATTTTTTTAAAAAATAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTAtt          <  3:270908/137‑1 (MQ=255)
                                                                                                                                gcggACAAATTTAAATTT‑AATGTTTTAATAAAGAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGcc       <  6:92478/137‑1 (MQ=255)
                                                                                                                                     cAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCaaaaa  >  5:110258/1‑137 (MQ=255)
                                                                                                                                     |                                                                                                                                         
TCAGCATGATGCCGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTG‑AATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAA  >  ACB143/4028295‑4028564

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