Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,028,4510TA26.7% 71.1 / 10.6 45intergenic (‑120/‑127)CJ019_03802/mcpU_3hypothetical protein/Methyl‑accepting chemotaxis protein McpU
Reads supporting (aligned to +/- strand):  ref base T (25/8);  new base A (0/12);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 4.38e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGATCCGCCCCCATGCCTTTG  >  ACB143/4028315‑4028586
                                                                                                                                         |                                                                                                                                       
gcgATCACTGCATTCATGGTTGGCTAACCTCATAACGACATAAAAAAGCGCTAACTCTGAAGCAGTGGGCGGCAGATGTAAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTTAAt                                                                                                                                         <  2:294665/137‑1 (MQ=255)
  gATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATaa                                                                                                                                       >  4:217201/1‑137 (MQ=255)
  gATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATaa                                                                                                                                       >  6:201726/1‑137 (MQ=255)
  gATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATaa                                                                                                                                       >  2:66796/1‑137 (MQ=255)
      aCTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGaa                                                                                                                                   >  8:330351/1‑137 (MQ=255)
                   ttGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTATGAAGCAGTGGGCTGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAC‑CGGCGGACAAATTTAAATTGAAGGTTTGAATAAAGAAAAAACGGCAgggg                                                                                                                      >  2:344112/1‑137 (MQ=255)
                     ggCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGcc                                                                                                                    >  4:286718/1‑137 (MQ=255)
                           actctTTACACCAGAAAAAAGCCCGAACTCTGACGCAGGGGGGGTCATATTTAAAACGGGCGTTGATGCGGCGGGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAACACGGCAGGGGCCGATAtt                                                                                                              <  6:280041/136‑1 (MQ=255)
                              cATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATtgtg                                                                                                           >  8:226212/1‑137 (MQ=255)
                              cATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATtgtg                                                                                                           >  4:80945/1‑137 (MQ=255)
                                            aaaGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGt                                                                                             >  2:307015/1‑137 (MQ=255)
                                                     aCTCTCACGCCGTGGGCGGCATTTCTCAAACGCCCCTTGATGCGGCGCTTAAA‑CGGCGAACAATTTTAAATTTAATGTTTGAATAAACAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCg                                                                                    <  7:89583/137‑1 (MQ=255)
                                                              gtggggggggCATATGTAAAACGCGCCTGGATGGGGCGCTTAAA‑CGGCGGACAAATTTAAATTTAATTTTTCAATAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCttt                                                                           <  6:186937/133‑1 (MQ=255)
                                                              ccgggggCGGAAGTTTTAAAACGCGCGTTGATGCGGCG‑GTTAACCGCCTAAAAAATATAAATTTATTTTTTCAAAAAAAAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCttt                                                                           <  4:139066/133‑1 (MQ=255)
                                                               tgTGGGCGGCAGATGTAAACAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGAAAAATTTAAATTGAATGTTTCAATAAATAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTg                                                                          <  7:253100/136‑1 (MQ=255)
                                                                gggggCGGCATATTACAAAAGCGCCTTGATGCGGCGCTTAAA‑CGCCGTAAAAATTAAAATTTATTTTTTCAAAAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTgg                                                                         <  4:315162/135‑1 (MQ=255)
                                                                gggggCGGCATATGACAACCGCGCCTTGAGGCGGCGCTTAAA‑CGGCGCACAAATTTAATTTTAATGTTTTAAAAAAAAAACAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTgg                                                                         <  3:164141/135‑1 (MQ=255)
                                                                gggggCGGCAGTGGACAAACGCGCCTTGAGGCGGCGCGTTAA‑GGGCGCACAAATTAAAATTTAATGTTTCAATAAACAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTgg                                                                         <  4:309501/135‑1 (MQ=255)
                                                                    gcggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCCGCGAGGTGCCTGGCCGGAGCGCTTTGGCAcc                                                                     >  6:234853/1‑137 (MQ=255)
                                                                    gcggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCAcc                                                                     >  1:85718/1‑137 (MQ=255)
                                                                    gcggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCAcc                                                                     >  8:265788/1‑137 (MQ=255)
                                                                    gcggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCAcc                                                                     >  3:49347/1‑137 (MQ=255)
                                                                    gcggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCAcc                                                                     >  5:343907/1‑137 (MQ=255)
                                                                    gcggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCAcc                                                                     >  3:326038/1‑137 (MQ=255)
                                                                      ggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACctc                                                                   >  4:172502/1‑137 (MQ=255)
                                                                         aGATGCCAAAAGCGCCTTGATGCGGCGCTTTAA‑CGGCGCACAAATTTAAATTGAATGTTTTAAAAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCg                                                                <  6:195115/137‑1 (MQ=255)
                                                                         aGATGCAAAACGCGCCTTGATGCGGCGCTTAAA‑CGGCGGAAAAATTTAAATTTAATGTTTTAAAAAAGAAACAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCg                                                                <  8:162778/137‑1 (MQ=255)
                                                                           aTGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGc                                                              >  6:364821/1‑137 (MQ=255)
                                                                               cAAAAGCGCGTTGATGCGGCGGTTTAA‑CGCCGTAAAATTTTAAATTGACTTTTTCAAAACACAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCa                                                          <  4:183446/137‑1 (MQ=255)
                                                                               aaaaaGCGCCTTGAGGCGGCGCGTAAA‑CGGCGGACAAATTTAAATTGAATGTTTTAAAAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCa                                                          <  6:137061/136‑1 (MQ=255)
                                                                                     cgcCTTTATGGGGCGTGTTAC‑CGGCTAAAAAATTTAATTTTATTGTTTAAAAAAACAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAg                                                    <  3:217201/137‑1 (MQ=255)
                                                                                     cgcCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAg                                                    >  6:209461/1‑137 (MQ=255)
                                                                                        cTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATa                                                 >  1:132485/1‑137 (MQ=255)
                                                                                          tGATGCGGCGCGTTAC‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGATTTGGCACCTCTCGGCCCCACCGAAGATAcc                                               >  4:10237/1‑136 (MQ=255)
                                                                                          tGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACa                                               >  8:7583/1‑137 (MQ=255)
                                                                                          tGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATCGTGGGAGGCAGCTAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACa                                               >  4:84101/1‑137 (MQ=255)
                                                                                            ttGCGGCGTGTAAA‑GGGCGCAAAAATTTAATTTTATTGTTTCAAAAAACAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGc                                             <  1:58265/136‑1 (MQ=255)
                                                                                                     tttAA‑CGGCGCAAAATTTTAAATTTAATGTTTCAATAAAAAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCt                                    <  5:355981/136‑1 (MQ=255)
                                                                                                                      tttAATTTGAATTTTTAAAAAAAAAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGa                    <  5:336159/137‑1 (MQ=255)
                                                                                                                      tttAAATTTAATGTTTTAAAAAAGAAAACACGGCAGGGTCCGATATTGTGGGAGGCAGCGCGGGGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGa                    <  8:338380/137‑1 (MQ=255)
                                                                                                                      tttAAATTGAATGTTTCAAAAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGa                    <  8:242377/137‑1 (MQ=255)
                                                                                                                               aaTGTTTTAATAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGATCCGccccc           <  1:89632/137‑1 (MQ=255)
                                                                                                                                  gTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCTATAACCCTTATTGCCAAAAAGGGCTCCGCCCCCATg        >  3:8059/1‑137 (MQ=255)
                                                                                                                                  gTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGATCCGCCCCCATg        >  7:365191/1‑137 (MQ=255)
                                                                                                                                  gTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGATCCGCCCCCATg        >  1:95901/1‑137 (MQ=255)
                                                                                                                                        aTAAAGAAAAAACGGCAGGGGCAGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGATCCGCCCCCATGCCTTTg  <  5:364821/137‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
GCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAA‑CGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGATCCGCCCCCATGCCTTTG  >  ACB143/4028315‑4028586

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: ‑