Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,014,0230TG47.4% 56.1 / 10.3 39G255G (GGT→GGG) pqiBparaquat‑inducible, SoxRS‑regulated MCE domain protein
Reads supporting (aligned to +/- strand):  ref base T (1/19);  new base G (18/0);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 1.13e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGG  >  NC_000913/1013874‑1014162
                                                                                                                                                     |                                                                                                                                            
gCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGt                                                                                                                                              <  1:903365/150‑1 (MQ=255)
         gACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAg                                                                                                                                   <  1:26331/151‑1 (MQ=255)
          aCCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGGGCTGGA‑ttt                                                                                                                                     >  2:257117/1‑148 (MQ=255)
              aCTGGTGACCAACAATGTTCGCTGCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAg                                                                                                                              <  2:501304/151‑1 (MQ=255)
               cTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGGGGCAGCTTTTATGTGGCGGAAGGGCTGGA‑TTTACGGc                                                                                                                                >  1:360175/1‑148 (MQ=255)
                         acaaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCAcc                                                                                                                   <  2:218743/151‑1 (MQ=255)
                             tGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGCGGCCGCTTTGATGGGCCGGAAGGACGGGG‑TTTAGGGCAGCCAGGGGCACCGaa                                                                                                                >  1:780950/1‑150 (MQ=255)
                              gTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGGCGGGAGGGCTGGA‑TTTAGGGCAGCCAGTGGGACCGaaa                                                                                                               >  1:875842/1‑150 (MQ=255)
                              gTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGaaaa                                                                                                              <  2:494735/151‑1 (MQ=255)
                                    ttCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGGGCTGGA‑TTTAGGGCAGCCAGGGGCACCGAAAAAAGCtt                                                                                                        >  2:322414/1‑151 (MQ=255)
                                     tCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGTTGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCATTGTATTGCCGGACAGGCTGGG‑TTTTGGGTAGCCTGGGGGACCGaaaaaaacgt                                                                                                        >  1:375813/1‑144 (MQ=255)
                                      cTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGAAGGGGTGGG‑GGTTGGGGAGCCAGTGGCACCGAAAAAAACtt                                                                                                        >  2:462453/1‑149 (MQ=255)
                                         gAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGGCGGAAGGGCTGGG‑TTTAGGGGAGCCCGTGGCACCGAAAAAAACTTTTGtt                                                                                                   >  1:881714/1‑151 (MQ=255)
                                           aaGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGGTTTGTTGTGCCGGGAGGGCTGGG‑TTTTGGGGAGCCAGTGGCACCGAAAACAGCCTTCGtttt                                                                                                 >  2:497152/1‑151 (MQ=255)
                                           aaGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGtttt                                                                                                 <  1:579039/151‑1 (MQ=255)
                                              aTAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTAAGCCTTGTGGTGGCGGGAAGGCTGGG‑GTTTGGGCGGGCAGGGGCCCCCGAAACAGCTtttgttttgt                                                                                               >  2:105686/1‑150 (MQ=255)
                                                              tGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCCGCTTTGATTTGGCGGAAGGGCTGGG‑TTTTGGGGAGCCAGGGGGACCGCAAACAGCTTTCGGTTTGTTTGATGATaagaag                                                                                 >  2:724130/1‑148 (MQ=255)
                                                               gATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGc                                                                              <  1:415916/150‑1 (MQ=255)
                                                                       gTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGAAGGGCTGGG‑TTTTGGGCAGCCAGGGGCACCGAAAACAGCTTTCGTTTTTTTTGATGATCAGAAGAGCATTTAgg                                                                       >  1:817357/1‑149 (MQ=255)
                                                                       gTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGAt                                                                     <  2:415704/151‑1 (MQ=255)
                                                                        tCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCCGCTTTTATGTGGCGGGAGGGCTGGG‑TTTTGGGGAGCCCGTGGGACCGGCAAAAACCTTTGTTTTTTTTGATGATCAAAaaaggaattagggtt                                                                    >  2:847285/1‑140 (MQ=255)
                                                                              gggATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGt                                                              <  1:952573/151‑1 (MQ=255)
                                                                               ggATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTCTCCGATTTTAGTTGGGGGGAGGGCCGGAGTTTTGGGGAGCC‑GGGGGACCGGAAACAACTTTCGGTTTGTTTTTTGGTTCGAAAAGCCGTCAGGAGTCGtt                                                                >  2:575298/1‑148 (MQ=255)
                                                                                       gtgGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGAAGGGCTGGG‑TTTAGGGGAGCCAGTGGCACCCAAAACAGCTTTTGTTTTGTATGATGATCAGAAGAGCATTCAGGAATCGTTGTACACCgct                                                      >  1:260052/1‑148 (MQ=255)
                                                                                               gggCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGa                                              <  2:115499/150‑1 (MQ=255)
                                                                                               gggCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGAt                                             <  1:220596/151‑1 (MQ=255)
                                                                                                  cTCATTGACAACGCTGCTGAGTGGCGGTGGCAGCTTTGTTGTGCCGGAAGGGCTGGG‑TTTTGGGGAGGCAGTGGGACCGGAAAACGCTTTCGTTTTGTATTGTGATTAGAAGAGGATTCAAGATTCGTTGTTCACCGATTACATTGAt                                             >  2:705421/1‑148 (MQ=255)
                                                                                                   tCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGGCGGGAGGGCTGGA‑TTTTTGGGAGCCCGGGGGACCGAAAAAAGCTTTCGTTTTGTGTGATGATTAGAAAAGCATTCCGGATTCGTTGTACACCGATCACAtttttt                                            >  2:187539/1‑144 (MQ=255)
                                                                                                   tCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGGCTGGG‑TTTAGGGGAGCCAGTGGGACCGAAAACAGCTTTTGTTTTTTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACAACGATCACATTGAtt                                            >  2:492459/1‑148 (MQ=255)
                                                                                                          cAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGGGCTGGG‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGtt                                  >  1:800331/1‑151 (MQ=255)
                                                                                                             cGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCtt                               <  2:505999/151‑1 (MQ=255)
                                                                                                              gctgctGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCttt                              <  1:588156/151‑1 (MQ=255)
                                                                                                              gctgctGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCttt                              <  2:759643/151‑1 (MQ=255)
                                                                                                                    gAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGAt                        <  2:87576/151‑1 (MQ=255)
                                                                                                                      gTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTc                      <  1:518104/151‑1 (MQ=255)
                                                                                                                              gtCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTAcgc               <  1:226620/150‑1 (MQ=255)
                                                                                                                               tCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTAcgcg              <  2:609436/150‑1 (MQ=255)
                                                                                                                                  gCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCTATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCt          >  1:742007/1‑151 (MQ=255)
                                                                                                                                        aTGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAAcc    <  2:380890/151‑1 (MQ=255)
                                                                                                                                          gtgCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCgg  <  1:422008/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGA‑TTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGG  >  NC_000913/1013874‑1014162

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑