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→GGGpqiBparaquat‑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

GCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGG  >  NC_000913/1013874‑1014162
                                                                                                                                                     |                                                                                                                                            
gCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGt                                                                                                                                              <  1:903365/150‑1 (MQ=255)
         gACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAg                                                                                                                                   <  1:26331/151‑1 (MQ=255)
          aCCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGGGCTGGAttt                                                                                                                                     >  2:257117/1‑148 (MQ=255)
              aCTGGTGACCAACAATGTTCGCTGCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAg                                                                                                                              <  2:501304/151‑1 (MQ=255)
               cTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGGGGCAGCTTTTATGTGGCGGAAGGGCTGGATTTACGGc                                                                                                                                >  1:360175/1‑148 (MQ=255)
                         acaaTGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCAcc                                                                                                                   <  2:218743/151‑1 (MQ=255)
                             tGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGCGGCCGCTTTGATGGGCCGGAAGGACGGGGTTTAGGGCAGCCAGGGGCACCGaa                                                                                                                >  1:780950/1‑150 (MQ=255)
                              gTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGGCGGGAGGGCTGGATTTAGGGCAGCCAGTGGGACCGaaa                                                                                                               >  1:875842/1‑150 (MQ=255)
                              gTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGaaaa                                                                                                              <  2:494735/151‑1 (MQ=255)
                                    ttCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGGGCTGGATTTAGGGCAGCCAGGGGCACCGAAAAAAGCtt                                                                                                        >  2:322414/1‑151 (MQ=255)
                                     tCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGTTGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCATTGTATTGCCGGACAGGCTGGGTTTTGGGTAGCCTGGGGGACCGaaaaaaacgt                                                                                                        >  1:375813/1‑144 (MQ=255)
                                      cTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGAAGGGGTGGGGGTTGGGGAGCCAGTGGCACCGAAAAAAACtt                                                                                                        >  2:462453/1‑149 (MQ=255)
                                         gAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGGCGGAAGGGCTGGGTTTAGGGGAGCCCGTGGCACCGAAAAAAACTTTTGtt                                                                                                   >  1:881714/1‑151 (MQ=255)
                                           aaGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGGTTTGTTGTGCCGGGAGGGCTGGGTTTTGGGGAGCCAGTGGCACCGAAAACAGCCTTCGtttt                                                                                                 >  2:497152/1‑151 (MQ=255)
                                           aaGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGtttt                                                                                                 <  1:579039/151‑1 (MQ=255)
                                              aTAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTAAGCCTTGTGGTGGCGGGAAGGCTGGGGTTTGGGCGGGCAGGGGCCCCCGAAACAGCTtttgttttgt                                                                                               >  2:105686/1‑150 (MQ=255)
                                                              tGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCCGCTTTGATTTGGCGGAAGGGCTGGGTTTTGGGGAGCCAGGGGGACCGCAAACAGCTTTCGGTTTGTTTGATGATaagaag                                                                                 >  2:724130/1‑148 (MQ=255)
                                                               gATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGc                                                                              <  1:415916/150‑1 (MQ=255)
                                                                       gTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGAAGGGCTGGGTTTTGGGCAGCCAGGGGCACCGAAAACAGCTTTCGTTTTTTTTGATGATCAGAAGAGCATTTAgg                                                                       >  1:817357/1‑149 (MQ=255)
                                                                       gTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGAt                                                                     <  2:415704/151‑1 (MQ=255)
                                                                        tCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCCGCTTTTATGTGGCGGGAGGGCTGGGTTTTGGGGAGCCCGTGGGACCGGCAAAAACCTTTGTTTTTTTTGATGATCAAAaaaggaattagggtt                                                                    >  2:847285/1‑140 (MQ=255)
                                                                              gggATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGt                                                              <  1:952573/151‑1 (MQ=255)
                                                                               ggATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTCTCCGATTTTAGTTGGGGGGAGGGCCGGAGTTTTGGGGAGCCGGGGGACCGGAAACAACTTTCGGTTTGTTTTTTGGTTCGAAAAGCCGTCAGGAGTCGtt                                                                >  2:575298/1‑148 (MQ=255)
                                                                                       gtgGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGAAGGGCTGGGTTTAGGGGAGCCAGTGGCACCCAAAACAGCTTTTGTTTTGTATGATGATCAGAAGAGCATTCAGGAATCGTTGTACACCgct                                                      >  1:260052/1‑148 (MQ=255)
                                                                                               gggCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGa                                              <  2:115499/150‑1 (MQ=255)
                                                                                               gggCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGAt                                             <  1:220596/151‑1 (MQ=255)
                                                                                                  cTCATTGACAACGCTGCTGAGTGGCGGTGGCAGCTTTGTTGTGCCGGAAGGGCTGGGTTTTGGGGAGGCAGTGGGACCGGAAAACGCTTTCGTTTTGTATTGTGATTAGAAGAGGATTCAAGATTCGTTGTTCACCGATTACATTGAt                                             >  2:705421/1‑148 (MQ=255)
                                                                                                   tCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGGCGGGAGGGCTGGATTTTTGGGAGCCCGGGGGACCGAAAAAAGCTTTCGTTTTGTGTGATGATTAGAAAAGCATTCCGGATTCGTTGTACACCGATCACAtttttt                                            >  2:187539/1‑144 (MQ=255)
                                                                                                   tCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGGCTGGGTTTAGGGGAGCCAGTGGGACCGAAAACAGCTTTTGTTTTTTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACAACGATCACATTGAtt                                            >  2:492459/1‑148 (MQ=255)
                                                                                                          cAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGGGCTGGGTTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGtt                                  >  1:800331/1‑151 (MQ=255)
                                                                                                             cGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCtt                               <  2:505999/151‑1 (MQ=255)
                                                                                                              gctgctGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCttt                              <  1:588156/151‑1 (MQ=255)
                                                                                                              gctgctGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCttt                              <  2:759643/151‑1 (MQ=255)
                                                                                                                    gAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGAt                        <  2:87576/151‑1 (MQ=255)
                                                                                                                      gTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTc                      <  1:518104/151‑1 (MQ=255)
                                                                                                                              gtCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTAcgc               <  1:226620/150‑1 (MQ=255)
                                                                                                                               tCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTAcgcg              <  2:609436/150‑1 (MQ=255)
                                                                                                                                  gCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCTATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCt          >  1:742007/1‑151 (MQ=255)
                                                                                                                                        aTGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAAcc    <  2:380890/151‑1 (MQ=255)
                                                                                                                                          gtgCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCgg  <  1:422008/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGG  >  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: