Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,560,7660TG42.9% 21.7 / 10.7 28Q55P (CAG→CCG) ddpBD,D‑dipeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (7/9);  new base G (0/12);  total (7/21)
Fisher's exact test for biased strand distribution p-value = 1.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

AAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAA  >  NC_000913/1560617‑1560912
                                                                                                                                                     |                                                                                                                                                  
aaGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTg                                                                                                                                                   <  1:784869/151‑1 (MQ=255)
      gCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGaa                                                                                                                                              >  2:651658/1‑150 (MQ=255)
          cAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATAttt                                                                                                                                         <  1:461081/151‑1 (MQ=255)
                cGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACg                                                                                                                                   <  1:719965/151‑1 (MQ=255)
                      gAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCa                                                                                                                             <  1:432439/151‑1 (MQ=255)
                          actataCGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGGTGCTGGGGGATATTTTTCACGATa                                                                                                                                >  1:791354/3‑144 (MQ=255)
                                cGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGgcgc                                                                                                                   <  2:264543/151‑1 (MQ=255)
                                         tCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCAccc                                                                                                           >  1:48038/1‑150 (MQ=255)
                                                                aaTGGATTTCCCCAGGCCCCAATGAAAAAGGTCGCTGCCGTAGGGTTAAACCTTTACATACAGTGGCTGTCCCGGCCCCAGTTGCGGGCAAATATTTCCCACGATACCCTCGCGGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCa                                                                                   <  2:671211/151‑1 (MQ=255)
                                                                   ggATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAggg                                                                                >  2:248475/1‑151 (MQ=255)
                                                                   ggATGTCCCCGGGCCCCCATAAAACGGGTCGCTGACGTGGCGGAAAAACTGTACATACAGGGGCGGGCCCAGCCCCGGTTGCGGGCGAATTTTTCCCACGATAGCACCGCGGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAg                                                                                  <  1:682119/149‑1 (MQ=255)
                                                                             aGGTCACCATAAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCGGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGc                                                                       <  2:45010/150‑1 (MQ=255)
                                                                             aGGTCACCAGGAAACAGGTCGCTGAGGTTGCGGTAAAACGGTACATACATTGGCTGTTCCAGTCCCAGTTGCTGGCGAATTTTTTCCACAATAGCTTCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCg                                                                      <  2:666211/151‑1 (MQ=255)
                                                                                   ccagaaaaagggTCGCTGACGTGGGGGAAAAACTGTACATACAGTGGCTGGTCCGGCCCCAGTGGCGGGCAAATATTTCCCACGATAGCATCGCGGGCGGGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGata                                                                 <  1:717279/140‑1 (MQ=255)
                                                                                       aaaaCAGGCCGCTGACGTGGCGGTAAAACTGTACATACAGTGGCTGGTCCAGCCCCGGTTGCGGGCGAGTATTTCCCACGATAGCATCGCTGGCGCGTTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATaaa                                                            <  2:798544/150‑1 (MQ=255)
                                                                                        aagcAGGCCGTTGAGGTGGCGATAAAACTGTCCATACAGTGGCTGGTCCGGCCCCAGTTGCGGGCGAATTTTTCCCACGATGCCACCGCTGGGGCGGCCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAAc                                                           <  1:17323/148‑1 (MQ=255)
                                                                                         aaCAGGTCGCTGTCGTGGTGGAAAACCTGTAAATACTGTGTCTTGCCCAGTCCCAGTTTTGGGCGATTTTTTCCCGTTATAGCTTCGTTGGCGCGTTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACg                                                          <  2:754147/151‑1 (MQ=255)
                                                                                                     aCGTGGCGGTAAAACTGTACATACGGGGGCTGGCCCAGCCCCGGTTGCGGGCAAATATTTCCCACGATAGCATCGCGGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCa                                                <  1:829740/149‑1 (MQ=255)
                                                                                                     aCGGAGCGAAAAAACTTTACATACAGTGGTGGGTCCAGTCCCAGTTGCGGGCAAATATTTCCCAGGATAGCATCGCGGGCGGGTACACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCa                                                    <  1:726987/145‑1 (MQ=255)
                                                                                                          gCGGAAAAACTGTACAAACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTCCCACGATAGCATCGCTGGCGCGGCCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAcgccc                                         <  2:51720/151‑1 (MQ=255)
                                                                                                               aaacCTGTACATACGGGGGCTGGCCCAGCCCCAGTTGGGGGCAAATATTTTCCACGAAGCCATCGCGGGCGCGTCCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAcgccc                                         <  1:883100/146‑1 (MQ=255)
                                                                                                               aaaaCTGTAAATACAGTGGCGGGTCCAGCCCCGGTTGCTGGCAATTTTTTTCCCGGATAGCTCCGGTGGCGGGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAc                                             <  2:341346/142‑1 (MQ=255)
                                                                                                                         atacaGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCAcc                          >  2:391480/1‑151 (MQ=255)
                                                                                                                           acaGTGGCTGGCCCACCCCCGGTGGAGGGCAAATTTTTCCCAGGATAGCACCGCGGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAg                        <  1:481893/151‑1 (MQ=255)
                                                                                                                                  cTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGccccc                   >  1:545163/1‑149 (MQ=255)
                                                                                                                                       ccATTCCCAGTTGTGGGCAAATTTTTCCCACGATAGCATCGCGGCCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCcagca               <  2:693563/148‑1 (MQ=255)
                                                                                                                                             ccAGTTGCTGGCGAATATTTTCCACGATAGCCTCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACCAGCACCAGCCCCCAGCAGCGTTGGc       >  1:322187/1‑150 (MQ=255)
                                                                                                                                                  tGCTGGCGATTTTTTCCCACGATGGCATCGCGGGCGGGGCCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTaaa  <  2:417206/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAA  >  NC_000913/1560617‑1560912

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