Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9110GA35.3% 42.9 / 14.0 34intergenic (‑122/+1)fliA/fliXRNA polymerase sigma factor FliA/small RNA FliX
Reads supporting (aligned to +/- strand):  ref base G (16/6);  new base A (0/12);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 3.72e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCT  >  NC_000913/2001664‑2002155
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     
gTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaa                                                                                                                                                                                                                                                     >  2:427273/1‑250 (MQ=255)
     gTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTTTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAc                                                                                                                                                                                                                                               <  1:272367/251‑1 (MQ=255)
         cACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccg                                                                                                                                                                                                                                          >  1:121485/1‑251 (MQ=255)
           cGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgcc                                                                                                                                                                                                                                        <  1:362828/251‑1 (MQ=255)
                  aggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATATCGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTgg                                                                                                                                                                                                                                   >  2:47252/2‑249 (MQ=255)
                                           gcttcgggcccccccaggggaaCATAACGCTGCACCAGCGGTTTTTTTCCCATTACCCCTTCGGGGGTAAAGAGTGAATTCACAATAAACACCCCGGCTTTATATGAGTTATCGGCATTATTTCCCGTTTCGGCAGGGTTTTTAATCGGACTATTGGTGGGGAAAATGAGGGGTTATTGGGGGGTTAAAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                                                                                                                        <  1:245788/230‑1 (MQ=255)
                                              gcgggACGCCCCAGCGGGACATAACGTTGCCACAGGGATTTTTTACCCTTTCCCCTTTCAGGGGTATAGAGTGAATACACGAAAAACAGCCCTGTTTTATATGATTTACCGGCATGATTATCCTTTTCTGCGGGTTTTTAAAGCGGACGATTAGGGGGAGAAAGGGGGGTTTATTGGGGGGTTAGGGGTAAATTCAGGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc                                                                                                                                                                                                      <  1:321653/246‑1 (MQ=255)
                                                    cGCACCAGGGGGACATACCGCTGCCCCAGCGATTTTTTACCCATTCCCCCTTCGGGGTTATAGAGTAAATTCACGAAAAACAGCCCTGCTTTATATGAGTTATCGGAATGTTTTTCCTTTTCTGCAGGGTTTTTAATCGGACGATTGGGGGGTGAAATGGGGGTTTATTGGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                                                                                                               <  1:7369/251‑1 (MQ=255)
                                                     gCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                                                                                                              >  2:341042/1‑251 (MQ=255)
                                                              gaacataccggtgcaaagaggtttttttACCCATTACCCCTCCAGCGGTATAGGGTAAATTCACGATAAACAGCCCGGCTTTATATGTGTTATGGGAATGTTTATCCGTTTCGGCAGGGTTTTAAATCGGACTATTAGTGGGTGAAATGGGGGTTTTTTTGGGGGTTACAGGAAAATTCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                                                                                                                     <  2:97911/226‑1 (MQ=255)
                                                                           ccACAGCGAGGGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGGTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                                                                                                                        >  2:58856/1‑251 (MQ=255)
                                                                            cacaGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                                                                                                                       >  1:361517/1‑251 (MQ=255)
                                                                                 gaggttttttttcccttcctccTTCAGCGTTAGAGAGTGATTTCACAATAAACAGCCCTGCGTTATATTATTTATCGGCATGATTATCGTTTGTTGCGGGGTTTTTAATGGGAGGATGGGGGGGGGAAAGGGGGGTTTTTTGGGGGGTAACTGTTAAATCCCAGGCAAAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                                                   <  2:473693/231‑1 (MQ=255)
                                                                                      tttttcccATAACACCTTCGGCGGTGAAGAGTCAATTCAGGATAAACAGCCTGGGGTTATTTGGGTTAGCGGTATGATTCTCCGTTCCTGAAGGGTTTTAATTCGAAGGATTGGTGGGGAAATTGGGGGTTTTTTGGGGGGTTACAGTAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                                                                                              <  1:369329/244‑1 (MQ=255)
                                                                                                        gcgggggtttagagtaaagacaggATAAACAGCCTTGTGTTTTATGATTTAGGGGCATGTTTATCGGTTTCTGAGGGTTTTTTAATGGGACGGTTAGTGGGTAAAAGGGGGGTTTTTTTGGGGTTTAGAAAAAAATTCCCGCCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                                                                                                                            <  2:93291/228‑1 (MQ=255)
                                                                                                                   gagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGGAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATTTGCAAATTCAATAAATTGc                                                                                                                                  >  2:402959/1‑250 (MQ=255)
                                                                                                                       gatttCACGAAAAACGGCCCTGGTTTATATGATTTATCGGCATGTTTATCCTTTTCTGCAGGGTTTTTAATCGGACGTTTAGTGGGTGAAATGAGGGGTTTTTGGGGGGTTAAAGTAAAATTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                                                                                                              <  1:261562/247‑1 (MQ=255)
                                                                                                                       aaaTTCACAAAAAACAGCCCGGGGTTAAATGAGTTATCGGCATGTTTGTCCGTTTCGGAGGGTTTTTAAATGGGAGGTTTAGGGGGTGAAATGGGGGGTTATTGGGGGTTTACAGGAAAATCCCAGACGAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                                                                                                             <  1:141459/250‑1 (MQ=255)
                                                                                                                            cACGATAAACAGCCCTGCTTTATAGGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGAGAAATGAGGGTTTTTTGGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt                                                                                                                         <  1:625093/250‑1 (MQ=255)
                                                                                                                                     cAGCCCTGCGTTATATGCGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGGGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                                                 >  2:172066/1‑201 (MQ=255)
                                                                                                                                          cTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                                                                                                            >  2:464890/1‑249 (MQ=255)
                                                                                                                                                  tatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc                                                                                                  >  1:617917/1‑251 (MQ=255)
                                                                                                                                                     gAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc                                                                                               >  1:623709/1‑251 (MQ=255)
                                                                                                                                                     gAGTTATCGGCATGATTATCCGTTTCTGCAGGGGTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc                                                                                               >  2:198527/1‑251 (MQ=255)
                                                                                                                                                         tATCGGCATGATTATCCTTTTCTGCGGGTTTTTTAATCGGAGAATTGGTGGGTAAAATGAGGGTTTTTTTGGGGTTTACAGAAAAATTCCAGACAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgn                                                                                            <  2:285839/250‑2 (MQ=37)
                                                                                                                                                           tCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                                          >  2:297372/1‑250 (MQ=255)
                                                                                                                                                               cATGATTATCCGATTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCCGCGCATCAGGCAATTTg                                                                                     >  2:554571/1‑251 (MQ=255)
                                                                                                                                                                     tttccttttCTGAGGGTTTTTTAATCGGACGTTTAGTGGGTGAAAGGGGGGTTTTTTGGGGGGTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGtt                                                                                <  1:364587/244‑1 (MQ=25)
                                                                                                                                                                          gTTTCTGCAGGGTTTTTAATCGGACGATTAGGGGGTGAACTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTAGCCGGATGCGGCGTAAACTCCTTATCCGGCCTACAAAAATGCGCAAATTCAATAAATTGCAATTCCACTTGTCGGCCTGATAAGCGCAGTGCATCAGGCAATTAGGCGTTGCCGTc                                                                          >  2:452216/1‑251 (MQ=255)
                                                                                                                                                                             tttggCAGGTTTTTTAAGCGGACGATGAGTGGGTAAAAGGAGGGTTTATTTGGGGATGAAAGGTAAATCCCAGGAAAAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCa                                                                         <  2:95553/246‑1 (MQ=37)
                                                                                                                                                                                 ggaggttttttAAGCGGATGATTGGTGGGTAAAAGGGGGTGTTTTTTGGGTGTTACGGGAAATTTCGAGGAAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGtctc                                                                    <  2:190789/245‑1 (MQ=25)
                                                                                                                                                                                                                 gAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGGTGCTGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCACTTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAAccccgc                                      >  1:159732/1‑248 (MQ=255)
                                                                                                                                                                                                                            tttttggggggTTACAGGAAAATTCCGGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACa                            <  1:427273/244‑1 (MQ=255)
                                                                                                                                                                                                                                                      cAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCt  <  1:341042/251‑1 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     
GTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCT  >  NC_000913/2001664‑2002155

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: