Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,2130AT36.0% 55.4 / 12.3 51intergenic (‑356/+139)lysO/aqpZL‑lysine exporter/water channel AqpZ
Reads supporting (aligned to +/- strand):  ref base A (27/5);  new base T (0/18);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.86e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTC  >  NC_000913/915081‑915343
                                                                                                                                     |                                                                                                                                  
gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGgccgc                                                                                                                            >  2:167636/1‑141 (MQ=32)
 aTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGgccgcc                                                                                                                           >  1:664053/1‑141 (MQ=32)
   gtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGcccc                                                                                                                         >  1:1434897/1‑141 (MQ=32)
     gggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCt                                                                                                                       >  1:1670388/1‑141 (MQ=32)
      ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGcccccag                                                                                                                      >  1:328737/1‑139 (MQ=17)
      ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTg                                                                                                                      >  1:328676/1‑141 (MQ=32)
          gAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTgcgcg                                                                                                                  >  1:504034/1‑141 (MQ=32)
          gAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTgcgcg                                                                                                                  >  1:1843815/1‑141 (MQ=32)
          gAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTgcgcg                                                                                                                  >  1:1944857/1‑141 (MQ=32)
          gAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTgcgcg                                                                                                                  >  1:1944990/1‑141 (MQ=32)
                       accgggATGTAACCCCGCGCCGGACCCGGTAGCCTTACGGTTAACGTCGAGGCCCCCGCGCTGCGGCTATTTTGTACGCCAGTGCCCGGGGGTGCAGGGGGCGGCGTCG‑ATTGCCCCCCCCCTGCGCGCTCCTTGCGCCAg                                                                                                     <  2:952531/138‑1 (MQ=2)
                            gATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCa                                                                                                >  1:1551811/1‑141 (MQ=34)
                             ttGTTACCGCGCGCCCGACCGGGTAGCCTTACGGTATACGGGGTGGCCACCGCGCAGGGGCTATTTTGTTCGCCAGCGCCGGGGGGTGCGGGGGGCGGCGGCT‑TTGGCCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCaa                                                                                               <  2:1140952/140‑1 (MQ=255)
                               gTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAata                                                                                             >  1:1199739/1‑141 (MQ=37)
                               gTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAata                                                                                             >  2:1232450/1‑141 (MQ=37)
                                    cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                        >  2:1703672/1‑141 (MQ=40)
                                    cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                        >  2:273000/1‑141 (MQ=40)
                                                      gCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCCGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCCGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                      >  2:1016432/1‑141 (MQ=255)
                                                                  aaGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCAGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatg                                                          >  1:1215138/1‑141 (MQ=255)
                                                                         gggCTCCGCGCCGCGGCGATGGTGTTCGCCGGCGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTGGCGCCAGTGGCAATATGTTGCTTAGCTCCTGGCGGGAGCGCAACAAGAGGATCAATc                                                   >  2:1027486/1‑141 (MQ=17)
                                                                         gggCACCGCGCCGGGGGTTTTTTGTGCGCCAGCGCCGGGGGGTGCGGGGGGGGGCGGCG‑TTGGGCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATc                                                   <  2:1210470/141‑1 (MQ=18)
                                                                         gggCACCGCGCAGGGGGGTTTTTTTTCGCCAGAGCCGGGGGGGGCAGGGGGCGGCGGCT‑TTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATc                                                   <  2:1503304/141‑1 (MQ=18)
                                                                         gggCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCCCCACTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGAGGAATc                                                   >  2:1027490/1‑141 (MQ=255)
                                                                            cACCGCGCGGGGGCTTTTTTGGGCGCCATCGCGGGGGGGGGGGGGGGGGGGGGGCT‑TTGGCCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAAc                                                <  1:362135/141‑1 (MQ=255)
                                                                              ccgcgccgcgGCTTTTTTTTTCGCCAGAGCCGGGGGGTGCTGGGGGCGGGGGCT‑TTTGGCCCCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                              <  2:1944857/141‑1 (MQ=21)
                                                                              ccGCGCTGGGGCTTTTTTGTTCGCCTGAGCCCGGGGGTGCTGGGGGGGGCGGCT‑TTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                              <  2:504034/141‑1 (MQ=21)
                                                                              ccGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑CTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGGGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                              >  2:1070773/1‑141 (MQ=255)
                                                                              ccGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                              >  2:1070762/1‑141 (MQ=255)
                                                                              ccGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                              >  1:1421379/1‑141 (MQ=255)
                                                                                      ccgcgTTTGTGTTCGCCACCGCCCCGGGGTGCCGGGGGCGGCGCCT‑TTGGGCCGCCCCCCCCGCGCTCCTTGCGCCAGTGGCACTATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGCTGATGAATCAACATATaacaac                                      <  2:1183379/139‑1 (MQ=255)
                                                                                          ttttttGTACGCCAGAGCCCGGGGGTGCTGGGGGCGGGGGCT‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAACAACATATAACAACATCt                                  <  2:1484099/139‑1 (MQ=21)
                                                                                          gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCt                                  >  2:693227/1‑141 (MQ=255)
                                                                                          gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCAGCGCGCTCCTTGCGCCAGTGGCAATATTTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACCTATAACAACATCt                                  >  1:367870/1‑141 (MQ=255)
                                                                                          gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCAGCGCGCTCCTTGCGCCACTGGCAATATTTTTCTTAGCTCATGAAAGGAGCGCAACAAAATGATGTATCAACAAAAAacaacaaca                                  >  1:1589277/1‑138 (MQ=21)
                                                                                                       tgcgCCGGGGGGTGCGGGGGGCGGCGGCG‑TTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCt                     <  2:1551811/138‑1 (MQ=255)
                                                                                                       cgcgcgggggggTGCTGGGGGCGGCGGCG‑TTTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCt                     <  2:27609/134‑1 (MQ=37)
                                                                                                       agagCCGGGGGGTGCGGGGGGCGGCGGCT‑TTGGCCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCt                     <  2:1779027/141‑1 (MQ=37)
                                                                                                        gcgcgggggggTGCTGGGGGGGGGGGCG‑TTGGGCCCCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTg                    <  1:1941357/135‑1 (MQ=255)
                                                                                                        gcgccggggggggccgggggggGGGGCG‑TTGGCCCCCCCCCCCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTg                    <  2:712503/139‑1 (MQ=255)
                                                                                                        gcgCCGGGGGGTGCTGGGGGCGGCGGGT‑TTGGGCCCCCCCCGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTg                    <  2:1386361/139‑1 (MQ=255)
                                                                                                        gcgCCCGGCGGTCCAGGGGGCGCCGGCC‑ATTGGCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCACATGAAAGGAGCGCAACAAGCTGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTg                    <  2:583210/139‑1 (MQ=255)
                                                                                                        gagcgggggggTGGTGGGGGGGGGGGGG‑TTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTg                    <  1:459839/135‑1 (MQ=37)
                                                                                                        gagCCGGGGGGTGCTGGGGGGGGCGGCG‑ATTGGCCGCCCCATGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTg                    <  2:1386983/141‑1 (MQ=255)
                                                                                                            gggggggTGCTGGGGG‑GGGGGCGTTTTGCCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAt                <  1:1450153/139‑1 (MQ=255)
                                                                                                            cgggtggTGCTGGGGGCGGCGGCG‑TTTGGCCGCCCCGGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGCTGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAt                <  2:1100830/136‑1 (MQ=37)
                                                                                                            cggggggTGGTGGGGGGGGGGGCT‑TTTGGCCCCCCCCCGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAt                <  2:147101/139‑1 (MQ=37)
                                                                                                            ccGGGGGTGCTGGGGGCGGCGGCT‑ATTGGCCGCCCCGTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAt                <  2:315979/141‑1 (MQ=255)
                                                                                                                  gTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAATATGATGAAACAACATATAAAAACATATTAAAAAAAAGCCTTACAATAAGcc          >  1:2927/1‑141 (MQ=25)
                                                                                                                    gCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAGTATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAg        >  2:1836676/1‑141 (MQ=255)
                                                                                                                    gCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAACGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCAGACATTACGCCAg        >  2:990718/1‑141 (MQ=255)
                                                                                                                          ggggcgggcg‑TTTGCCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGACCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTc  <  2:159821/135‑1 (MQ=39)
                                                                                                                                     |                                                                                                                                  
GATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCG‑ATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTC  >  NC_000913/915081‑915343

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