Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,559,325 C→A 21.6% A196S (GCG→TCG)  ddpC ← putative D,D‑dipeptide ABC transporter membrane subunit DdpC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,559,3250CA21.6% 37.1 / 10.3 28A196S (GCG→TCG) ddpCputative D,D‑dipeptide ABC transporter membrane subunit DdpC
Reads supporting (aligned to +/- strand):  ref base C (7/15);  new base A (4/2);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-01

GTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGC  >  NC_000913/1559190‑1559466
                                                                                                                                       |                                                                                                                                             
gtttcgCACCCAGGCCAATAAGTGCTGTTGTGGCGGCCATTAGAACCGCGCTAGCGATGCCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATAGGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAg                                                                                                                                        <  2:650754/138‑1 (MQ=255)
 ttgAGCACCCACGCCACTAACTCCCGACGTCTCGGCCATTAAAATCGCGCTACCGATATCCAGCGATTCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAgg                                                                                                                                       <  2:164861/143‑1 (MQ=255)
   gAGCCCCCAGGCCCATAGATCCCCACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAGGtt                                                                                                                                     <  2:154749/143‑1 (MQ=255)
      cACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGGCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAGGTTTTc                                                                                                                                  <  1:427279/143‑1 (MQ=255)
       aCCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGAGCCAAAGGTTTTcg                                                                                                                                 >  1:592416/1‑143 (MQ=255)
       aCCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGAGCCAAAGGTTTTcg                                                                                                                                 >  1:594108/1‑143 (MQ=255)
          cAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCCGCGCTGCCTGCACGGTCCGCGGCGGTAGGGAGTTACGTCAAATATGCCAGTTGATCAGACGCCACCGCGACGCGCCCACGGCTTccgccg                                                                                                                              >  2:625654/1‑143 (MQ=255)
                aaTAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGCGCCGAAGGCTTTCGCCCCCGGaa                                                                                                                        >  2:220783/1‑143 (MQ=255)
                aaTAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGAGCCAAAGGTTTTCGCCGCCTGaa                                                                                                                        >  1:222765/1‑143 (MQ=255)
                      tCCCAACGTGGCGGCCATTAAAATCGCGCTACCGGTATCCCGCGATGCCTGCACGAGCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAgg                                                                                                                  <  2:404346/143‑1 (MQ=255)
                                   gctaTTCAAATCGCGCTACCGATCTCCACCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATACGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTAc                                                                                                     <  1:457847/140‑1 (MQ=255)
                                     cATTAAAATCGCGCTACCGATCTCCAGCGATGCCTGCACGAGCAGCGGCGGTAGGGAGTTACGTAAAATAGGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACaa                                                                                                   <  2:231237/143‑1 (MQ=255)
                                         aaaaTGGCGCTCCAGTGATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGAAGTTACGTAACATATGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTaa                                                                                               <  2:294445/143‑1 (MQ=255)
                                            aTCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGCTCCGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGGACATAGGTATACTGGCGTACAACTAATGc                                                                                            >  1:651417/1‑143 (MQ=255)
                                             tCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAACATATGCCAGTTGATCAGATGCCAACGAGACGCGCCCAAGGTCTTCGGCGCCCGCACATAGGTACACTGGCGCCCAACCAATGcc                                                                                           >  2:570368/1‑142 (MQ=255)
                                                                       tGCACGATCAGCGGCGGTAGGGCGTTACGTAAAATATGCCAGTTGATCAGCTGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAAAAAAGGGATCCTGGCGTACCACTAATGCTTCACCCCGCGCCCGGCACACATcaaa                                                                 >  2:445092/1‑143 (MQ=255)
                                                                           cGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGAGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGa                                                             <  1:658138/143‑1 (MQ=255)
                                                                                 gcggcgGTAGGCAGTTACGTAAAATATGCCAGTTGATCAGAGGCCAACGGTACGCGCCAAAGGGTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGc                                                       <  2:77444/143‑1 (MQ=255)
                                                                                     cggTAGGGAGTTACGTAAAATATGCCAGTTGATCAGCTGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGCACATAGGTATACTGGCGTACGACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCAcaa                                                   >  1:176434/1‑143 (MQ=255)
                                                                                           ggAGTTACGTAAAATATGCGAGTTGATCAGATGCCAACGAGACGAGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCAcaatagcaa                                             <  1:505082/143‑1 (MQ=255)
                                                                                                 aCGTAAAATATGCCAGTGGATCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCa                                       <  2:171675/143‑1 (MQ=255)
                                                                                                  cGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAg                                      <  1:567384/143‑1 (MQ=255)
                                                                                                           tGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAGGACTTCGCCGCCTGAACATAGCTATACCGGCGTACAACTAATGCTCGCCCCCGCGCCAGGCGGACAGGCAAGAGAACACGCACCACCGCCAGAGCCAGCACggggca                             >  2:64465/1‑139 (MQ=255)
                                                                                                            gCCTGTGGATCAGAGGCCGACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTa                            <  2:71254/143‑1 (MQ=255)
                                                                                                            gCCAGTTGATCAGATGCCAACGAGACGAGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAAAGAGCATTCTTCGGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATCGCCAGCATGGCGGTa                            >  1:494796/1‑143 (MQ=255)
                                                                                                                   gATCAGAGGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTaaacaaac                     <  2:513136/143‑1 (MQ=255)
                                                                                                                     tccgATGCCACCGAGACGCGCCAAAGGTGTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTc                   <  1:635372/140‑1 (MQ=255)
                                                                                                                       cgagtccAAGGATACGCGCCAAAGTCTGTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCgg                 <  1:668954/138‑1 (MQ=255)
                                                                                                                                      gcgcCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATCGCCAGCATGGCGTTACACAAACTCGGCCCGAGACCGCCAGc  >  2:124802/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                             
GTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAACTCGGCCCGAGAGCGGCAGC  >  NC_000913/1559190‑1559466

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