Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913300,7210GT42.1% 61.3 / 16.7 38Q72K (CAA→AAA) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base G (5/17);  new base T (16/0);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.04e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGT  >  NC_000913/300586‑300858
                                                                                                                                       |                                                                                                                                           
ccACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCCGTACGCCCGGCGCTTTTTTCgcgg                                                                                                                                        >  2:179504/1‑141 (MQ=255)
          cgcAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGAAAGCCCGTACGCCCCGCGCCTTTTTCGCGGCGTcccggtcaaaggggggg                                                                                                                    >  2:465434/1‑136 (MQ=255)
           gcAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTg                                                                                                                   <  1:419625/150‑1 (MQ=255)
             aaTGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGAAAGCCAGTACCCCCCGCGCCTTTTTTGCGGCCGCCGGGTC‑AAggggggggggg                                                                                                                >  2:726010/1‑146 (MQ=255)
              aTGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAACCCCTAATCCCCGCGCTTTTTTTGCGGCGGCCGTATC‑AAGGGGGGTggggg                                                                                                               >  1:152072/1‑150 (MQ=255)
                  ccTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCCGTACACCCGGCGCTTTTTTTGGGGGGGCCGGATT‑AAgggggggggggggggc                                                                                                           >  2:772253/1‑141 (MQ=255)
                     gatgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCttt                                                                                                        <  2:172470/151‑1 (MQ=255)
                     gatgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATTAAAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGTCCGTATC‑CAAGGCggggggggggccctt                                                                                                        >  2:344806/1‑138 (MQ=255)
                       agatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAAGAAAAACAGTACCCCCGGGGCCTTTTTTGCGGCGGCCCGTTC‑AAAGGGGGGGAgggggccctttg                                                                                                      >  2:598567/2‑141 (MQ=255)
                               cTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATTAAAGCCCATACGCCCCGCGCTTTTTTCGCGGGGGCCGTATC‑AAGGGGGGGGGGGCGGCCTTTTGCAATgg                                                                                                >  2:753099/1‑149 (MQ=255)
                                 cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTTATGAAAACCACTACGCCCCGCGCTTTTTTCGCGGCGCCCCCATC‑AAGGGGGGGGGGGCGTCCttttggaacgg                                                                                                >  2:150075/1‑139 (MQ=255)
                                 cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCg                                                                                              <  2:812219/149‑1 (MQ=255)
                                 cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCg                                                                                              <  1:307429/149‑1 (MQ=255)
                                 cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGACAACCCGGAAGCCCCGCGCTTTTTTCGCGGCGGCCGGATCAAAGGGGGGGGGGGGGCCCTTTGGCatggggg                                                                                              >  2:999336/1‑144 (MQ=255)
                                   aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAAc                                                                                          <  1:839301/151‑1 (MQ=255)
                                   aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGGACGCCCGGCGCTTTTTTTGCGGCGGCCGGATC‑AAGGGCGGGGGGGCGGCCTTTGGCGAAGGCggga                                                                                           >  2:146899/1‑148 (MQ=255)
                                             cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGAtata                                                                                <  2:476200/151‑1 (MQ=255)
                                                       ccTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑ggcg                                                                     <  2:106857/151‑1 (MQ=255)
                                                          ggCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑Ggcggcg                                                                  <  1:140086/151‑1 (MQ=255)
                                                                    ttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGct                                                         <  2:318471/150‑1 (MQ=255)
                                                                           cTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGGGGCCGGATC‑AAAGGCGGGGAGGCGGCCTTTGGCAATGGCGGGAACGCCGATATAACCATAGGGGGGGGTTGGGGGctctt                                                      >  2:713620/1‑144 (MQ=255)
                                                                            tttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAGCCAGTAAGCCCGGCGCCTTTTTCGCGGGGGCCGGATC‑AAAGGGGGGGGGGGGGCCTTTGGGAATGGGGGGAACGGCGATATAGGCCTA‑GGGGGGGGTgggggggtttttat                                                  >  2:296711/1‑140 (MQ=255)
                                                                             ttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAAGCCCGGCGCCTTTTGCGCGGCGTCCGTATC‑AAGGGCGGGGGGGCGGCCTTTTGCAATGGCGGGACCGACGATATAGCCATA‑GGCGGCGTTGGGGGGTTTTTCATGcc                                               >  2:773890/1‑151 (MQ=255)
                                                                              tGCCGAGTGCCCCGGCGTTACTGGCGGTAAATACAACCCCCACGCCCGGCGCCTTTTTCGCGGGGTCCGGTTC‑CAGGGGGGTGGGGGGGCCTCTTGCAAAGGGGGGACCGACGCAATAACCCTA‑AGCGGGGGTgggggggtttt                                                     >  1:1010131/1‑138 (MQ=255)
                                                                                ccGAGTGCCCCGGCGTTACTGGCGGTAAAGACAACCCGGACCCCCCGCGCTTTTTTCGCGGCGGCCGGATC‑AAGGGGGGGGGGGGGTCCTTTGGCAAAGGCGGGACCGCCGGAATAGCCCTA‑GGCGGCGTTGGGGGGGTCTTTATGCCCtt                                            >  1:1011571/1‑151 (MQ=255)
                                                                                 cGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAAGCCCCGCGCTTTTTTCGCGGCGGCCGGATC‑AAGGGGGGGGAGGCGGCCTTTTGGCATGGCGGGAACGACGGGATAACCATA‑GGGGGGGGTGGGGGGGTCTTTATTCCCtt                                            >  1:666868/1‑150 (MQ=255)
                                                                                      gCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAgg                                      <  1:676238/151‑1 (MQ=255)
                                                                                       ccccGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAgg                                      <  1:365051/150‑1 (MQ=255)
                                                                                       ccccGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAgg                                      <  2:835058/150‑1 (MQ=255)
                                                                                         ccGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGt                                   <  2:43351/151‑1 (MQ=255)
                                                                                             cGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcg                               <  2:357044/151‑1 (MQ=255)
                                                                                                    ggCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtgccg                          <  1:976213/149‑1 (MQ=255)
                                                                                                      cGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCgtgccgtgcc                      >  2:1015088/1‑151 (MQ=255)
                                                                                                        gTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGGAGCGTGCCGTGCCgg                    >  2:735865/1‑151 (MQ=255)
                                                                                                         tAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCgg                    >  1:444341/1‑150 (MQ=255)
                                                                                                           aTGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAg                 <  1:636879/151‑1 (MQ=255)
                                                                                                                         gCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCg   >  1:12011/1‑151 (MQ=255)
                                                                                                                          cccGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGt  <  2:783652/151‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                           
CCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGT  >  NC_000913/300586‑300858

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