Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913300,7110GC40.5% 68.1 / 20.2 42P75R (CCG→CGG) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base G (5/20);  new base C (17/0);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.08e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTC  >  NC_000913/300565‑300859
                                                                                                                                                  |                                                                                                                                                     
gcgcCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGgcgc                                                                                                                                                   <  1:716852/151‑1 (MQ=255)
  gcCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCtt                                                                                                                                                 <  1:610509/151‑1 (MQ=255)
        tCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCCCTACCGCCCGCGCtttttt                                                                                                                                             >  2:303104/1‑148 (MQ=255)
          gAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATTACAACCCGAACGCCCCGCGCctttt                                                                                                                                              >  2:946356/1‑146 (MQ=255)
                   ggCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAAGAAAGCCCGGACGCCCCGCGCTTTTTTCgcgggggcc                                                                                                                                   >  1:260521/1‑143 (MQ=255)
                    gCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGAAAACGAGTACGCCCCGCGCCCTTTTCGCGGGGGCCGGAt                                                                                                                               >  2:352240/1‑151 (MQ=255)
                       aCTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGTATTACAACCCGGACGCCCCGCGCCTTTTTCgcgggggcc                                                                                                                                   >  1:680372/1‑139 (MQ=255)
                        cTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑aa                                                                                                                           <  1:904393/150‑1 (MQ=255)
                            cAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGGAAGCCCCGCGGCTTTTTTGCGGGGGCCGTATC‑AAagggg                                                                                                                      >  1:6446/1‑149 (MQ=255)
                            cAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCg                                                                                                                      <  1:647686/151‑1 (MQ=255)
                                        cTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGGACGCCCCGCGCTTTTTTTGCGGCGTCCGTATC‑AAGGGCGGGGGGGCGGcc                                                                                                           >  1:123883/1‑150 (MQ=255)
                                         tgatgaTAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAACCAGTACGCCCCGCGCTTTTTTTGCGGGGGCCCGATC‑AAGGGGGGGGGGGGGTcc                                                                                                           >  1:515919/1‑149 (MQ=255)
                                          gatgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAATAAAAACAATACGCCCCGGGCCTTTTTGGCGGGGGCCCGATCAAAGGGGGGGGGGGGGTCCtt                                                                                                         >  1:535235/1‑151 (MQ=255)
                                           atgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTg                                                                                                       <  2:748592/151‑1 (MQ=255)
                                            tgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATTACAGCCCGTACGCCCCGCGCCTTTTTTGCGGCGCGCGGATC‑AAAGGCGGGGGGGCGGcc                                                                                                           >  2:105424/1‑146 (MQ=255)
                                            tgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGACAACCACAACGCCCCGCGCTTTTTTCGCGGCGGCCCTATC‑CAGGGCGGGGAGGGGGCCttt                                                                                                        >  1:618570/1‑149 (MQ=255)
                                              atAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAGCCAGTACGCCCCGCGCCTTTTTCGCGGCGGCCCGATC‑AAGGGGGGGGGGGGGGCCttttg                                                                                                      >  1:501377/1‑147 (MQ=255)
                                                    cTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAATACAGCCCGTACCCCCCGCGCGTTTTTTGCGGCGGCCGGATT‑CAGGGGGGGGGGGCGGCCTTTTGTAATGGCg                                                                                              >  2:998686/1‑151 (MQ=255)
                                                          aGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcg                                                                                        <  1:155497/151‑1 (MQ=255)
                                                              gggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgacga                                                                                    <  1:365500/151‑1 (MQ=255)
                                                                    ccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATTACAGCCAGTTCGCCCCGCGCCTTTTTCGCGGGGGCCGGATT‑AAAGGCGGGGGGGCGGCCCTTTGCGAAGGGGGGACcggcggaaaa                                                                                >  2:98920/1‑144 (MQ=255)
                                                                               ggCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggcggc                                                                    <  1:265844/150‑1 (MQ=255)
                                                                                    tgtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAGCCAGTACGCCCGGCGCTTTTTTCGCGGCGTCCGTATC‑AAGGGGGGGGGGGCGTCCTTTGGCAATGGCGGGACCGACGGTATAACCATAGGCGGCGGTggg                                                              >  1:829094/1‑151 (MQ=255)
                                                                                      tttttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACACCCAGTACCCCCCGCGCTTTTTTCGCGGCGTCCGTATC‑AAGGGGGGGGGGGGGTCCTTTTGCAATGGCGGGAACGACGATATAGCCCTAGGGGGGGTTggggg                                                            >  1:937345/1‑151 (MQ=255)
                                                                                           gTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGctg                                                         <  2:811314/149‑2 (MQ=255)
                                                                                             cGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGcttctt                                                      <  1:666317/150‑1 (MQ=255)
                                                                                             cGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGcttctt                                                      <  1:190948/150‑1 (MQ=255)
                                                                                                    gCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCa                                               <  2:359427/150‑1 (MQ=255)
                                                                                                       gAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACCCCCCGCGCTTTTTTCGCGGCGTCCGTATC‑AAGGGCGGGGAGGCGTCCCTTGGCAATGGCGGGACCGACGATATAGCCCTAGGCGGCGGTGGGGGGTTCTTCATGCCCtt                                             >  1:53233/1‑149 (MQ=255)
                                                                                                         gTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGt                                          <  2:472252/150‑1 (MQ=255)
                                                                                                         gTGCCCCGGCGTTACTGGCGGTAAAGACAGCCAACAAGCCCCGCGCTTTTTTCGCGGCGGCCGGATC‑AAGGGGGGGGGGGGGGCCTTTTGTAAAAGCGGAAACCGCCAAATAGCACTAAGCGGCGGTGGGGGGGTCTTCATGcccattgt                                          >  2:832033/1‑144 (MQ=255)
                                                                                                                 gCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGc                                 <  1:491176/151‑1 (MQ=255)
                                                                                                                               aaTGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTgccgggccgg                     >  2:57419/1‑149 (MQ=255)
                                                                                                                               aaTGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCgg                     >  1:966780/1‑149 (MQ=255)
                                                                                                                                 tGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAg                  >  2:516035/1‑150 (MQ=255)
                                                                                                                                    cAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGtttt              <  2:260481/151‑1 (MQ=255)
                                                                                                                                    cAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGttt               <  2:256662/150‑1 (MQ=255)
                                                                                                                                      gCCAGTACGCCGGCCGCCTTTTGCGCGGCGTCCGTCTC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGCCGATATAGCCACAGCCCGCGGTCGGGGCTGCTTCATTCCATTCGTACGCGTAGCGTGCCGTGCCGGTCGGTTTc             >  2:425516/1‑150 (MQ=255)
                                                                                                                                      gCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCa            <  1:588488/151‑1 (MQ=255)
                                                                                                                                        cAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGAGGATATAGCCATAGGCTGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGc          <  2:824687/151‑1 (MQ=255)
                                                                                                                                               cccGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATGGGCGGCGTTGGGTGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCTGt   <  2:176986/151‑1 (MQ=255)
                                                                                                                                                ccGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTc  <  2:515789/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTC  >  NC_000913/300565‑300859

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