Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,6610TC50.0% 39.9 / 18.8 48I44V (ATC→GTC) dppBdipeptide/heme ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (18/6);  new base C (0/24);  total (18/30)
Fisher's exact test for biased strand distribution p-value = 3.68e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAA‑TGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTAT  >  NC_000913/3705541‑3705806
                                                                                                                         |                                                                                                                                                 
ttttCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAA‑TGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCgg                                                                                                                     >  1:520010/1‑151 (MQ=255)
ttttCATTGAAATGCCTAGATCGCCATGCATAACGCCCAAAA‑TGTAAGGGAGATACTGCTGCCACATCGTTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCg                                                                                                                      <  2:675501/150‑1 (MQ=255)
tttgCATTAAAAGCCCTAGACCGCCTGGCA‑AACCCCCCAAAGTTTAGTGGAAATCCGGCGCCCACATTGTTTTATCTAACCCGAGTTCCGCGAGGAGCTGCGCGGGACCCTGTGGGGAGTTCCCCCTTCCGCCCGCCATGATCATCAcc                                                                                                                       <  1:890162/149‑1 (MQ=255)
      ttaaaagcccTGGATCCCCATGCATACCCCCCAAAA‑TGTAATGGAGATCCGGCTCCAACATCGTTTTATCAAAGCCGAGTCCAGCCGGCGGCTGCGCGTGCCGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATc                                                                                                                  <  2:54082/140‑1 (MQ=255)
          aaTGCCTAGCTCGCAATGCATACCCCCCAAAA‑TGTAAGGGAGATTCTGCTCCCACATCGTTTTATCTAAGCCGATTTCAGCCAGCAGCTGCGCGTGCCGCTCGGGGGAGACCCCACTTTCGCCCGCCATGATCATCACCGGATCGCCCggg                                                                                                           <  2:760468/151‑1 (MQ=255)
            tGCCTAGATCGCCATGCATAACGCCCCAAA‑TGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGAt                                                                                                         >  2:943077/1‑151 (MQ=255)
              ccTAGATCGCCATGCATAACGCCCCAAA‑TGTAATGGAGATACTGCTGCCACTTCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATc                                                                                                        >  2:528328/1‑150 (MQ=255)
              ccTAGACCGCCATGCATAACCCCCAAAA‑TGTAATGGAGATCCGGCTCCCACATCGGTTTATCAAAGCCGAGTTCAGCCAGCAGCGGCGGGGGACGCTCGGGGGAGCCCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATc                                                                                                        <  2:916601/150‑1 (MQ=255)
               cTAGATCGCCATGCATAACGCCCCAAA‑TGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATTCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAt                                                                                                      >  2:777601/1‑151 (MQ=255)
                 aGATCGCCATCCATACCCCCCAAAA‑TGTAAGGGAGATCCTGCTCCCACTTCGTTTTATCTAAGCCGAGTTCACCCAGCGGCTGCGCGGGCCGCTCTGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtgt                                                                                                    <  2:348831/151‑1 (MQ=255)
                         aTGCATAACGCCCCAAA‑TGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAACTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACaa                                                                                              >  2:212251/1‑149 (MQ=255)
                            catacccccaaaaa‑TGAAAGGGATATACGGCTGCCACACGGTTTTACCAACGCCGAGTGCGGCCAGCAGCGGCGCGTGCCGTTCGGGGAAGATCCCCCTTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCAtgtg                                                                                                   <  2:588223/130‑1 (MQ=255)
                            cATAACCCCCCAAT‑TTTAAGGGAGATACGGCTCCCCCATCGTTTTATCAAGGCCGAGTCCACCCAGCAGCTGCGCGGGACGCTGGGGGGAAACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGc                                                                                          <  2:80630/150‑1 (MQ=255)
                               accccccaAAA‑TGTAAGGGAGATACTGCTCCCACATCGTTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCGGGGGAGACCCCCCTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCa                                                                                         <  2:693185/144‑1 (MQ=255)
                               accGCCCAAAA‑TGTAAGGGAGATACTGCTCCCACATCGTTTTATCTAAGCCGATTTCAGCCAGCAGCTGCGCGTGACGCTCGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAgg                                                                                           <  2:942552/144‑1 (MQ=255)
                                accccccAAATTGAAAGGGAGATCCTGCTGCCACAGGGTTTTATCAAAGCCGATTTCAGCCGGCAGCGGCGCGGGCCGCTGGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCa                                                                                         <  1:402904/145‑1 (MQ=255)
                                  ccccaAAA‑TGTAATGGAAATCCTGCTGCCACTTCGTTTTATCAAAGCCGATTTCAGCCAGCAGCGGCGGGTGCCGCTGGGGGGAGCCCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtgtg                                                                                   <  2:603517/150‑1 (MQ=255)
                                       aaa‑TGTAATGGAAATCCTGCGCCCACATCGTTTTATCAAAGCCGAGTTCAGCCAGCAGCTGGGCGGGCCGCTCGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAg                                                                              <  1:156201/151‑1 (MQ=255)
                                                tGGAGATACTGCTGCCACATCGGTTTATCTAATCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGACCGCCCGGTATCATGTGCACAAAGGCAAATGTGAGAAGGGTAATAc                                                                      >  1:617643/1‑151 (MQ=255)
                                                   agaTACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCg                                                                    >  1:664605/1‑150 (MQ=255)
                                                           cgcccACATGGGTTTATCTAAGCCGATTTCAGCCAGCAGCTGCGCGGGCCGCTCGGGGGAGATCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGAt                                                                  <  1:437026/141‑1 (MQ=255)
                                                              ccACATCGTTTTATCTAAGCCGAGTTCAGCCAGCGGCGGCGGGTGACGCTCGGGGGAGACCCCACGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTgggg                                                        <  2:85291/151‑1 (MQ=255)
                                                               cacaTGGTTTTATCTAAGCCGAGTTCAGCCGGCAGCTGCGCGTGCCGCTCTGGGGAGACCCCCCTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGa                                                       <  2:797839/151‑1 (MQ=255)
                                                                acaGCGGTTTATCAAAGCCGATTTCAGCCAGCAGCTGCGCGTGACGCTGGGGGGAGCCCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGa                                                       <  2:477448/150‑1 (MQ=255)
                                                                       tttATCTAAGCCGAGTTCAGCCACCAGCGGCGGGTGCCCCTGGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGAt                                                      <  1:688181/144‑1 (MQ=255)
                                                                          aTCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTccc                                            >  2:465610/1‑151 (MQ=255)
                                                                            cTAAGCCGAGTTCAGCCAGGCGCTGCGCGTGGCGCTGGGGGGAGCTCCCTCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTc                                              <  1:912771/147‑1 (MQ=255)
                                                                              aggCCGAGTCCAGCCACCCGGCGCGCGGCACGCTCGGGGGGAACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAAc                                        <  2:547715/149‑1 (MQ=255)
                                                                              aaGCCGAGTACAGCCAGCGGCTGCGCGTGCCGCTGTGGGGGGCTCCCGGGTGGGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCa                                           <  1:955796/148‑1 (MQ=255)
                                                                                 ccGAGTCCAGCCAGCAGCGGCGCGTGCCGCTGGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCa                                           <  1:633131/145‑1 (MQ=255)
                                                                                  cGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCg                                    >  1:658272/1‑151 (MQ=255)
                                                                                     gTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgg                                   >  2:534204/1‑149 (MQ=255)
                                                                                      tcccgccagcGGCGGCGCGTGCCGCTCGGGGGAGCCCCCCGGGTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGga                                  <  1:516800/145‑1 (MQ=255)
                                                                                       ccagccaggAGCGGCGCGTGCCCCTGGGGGGAGACCCCACGGCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGgaga                                <  1:446392/149‑1 (MQ=255)
                                                                                            cagcagCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAAc                          >  1:131448/1‑151 (MQ=255)
                                                                                             agcagcTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACt                         >  2:835688/1‑151 (MQ=255)
                                                                                                 gcggcgcGTGGCGCTGGGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTc                                     <  1:549997/132‑1 (MQ=255)
                                                                                                  cTGCGCGCGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGCCGAGTCCCAAACGTCGGAGAATAAACTGCaa                     >  1:375348/1‑150 (MQ=255)
                                                                                                    gggcgTGACGCTGGGGGGACACCCCACTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTg                        <  1:941373/143‑1 (MQ=255)
                                                                                                      gcgTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATa                 >  1:558977/1‑150 (MQ=255)
                                                                                                        gTGACGCTGGGGGGAGACCCCCCGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCaa                     <  2:6700/144‑1 (MQ=255)
                                                                                                        gTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAAcc              >  1:109049/1‑151 (MQ=255)
                                                                                                          gACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCg            >  2:479004/1‑151 (MQ=255)
                                                                                                             gCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGAt         >  1:308725/1‑151 (MQ=255)
                                                                                                             gCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGAt         >  1:281077/1‑151 (MQ=255)
                                                                                                                  ggggAGCCCCCCCGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAAccc             <  1:424451/142‑1 (MQ=255)
                                                                                                                  ggggAGACCCCCCGTTCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGAtt        <  2:352608/147‑1 (MQ=255)
                                                                                                                    ggAGACCCCCCGTTCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCtat  <  2:904657/151‑1 (MQ=255)
                                                                                                                         |                                                                                                                                                 
TTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAA‑TGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTAT  >  NC_000913/3705541‑3705806

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