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

CCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGT  >  NC_000913/300586‑300858
                                                                                                                                       |                                                                                                                                           
ccACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCCGTACGCCCGGCGCTTTTTTCgcgg                                                                                                                                        >  2:179504/1‑141 (MQ=255)
          cgcAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGAAAGCCCGTACGCCCCGCGCCTTTTTCGCGGCGTcccggtcaaaggggggg                                                                                                                    >  2:465434/1‑136 (MQ=255)
           gcAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTg                                                                                                                   <  1:419625/150‑1 (MQ=255)
             aaTGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGAAAGCCAGTACCCCCCGCGCCTTTTTTGCGGCCGCCGGGTCAAggggggggggg                                                                                                                >  2:726010/1‑146 (MQ=255)
              aTGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAACCCCTAATCCCCGCGCTTTTTTTGCGGCGGCCGTATCAAGGGGGGTggggg                                                                                                               >  1:152072/1‑150 (MQ=255)
                  ccTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCCGTACACCCGGCGCTTTTTTTGGGGGGGCCGGATTAAgggggggggggggggc                                                                                                           >  2:772253/1‑141 (MQ=255)
                     gatgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCttt                                                                                                        <  2:172470/151‑1 (MQ=255)
                     gatgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATTAAAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGTCCGTATCCAAGGCggggggggggccctt                                                                                                        >  2:344806/1‑138 (MQ=255)
                       agatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAAGAAAAACAGTACCCCCGGGGCCTTTTTTGCGGCGGCCCGTTCAAAGGGGGGGAgggggccctttg                                                                                                      >  2:598567/2‑141 (MQ=255)
                               cTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATTAAAGCCCATACGCCCCGCGCTTTTTTCGCGGGGGCCGTATCAAGGGGGGGGGGGCGGCCTTTTGCAATgg                                                                                                >  2:753099/1‑149 (MQ=255)
                                 cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTTATGAAAACCACTACGCCCCGCGCTTTTTTCGCGGCGCCCCCATCAAGGGGGGGGGGGCGTCCttttggaacgg                                                                                                >  2:150075/1‑139 (MQ=255)
                                 cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCg                                                                                              <  2:812219/149‑1 (MQ=255)
                                 cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCg                                                                                              <  1:307429/149‑1 (MQ=255)
                                 cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGACAACCCGGAAGCCCCGCGCTTTTTTCGCGGCGGCCGGATCAAAGGGGGGGGGGGGGCCCTTTGGCatggggg                                                                                              >  2:999336/1‑144 (MQ=255)
                                   aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAAc                                                                                          <  1:839301/151‑1 (MQ=255)
                                   aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGGACGCCCGGCGCTTTTTTTGCGGCGGCCGGATCAAGGGCGGGGGGGCGGCCTTTGGCGAAGGCggga                                                                                           >  2:146899/1‑148 (MQ=255)
                                             cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGAtata                                                                                <  2:476200/151‑1 (MQ=255)
                                                       ccTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggcg                                                                     <  2:106857/151‑1 (MQ=255)
                                                          ggCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGgcggcg                                                                  <  1:140086/151‑1 (MQ=255)
                                                                    ttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGct                                                         <  2:318471/150‑1 (MQ=255)
                                                                           cTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGGGGCCGGATCAAAGGCGGGGAGGCGGCCTTTGGCAATGGCGGGAACGCCGATATAACCATAGGGGGGGGTTGGGGGctctt                                                      >  2:713620/1‑144 (MQ=255)
                                                                            tttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAGCCAGTAAGCCCGGCGCCTTTTTCGCGGGGGCCGGATCAAAGGGGGGGGGGGGGCCTTTGGGAATGGGGGGAACGGCGATATAGGCCTAGGGGGGGGTgggggggtttttat                                                  >  2:296711/1‑140 (MQ=255)
                                                                             ttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAAGCCCGGCGCCTTTTGCGCGGCGTCCGTATCAAGGGCGGGGGGGCGGCCTTTTGCAATGGCGGGACCGACGATATAGCCATAGGCGGCGTTGGGGGGTTTTTCATGcc                                               >  2:773890/1‑151 (MQ=255)
                                                                              tGCCGAGTGCCCCGGCGTTACTGGCGGTAAATACAACCCCCACGCCCGGCGCCTTTTTCGCGGGGTCCGGTTCCAGGGGGGTGGGGGGGCCTCTTGCAAAGGGGGGACCGACGCAATAACCCTAAGCGGGGGTgggggggtttt                                                     >  1:1010131/1‑138 (MQ=255)
                                                                                ccGAGTGCCCCGGCGTTACTGGCGGTAAAGACAACCCGGACCCCCCGCGCTTTTTTCGCGGCGGCCGGATCAAGGGGGGGGGGGGGTCCTTTGGCAAAGGCGGGACCGCCGGAATAGCCCTAGGCGGCGTTGGGGGGGTCTTTATGCCCtt                                            >  1:1011571/1‑151 (MQ=255)
                                                                                 cGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTAAGCCCCGCGCTTTTTTCGCGGCGGCCGGATCAAGGGGGGGGAGGCGGCCTTTTGGCATGGCGGGAACGACGGGATAACCATAGGGGGGGGTGGGGGGGTCTTTATTCCCtt                                            >  1:666868/1‑150 (MQ=255)
                                                                                      gCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAgg                                      <  1:676238/151‑1 (MQ=255)
                                                                                       ccccGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAgg                                      <  1:365051/150‑1 (MQ=255)
                                                                                       ccccGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAgg                                      <  2:835058/150‑1 (MQ=255)
                                                                                         ccGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGt                                   <  2:43351/151‑1 (MQ=255)
                                                                                             cGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcg                               <  2:357044/151‑1 (MQ=255)
                                                                                                    ggCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtgccg                          <  1:976213/149‑1 (MQ=255)
                                                                                                      cGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCgtgccgtgcc                      >  2:1015088/1‑151 (MQ=255)
                                                                                                        gTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGGAGCGTGCCGTGCCgg                    >  2:735865/1‑151 (MQ=255)
                                                                                                         tAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCgg                    >  1:444341/1‑150 (MQ=255)
                                                                                                           aTGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAg                 <  1:636879/151‑1 (MQ=255)
                                                                                                                         gCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCg   >  1:12011/1‑151 (MQ=255)
                                                                                                                          cccGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGt  <  2:783652/151‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                           
CCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGT  >  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: