Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,039,7930CA34.4% 41.6 / 12.5 32intergenic (‑97/‑150)PP_3556/PP_3557membrane protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (16/5);  new base A (0/11);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 6.77e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCAGCATGATGCCGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATT  >  NC_002947/4039660‑4039921
                                                                                                                                     |                                                                                                                                
tCAGCATGATGCCGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACaa                                                                                                                                >  1:110302/1‑136 (MQ=255)
     atgatgCCGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTa                                                                                                                           >  1:558348/1‑136 (MQ=255)
           ccGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGa                                                                                                                     >  1:441342/1‑136 (MQ=255)
           ccGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGa                                                                                                                     >  1:374775/1‑136 (MQ=255)
           ccGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGa                                                                                                                     >  1:287597/1‑136 (MQ=255)
           ccGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGa                                                                                                                     >  1:347326/1‑136 (MQ=255)
                 gcggcgATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGttt                                                                                                               >  1:353939/1‑136 (MQ=255)
                                      gtgggCTTCCCTCATTCCGACAGAAAAACGCGCCAACTCTGACGCAGTGGGCGGCAGAGGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAgg                                                                                          <  2:386611/133‑1 (MQ=255)
                                      gTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAgg                                                                                          >  2:353350/1‑136 (MQ=255)
                                         ggCTTACCTCATTACCACATAAAAACGCGCGACCTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCTTTAACGGCGGACAAATTAAAATTGAATGTTTCAAAAAACAAAAAACGGCAGGGGc                                                                                       <  2:11323/136‑1 (MQ=255)
                                                  aaTTACCACATAAAAAAGCGCCAACTCTGAAGCAGTGGGCGGCAGATGTCAAACGCGCCTTGATGCGGCGCGTTAACGGCGCACAAATTTAAATTGAATGTTTTAATAAACAAAAAACGGCAGGGGCCGATATtgt                                                                              <  2:287597/135‑1 (MQ=255)
                                                                aaaGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAgg                                                                >  1:218871/1‑136 (MQ=255)
                                                                aaaGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAgg                                                                >  1:76364/1‑136 (MQ=255)
                                                                          ctctGAAGCAGGGGGCGGCAGATGTAAAAAGCCCGTTGATGCGGCGCTTAAACGGCGCAAAAATTAAATTTTAATGTTTCAAAAAACAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCg                                                      <  2:374775/136‑1 (MQ=255)
                                                                                  ttgTGGGCGGCATTTGTAAAACGCGCCTTGATGCGGCGCTTTAACGGCGTAAAAATTTAATTTTATTGTTTTAAAAAATAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCtt                                              <  2:441342/134‑1 (MQ=255)
                                                                                  ttgTGGGCGGCATATGAAAAATGCGCCTTTTTGCGGCGCTTTAACGGCTTAAAAATTTAAATTTAATTTTTCAATAAAAAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCtt                                              <  2:347326/134‑1 (MQ=255)
                                                                                   tgggggCGCCATATGTCAAACGCCCCTTGTTGCGGCGCGTAAACGGCGCAAAAATTTAAATTAAATTTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCttt                                             <  1:252205/133‑1 (MQ=255)
                                                                                    gTGGGCGTCATATGTAAAAAGCTCGTTCTTGCGGCGCTATAACGGCTTAAAAATTTAAATTTAATTTTTCAATAAACAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTg                                            <  2:193442/136‑1 (MQ=255)
                                                                                        gcggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCAc                                        >  2:49070/1‑136 (MQ=255)
                                                                                        gcggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCAc                                        >  1:292407/1‑136 (MQ=255)
                                                                                          ggcAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACct                                      >  1:101318/1‑136 (MQ=255)
                                                                                                    aaacGCGCCTTGATGCGGCGCTTAAACGCCGTAAAAATTAAAATTTAATGTTTTAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCa                            <  2:218871/136‑1 (MQ=255)
                                                                                                    aaaaGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGccccc                            >  2:195537/1‑135 (MQ=255)
                                                                                                    aaaaGCGCCTTGAGGGGGCGCTTAACCGGCGCAAAAATTAAAATTTATTTTTTTAAAAAAAAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCa                            <  1:353350/136‑1 (MQ=255)
                                                                                                       aGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCg                         >  2:370635/1‑136 (MQ=255)
                                                                                                          gcCTTGAGGCGGCGCTTTAACGGCGCAAAATTTAAAATTAAATTTTTCAAAAAACAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAg                      <  1:118614/136‑1 (MQ=255)
                                                                                                                 tGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGc               <  2:131563/136‑1 (MQ=255)
                                                                                                                 ggcggcgTTTTAACGGCTAAAAAATTTAAATTATTTTTTTTAATAAAAAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGc               <  1:495073/135‑1 (MQ=255)
                                                                                                                         ttaaaCGGCTAAAAAATTTAAATTTATTTTTTAAAAAAAAAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAAccc       <  1:45070/133‑1 (MQ=255)
                                                                                                                         gTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAAccc       >  1:169713/1‑136 (MQ=255)
                                                                                                                          ttAACGGCGTACAAATTTAAATTGAATGTTTGAATAAATAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCt      <  2:143151/136‑1 (MQ=255)
                                                                                                                              gggcgtAAAAATTAAATTTTAATGTTTTAATAAAAAAAACACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTAtt  <  1:427809/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                
TCAGCATGATGCCGACCGCGGCGATCACTGCATTCATGGTTGGCTTACCTCATTACGACAGAAAAAAGCGCGAACTCTGAAGCAGTGGGCGGCAGATGTCAAAAGCGCCTTGATGCGGCGCGTTAACGGCGGACAAATTTAAATTGAATGTTTGAATAAAGAAAAAACGGCAGGGGCCGATATTGTGGGAGGCAGCGAGGTGCCTGGCCGGAGCGCTTTGGCACCTCTCGGCCCCACCGCAGATACAGCCATAACCCTTATT  >  NC_002947/4039660‑4039921

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