Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,028,4550GC33.3% 40.8 / 14.0 36intergenic (‑124/‑123)CJ019_03802/mcpU_3hypothetical protein/Methyl‑accepting chemotaxis protein McpU
Reads supporting (aligned to +/- strand):  ref base G (16/2);  new base C (0/9);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 1.17e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGATCCGCCCCCATGCCTTTG  >  ACB143/4028326‑4028586
                                                                                                                                 |                                                                                                                                   
aTTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGaaaaaa                                                                                                                               >  5:91861/1‑136 (MQ=255)
       gTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAgg                                                                                                                        >  3:183030/1‑136 (MQ=255)
          ggCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGc                                                                                                                     >  5:275205/1‑136 (MQ=255)
               aCCTCTTTACTACAGAAAAATGCGCTAACTGTCATCCAGTGGGGGGCAGATGAAAAACGCGCGTTGTTGCGGGGGTTTAAGGGCGAAAAAATTTAAATTTAATTTTTCATTAAATAAAAAACGGCAGGGGCCg                                                                                                                   <  6:253464/133‑1 (MQ=255)
                 cTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGatatc                                                                                                              >  7:302135/1‑135 (MQ=255)
                 cTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGatatc                                                                                                              >  5:180178/1‑135 (MQ=255)
                         gACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAgg                                                                                                      >  3:114733/1‑136 (MQ=255)
                          aCAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAgg                                                                                                      >  6:277716/1‑135 (MQ=255)
                              aaaaaaGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCg                                                                                                 >  1:157654/1‑136 (MQ=255)
                                 aaaGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAgg                                                                                              >  5:72663/1‑136 (MQ=255)
                                                    agggggCGGCAGATGAAAAACGCGCCTTGATGCGGCGCTTAAACGGCGCACAATTTAAAATTTAATGTTTCAAAAAAAAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCttt                                                                           <  8:155956/133‑1 (MQ=255)
                                                     ttGGGGGGCAGATGACAAACGCGCCTTGTTGCGGCGCTTTAACGGCGTAAAATTTTAAATTGAATGTTTCAATAAATAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTg                                                                          <  8:24558/135‑1 (MQ=255)
                                                     gggggCGGCATATGAAAAACGCGCCTTGTGCCGGCGCTTTAACGCCTCAAAAATTAAAATTTAATTTTACAAAAAAAAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTg                                                                          <  4:136129/134‑1 (MQ=255)
                                                      ggggCGGCTCTTGACAAACGCGCCTTGAGGCGGCGCTTAAACGGGCCAAAATTTTAAATTTATTTTTTAAAAAAAAAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTgg                                                                         <  6:78544/135‑1 (MQ=255)
                                                         gcggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCAc                                                                      >  1:262162/1‑136 (MQ=255)
                                                         gcggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCAc                                                                      >  1:336615/1‑136 (MQ=255)
                                                         gcggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCAc                                                                      >  4:145555/1‑136 (MQ=255)
                                                                 tGTCAAATGCGCCTGCCGGCGGCGCTTTACCGGCGGACAAATTTAAATTGAATTTTTTAAAAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCgg                                                               <  6:96687/135‑1 (MQ=255)
                                                                 tGTCAAAAGCGCCTTGATGCGGCGCTTTAACGGCGGACAAATTTAAATTTAATGTTTGAAAAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCgg                                                               <  5:268076/135‑1 (MQ=255)
                                                                     aaacccgcCTTGAGGCGGCGCTTAACCGGCGCAAAAATTAAAATTGATTTTTTAAAAAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCa                                                          <  3:330414/131‑1 (MQ=255)
                                                                     aaacGCGCCGTGAGGGGGCGCGTAACCGGCGCACAATTTAAATTTTAATTTTTCAAAAAACAAAACCCGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCa                                                          <  6:193168/136‑1 (MQ=255)
                                                                     aaaaGCGCGTTGATGCGGCGCTTTAACGGCGTAAAATTTAAAATTCATTTTTTCAAAAAACAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCa                                                          <  3:317465/136‑1 (MQ=255)
                                                                          cgcCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCg                                                       >  7:243652/1‑134 (MQ=255)
                                                                           gcCTTTATGCGGCGGTTTAACGGCGCAAAATTTTAAATTGAATGTTTCAATAAATAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAg                                                    <  8:187551/136‑1 (MQ=255)
                                                                           gcCTTGCTGCGGCGGTTAAAGGGCGCAAAAATTAAATTTTATTGTTTAAAAACAAAAAACCCGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAg                                                    <  5:127118/136‑1 (MQ=255)
                                                                           gcCTTGATGCGGCGGTTTAACGGCGTAAAAATTTAAATTCAATGTTTTAAAAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAg                                                    <  2:4485/136‑1 (MQ=255)
                                                                            ccTTGAGGCGGCGGTTAAAGGGCGTAAAAATTTAATTTTAATGTTTCAATAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGa                                                   <  6:229400/136‑1 (MQ=255)
                                                                               tGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCCCCGCAGAtcc                                                >  2:325035/1‑134 (MQ=255)
                                                                                   gcggcgCGTTAACGCCGCAAAAATTAAAATTTATTGTTTCAAAAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGc                                             <  6:80085/135‑1 (MQ=255)
                                                                                          gtaaaCGGCGTAAAAATTTAAATTCAATTTTTTAAAAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAAccc                                     <  4:321807/133‑1 (MQ=255)
                                                                                                     aCAATTTAAAATTTATTTTTTTAAAAAAAAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCaa                           <  3:81079/135‑1 (MQ=255)
                                                                                                       aaaTTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCaaaag                        >  6:196925/1‑135 (MQ=255)
                                                                                                            tAAATTTATTTTTTTAATAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGa                    <  7:315792/135‑1 (MQ=255)
                                                                                                                   aaTGTTTCAATAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGATCCGcccc            <  6:309119/136‑1 (MQ=255)
                                                                                                                      gTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGATCCGCCCCCAt         >  7:168283/1‑136 (MQ=255)
                                                                                                                       tataaaTAAAAAAACCACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGATCCGCCCCCAt         <  8:71404/131‑1 (MQ=255)
                                                                                                                             aaaaaaaaaCACGGCCGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGATCCGCCCCCATGCCTTTg  <  1:88642/135‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                   
ATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATTGCCAAAAAGGGATCCGCCCCCATGCCTTTG  >  ACB143/4028326‑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: ‑