Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,6530TG48.9% 46.1 / 17.2 47P46P (CCA→CCC) dppBdipeptide/heme ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (19/5);  new base G (0/23);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 7.37e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAA‑TGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGAT  >  NC_000913/3705504‑3705799
                                                                                                                                                      |                                                                                                                                                  
aaGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAA‑TGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTg                                                                                                                                                   >  2:53578/1‑151 (MQ=255)
   cgcgGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAA‑TGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTgggg                                                                                                                                                <  1:399478/151‑1 (MQ=255)
                                     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)
                                                                                                                                                      |                                                                                                                                                  
AAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAA‑TGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGAT  >  NC_000913/3705504‑3705799

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: ‑