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

GAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATT  >  NC_000913/3705515‑3705800
                                                                                                                                                  |                                                                                                                                           
gAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATccc                                                                                                                                          >  1:296942/1‑150 (MQ=255)
          cAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGTCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTcgcccg                                                                                                                               >  1:190521/1‑151 (MQ=255)
          cAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTcgcccg                                                                                                                               >  1:371568/1‑151 (MQ=255)
              ccGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCAt                                                                                                                           >  2:356546/1‑151 (MQ=255)
                     gCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGatcatc                                                                                                                    >  2:230134/1‑151 (MQ=255)
                           tttaatgaaaaGGCCTAGATCCCCAGACACACCGCCCAAATTGTAAGGGAAATCCGGCTGCCACGTCGTTTTATCGAAGCCGAGTTCGGCCACCGCCTGCGCGTGCCGCGGGGGGGGAACCCCACGTGCCCCCGCCATGATCATCACCgg                                                                                                               <  2:161525/142‑1 (MQ=255)
                              cATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATc                                                                                                            >  1:77487/1‑150 (MQ=255)
                                ttGAAATGCCTAGACCGCCAGGCAAAACCCCCCAAATGTAAGGGAGATCCTGCTGCCACATCGGTTTATCTAGCCCGATTCCAGCCAGCAGCTGCGCGTGGCGCTCGGGGGAGACCCCGCTTTCGCCCGCCATGATCATCACCGGATCGc                                                                                                          <  2:50239/150‑1 (MQ=255)
                                   aaaTGCCTAGCCCCCCAGGAATACCCCCCCAAAGGAAAGGGAGATCCGGCGGCCACATGGTTTTAGCAAAGCCGAGTCCCCCCGGCGGCTGCGCGTGACGCGCGGGGGAGACCCCAGGTTCGCCCGCCATGATCATCACCGGATCGCCCgg                                                                                                      <  2:243394/151‑1 (MQ=255)
                                    aaTCCCTATATCCCCATGAATACCCCCCAAAATATAAGGGAGTTGCTGCCCCCCCTTCGTTTTATCTACGCGGTGTTCGCCCAGCACCTGCGCGTGACGTTCTGGGAAAACCCCCTTTCCGCCCGCCATGATCATCACCGGATCGCCCggg                                                                                                     <  1:28936/151‑1 (MQ=255)
                                            gtcccccATGCATACCCCCCAAAATGTAAGGGAGATCCTGCTGCCACATGGTTTTATCAAAGCCGATTTCCGCCAGCGGCTGCGCGGGCCGCTGTGGGGAGACCCCACTTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCAtg                                                                                               <  1:235124/144‑1 (MQ=255)
                                               ccccATGCATACCCCCCAAAATGTAATGGAGATCCTGCTCCCACATCGGTTTATCTAAGCCGAGTTCACCCAGCAGCGGCGGGTGACGCTCGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGAc                                                                                          <  1:160231/149‑1 (MQ=255)
                                                             ccaaaaaTGAAAGGGAGATCCTGCCCCCACATCTTTTTATCTAAGCCGAGTCCCGCCACCACCTGCGCGGGCCGCTGGGGGGAGACCCCCCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTga                                                                            <  2:296913/147‑1 (MQ=255)
                                                                  aatgaaaggggaatcctgctgCCCCTTCGTTTTATCAAGGCCGATTTCAGCCAGCAGCGGCGGGGGCCGCTCGGGGGAGACCCCCCGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAgg                                                                       <  1:96057/136‑1 (MQ=255)
                                                                  aaTGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAgg                                                                       >  1:618607/1‑151 (MQ=255)
                                                                     gTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTAGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGt                                                                     >  2:567415/1‑150 (MQ=255)
                                                                       aagggAGATGCTGCGCCCCCACCGGTTTATCTAAGCCGTTTTCCGCCAGCAGCTGCGCGTGACCCTCGGGGGAGACCCCTCGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTa                                                                    <  1:93266/146‑1 (MQ=255)
                                                                         tGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATAc                                                                >  2:320899/1‑151 (MQ=255)
                                                                          ggAGATACTGCTGCCACATCTTTTTATCAAAGCCGATTTCAGCAAGCAGCTGGGCGTGACGCTGGGGGGAGACCCCCCGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATa                                                                 <  1:611358/149‑1 (MQ=255)
                                                                            agaTCCTGCTGCCACATGGTTTTATCAAAGCCGACTGCAGCCGGCAGCTGCGGGTGGCGCTGGGGGGAGACCCCATGTGCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCg                                                              <  2:56962/150‑1 (MQ=255)
                                                                                     tgcCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTgg                                                    >  2:142855/1‑151 (MQ=255)
                                                                                       ccACACCGTTTTATCTACGCCGAGTTCAGCCAGCAGCGGGGCGTGCCGCTCGGGGGAGACCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTgg                                                    <  1:9541/149‑1 (MQ=255)
                                                                                             cGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACg                                            >  2:362304/1‑151 (MQ=255)
                                                                                               gTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACg                                            >  1:677203/1‑149 (MQ=255)
                                                                                                       aaGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAAc                                  >  2:569359/1‑151 (MQ=255)
                                                                                                          cgggttccccccggcagcTGCGGGGGCCGTTGGGGGGAGACCCCCCGTCCCCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTc                               <  1:219660/138‑1 (MQ=255)
                                                                                                                tCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGACGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGgaga                          >  2:474754/1‑150 (MQ=255)
                                                                                                                  ggccgcatgcggcgcGGGCCGGTCTGGGGGGATCCCACGTGGACCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTc                               <  2:58796/132‑1 (MQ=255)
                                                                                                                    cccgccgcTGGGCGGGACGCGCGGGGGAGACCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGgaga                          <  2:303357/140‑1 (MQ=255)
                                                                                                                         agcggggcgTGCCGCTCTGGGGAGCCCCCACGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACt                   <  2:508640/142‑1 (MQ=255)
                                                                                                                                  tGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAAccc       >  2:220714/1‑151 (MQ=255)
                                                                                                                                    cctctggggggACACCCCGGGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTTGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATagc         <  2:611319/140‑3 (MQ=255)
                                                                                                                                       ctctGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGAtt  >  1:600320/1‑151 (MQ=255)
                                                                                                                                          gggggAGGCCCCGCGTCCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGa    <  2:679291/145‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
GAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATT  >  NC_000913/3705515‑3705800

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