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

GCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCT  >  NC_000913/3705511‑3705804
                                                                                                                                                      |                                                                                                                                               
gCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGAt                                                                                                                                                 <  2:256993/151‑1 (MQ=255)
                        tGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGatcat                                                                                                                         >  2:626666/1‑151 (MQ=255)
                              ttttCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCgg                                                                                                                   >  1:82065/1‑151 (MQ=255)
                              ttttCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCgg                                                                                                                   >  1:516068/1‑151 (MQ=255)
                                  cATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCg                                                                                                               >  1:442741/1‑151 (MQ=255)
                                   aTTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGc                                                                                                              >  2:606403/1‑151 (MQ=255)
                                    taaaaaTCCCTGGATCCCCATGCAAACCCCCCAAAATGTAAGGGAGATCCTGCTCCCACATCGTTTTATCTAAGCCGAGTCCAGCCGGCAGCTGCGCGTGCCGCTCGGGGGAGACCCCCCTTTCGCCCGCCATGATCATCACCGGATCGc                                                                                                              <  1:378187/147‑1 (MQ=255)
                                           gCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATc                                                                                                      >  2:293830/1‑151 (MQ=255)
                                               aGATCGCCATGAATACCGCCAAAAATGTAAGGGAGATACTGCTGCCACATGGTTTTATCAAAGCCGAGTTCAGCCAGCAGCTGCGCGGGCCGCTCGGGGGAGACCCCCCTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtgt                                                                                                  <  1:549955/151‑1 (MQ=255)
                                               aGACCGCCATGCAAAACCCCCAAAATGTAAGGGAGATACTGCTCCCACATCGGTTTATCAAAGCCGAGTCCCGCCAGCAGCTGCGCGTGACGCTCGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAt                                                                                                    <  1:487253/149‑1 (MQ=255)
                                                      cATGCATAACGCCCAAAATGTAAGGGAGATACTGCTGCCACATGGGTTTATCAAACCCGAGTTCAGCCAGCAGCTGGGCGTGACGCTGGGGGGAGATCCCGCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACa                                                                                             <  1:44690/149‑1 (MQ=255)
                                                           aaacccccccAAATGAAAGGGAAATCCGGCCCCCACTGCGTTTTATCAAGGCGGAGTCCAGCCAGCGGCGGCGGGTGCCGCTGGGGGGAAACCCCACGTTCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAGGGCa                                                                                       <  1:116242/144‑1 (MQ=255)
                                                             accccccaAAATGTAAGGGAGATCCTGCTCCCACATCGTTTTATCAAAGCCGAGTTCAGCCAGCAGCTGCGCGTGCCGCTGGGGGGAGACCCCCCTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCaa                                                                                      <  2:400393/145‑1 (MQ=255)
                                                                  ccaaaatgaaatgGAGATACTGCTGCCACATCGTTTTATCAAAGCCGAGTTCAGCCAGCAGCGGCGCGTGACGCTCGGGGGAGACCCCACTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTgag                                                                               <  2:352508/148‑1 (MQ=255)
                                                                   ccAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAtgtg                                                                                 >  1:278945/1‑148 (MQ=255)
                                                                    cAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTgaga                                                                              >  1:296759/1‑150 (MQ=255)
                                                                     aagtttAAGGGAGATCCTGCTCCCACAGCGGTTTATCAAAGCCGAGTTCAGCCAGCAGCGGCGCGTGCCGCTGTGGGGAGACCCCCCGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTgag                                                                               <  2:508763/143‑1 (MQ=255)
                                                                     aaaTTTAAGGGAGATCCTGCTCCCACATCGTTTTATCAAGGCCGTGTTCAGCCGGCAGCGGCGCGTGACGCTCGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTgaga                                                                              <  1:411255/149‑1 (MQ=255)
                                                                     aaaTGTAATGGAGATCCTGCTCCCACATCGTTTTATCAAGGCGGATTTCAGCCGGCAGCGGCGCGTGCCGCTCGGGGGAGATCCCACGTTCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAg                                                                            <  2:290493/151‑1 (MQ=255)
                                                                     aaaTGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTgaga                                                                              >  1:626551/1‑149 (MQ=255)
                                                                        tGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGt                                                                         >  1:261674/1‑151 (MQ=255)
                                                                           aaTGGAGATCCTGCTCCCACCTCGTTTTATCTAAGCCGAGTCCAGCCAGCAGCTGCGCGTGCCGCTCGGGGGAGACCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAAt                                                                      <  1:79170/151‑1 (MQ=255)
                                                                           aaTGGAGATACTGCTGCCCCATCGTTTTATCTAAGCCGAGTTCACCCAGCAGCTGCGCGTGACGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTa                                                                        <  1:451352/149‑1 (MQ=255)
                                                                              ggAGATCCTGCTCCCACATGGTTTTATCTAAGCCGATTTCAGCCAGCGGCTGCGCGGGACGTTGTGGGGAGACCCCACTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTaa                                                                       <  2:264345/147‑1 (MQ=255)
                                                                              ggAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAcc                                                                   >  1:240120/1‑151 (MQ=255)
                                                                                    cggttcgccctttCGGTTTTTCTACGCCAGGTGCGCCGACCAGTTGCGGGAGGCGCTGGGGGGAGTCCCTGCTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATaa                                                              <  1:283865/138‑1 (MQ=255)
                                                                                           ccACATCGGTTTATCTAAGCCGATTCCAGCCAGCGGCGGCGCGGGCCGCTCGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTgggg                                                      <  1:645702/151‑1 (MQ=255)
                                                                                           ccACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTgggg                                                      >  1:310044/1‑151 (MQ=255)
                                                                                             acgtctgTTTATCAAACCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTGGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTgg                                                        <  1:643936/141‑1 (MQ=255)
                                                                                              catggttttATCTAAGCCGAGTTCAGCCAGCAGCGGCGCGGGCCGCTCTGGGGAGATCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATg                                                   <  1:421557/145‑1 (MQ=255)
                                                                                                          tAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCaa                                        >  1:327272/1‑150 (MQ=255)
                                                                                                            aGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACg                                     >  2:96855/1‑151 (MQ=255)
                                                                                                              ccGATTTCAGCCAGCAGCGGCGCGTGACGCTCGGGGGAGCCCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTcc                                           <  2:4706/143‑1 (MQ=255)
                                                                                                                     cagccagcAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAAt                            >  1:270234/1‑151 (MQ=255)
                                                                                                                         cagcagCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAAc                        >  1:261998/1‑151 (MQ=255)
                                                                                                                            cagcGGCCCTTGACGCTCGGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACt                       <  2:165992/149‑1 (MQ=255)
                                                                                                                             agcTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGc                     >  1:415390/1‑150 (MQ=255)
                                                                                                                                    cgTGAGGCTCTGGGGAGATCCCAAGTTCGCACGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATaa              >  1:59665/1‑150 (MQ=255)
                                                                                                                                               ggggAGCCCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCg          <  1:221060/143‑1 (MQ=255)
                                                                                                                                               ggggAGACCCCCCTTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATtct    <  2:497856/149‑1 (MQ=255)
                                                                                                                                                gggAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTctct  >  1:238632/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
GCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCT  >  NC_000913/3705511‑3705804

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