Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913298,9600AC42.6% 80.8 / 26.0 55Y659D (TAT→GAT) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (28/3);  new base C (0/23);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 2.39e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC  >  NC_000913/298813‑299107
                                                                                                                                                     |                                                                                                                                                   
tataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAcc                                                                                                                                                   <  1:412212/150‑1 (MQ=255)
    cacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGcc                                                                                                                                              >  1:69924/1‑151 (MQ=255)
      ccGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCat                                                                                                                                            >  1:748881/1‑151 (MQ=255)
                            cACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                      >  1:202172/1‑151 (MQ=255)
                             aCGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                      >  2:1033440/1‑150 (MQ=255)
                             aCGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                      >  1:587551/1‑150 (MQ=255)
                              cgccgcACGGCCCCAGCCCCCCGACCCCTTTGGCCTTCATGGGGGAGGTTAGGGGGTTGGTATCACCCAGGAAAACCCCCTCCTGTTTTGGGATATCCGCATGACCCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCaa                                                                                                                    <  1:621770/151‑1 (MQ=255)
                                ccgcACAGGCCCAGCCCCCCGACCCCTTTGCCCTTCATGGGGGAGGTTAGGGGGTCGGTATCACCCAGAAAAATCCCCTCCTGTTTGGGAATACCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGa                                                                                                                              <  2:926561/139‑1 (MQ=255)
                                ccccaccgcccccGCTCCCCCACCCCTTTGGCCTTCATCGGGGAGGATAGGGTGTCGGTATCATCCAGAAAAATCCCCTCCTGTTTTGGGATACCCGCATCAACCGACACCTCATTCCCCGCCATATCGTGATTAACGAAGTAGCCCaaa                                                                                                                   <  2:124338/141‑1 (MQ=255)
                                 cgcACAGGCCCAGCTCACCGACCCCTTTGGCCTTCAGCGGGGAGGATAGGGGGTGGGTATCACCCAGAAAAACCCCCTCCTGTTTTGGGATACCCGCATGAACCGGCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACg                                                                                                                 <  1:371704/151‑1 (MQ=255)
                                   cacaGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGt                                                                                                               >  1:944291/1‑151 (MQ=255)
                                       ggCCCAGCTCACCGACCCCTTTGGCCTTCACGGGGGAGGATATGGGGTGGGTATCACCCAGGAAAACCCCCTCTTTTTTTGGGATATCCGCAGGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt                                                                                                            <  1:70652/150‑1 (MQ=255)
                                               tCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGcc                                                                                                    >  2:362125/1‑150 (MQ=255)
                                                ccccGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGTTATCACCCAGAAAAATCCCCTCCTGTTTTGGGATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa                                                                                                   <  1:369093/148‑1 (MQ=255)
                                                cACCGCCCCCTTTGGCTTTCATCGGGGAGGATAGGGGGTGGGTATCACCCAGGAAAATCCCCTCCTTTTTTGGGATACCCGCATGACCCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCg                                                                                                      <  1:75221/147‑1 (MQ=255)
                                                  ccGACCCCTTTGGCCTTCATGGGGGAGGATATGGGGTCGGTATCACCCAGAAAAACCCCCCCCTGTTTTGGAATCACCGCAGAACCCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGct                                                                                                <  2:652917/151‑1 (MQ=255)
                                                   cGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctc                                                                                               >  2:918824/1‑151 (MQ=255)
                                                      cacCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcctc                                                                                            >  2:563079/1‑151 (MQ=255)
                                                          tttGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc                                                                                        >  2:631972/1‑151 (MQ=255)
                                                           ttGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc                                                                                        >  1:58036/1‑150 (MQ=255)
                                                               cctacacggggaAGGTTAGGGGGTCGTTTCCCCCAAGAAAACCCCCCCCTTTTTTTGGAATACCCGAATAACCCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc                                                                                   <  2:150369/143‑1 (MQ=255)
                                                                cTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGcc                                                                                  >  2:1058432/1‑151 (MQ=255)
                                                                 ttcaggggggggGATATGGGGTCGTTACCACCCGGAAAAACCCCCTCCTGTTTGGGGATACCCCCAGGACCCGCCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc                                                                                   <  1:354197/143‑1 (MQ=255)
                                                                     tCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTGACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcc                                                                             >  2:484921/1‑151 (MQ=255)
                                                                     tCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcc                                                                             >  2:173404/1‑151 (MQ=255)
                                                                       ggggAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccc                                                                            >  2:5595/1‑150 (MQ=255)
                                                                        gggAGGATGGGGGGTCGTTATCACCAAGAAAACCCCCCTCTTTTTTGGGGATCCCCACATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                          <  2:133382/151‑1 (MQ=255)
                                                                          gaggaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgc                                                                        >  1:511531/1‑151 (MQ=255)
                                                                             gaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccc                                                                      >  2:379773/1‑150 (MQ=255)
                                                                                  gggtgggttaccacccAGAAAAACCCCCGCCTTTTTGGGAATACCGGAAGAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                                <  2:216147/137‑1 (MQ=255)
                                                                                          tATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc                                                        >  1:662357/1‑151 (MQ=255)
                                                                                          gtaccccctagAAAATCCCCCCGTTTTTTGAGATATCTGAAGGAGCGGACACCCTCTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAACCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcg                                                          <  1:505619/139‑1 (MQ=255)
                                                                                            tcatcCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCa                                                       >  2:534512/1‑150 (MQ=255)
                                                                                            tcatcCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCa                                                       >  1:659281/1‑150 (MQ=255)
                                                                                                 cAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAcc                                                 >  2:463203/1‑151 (MQ=255)
                                                                                                   gaaaaaTCCCCTCCTTTTTGGGGATATCCGAATGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGa                                                   <  1:401903/145‑1 (MQ=255)
                                                                                                     aaaaTCCCCTCTTTTTTGGGAATACCCGAAGAAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt                                                <  1:429037/148‑1 (MQ=255)
                                                                                                       aaTCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg                                           >  2:1062577/1‑151 (MQ=255)
                                                                                                          cccctcttTTTTTGGGATCCCCGAAGGACCCGCCCCCCCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg                                         <  2:818136/148‑1 (MQ=255)
                                                                                                                 tGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt                                 >  1:877243/1‑151 (MQ=255)
                                                                                                                  gTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt                                >  1:743365/1‑151 (MQ=255)
                                                                                                                    ttttGGGATACCCGAATAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtt                                  <  1:1058751/148‑1 (MQ=255)
                                                                                                                     ttggggATACCCGAATGACCCGCCCCCACACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTc                               <  2:511380/147‑1 (MQ=255)
                                                                                                                      tggggATCCCCGAATAAACCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCg                              <  1:233364/148‑1 (MQ=255)
                                                                                                                             acccGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCa                        <  1:1033761/147‑1 (MQ=255)
                                                                                                                              tCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg                      >  1:286153/1‑150 (MQ=255)
                                                                                                                              cccGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc                         <  2:105024/146‑1 (MQ=255)
                                                                                                                               ccGCATAACCCGCCACCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAg                       <  1:191389/148‑1 (MQ=255)
                                                                                                                                 gCAGGACCCGACCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg                   <  1:1170893/150‑1 (MQ=255)
                                                                                                                                    tGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg                   <  2:351747/147‑1 (MQ=255)
                                                                                                                                         ccgcccccTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt             <  2:17030/143‑1 (MQ=255)
                                                                                                                                          cGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc           >  1:290913/1‑150 (MQ=255)
                                                                                                                                              aCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac      >  1:242859/1‑151 (MQ=255)
                                                                                                                                               ccTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcaca     <  1:432010/151‑1 (MQ=255)
                                                                                                                                                   aTACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacacac  >  1:879650/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC  >  NC_000913/298813‑299107

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