Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 469,207 Δ1 bp coding (337/1773 nt) mdlA → ABC transporter family protein MdlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913469,2070A.90.0% 59.7 / 4.4 20coding (337/1773 nt)mdlAABC transporter family protein MdlA
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base . (8/10);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

GTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAAT  >  NC_000913/469072‑469337
                                                                                                                                       |                                                                                                                                  
gTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGt                                                                                                                               <  1:330412/140‑1 (MQ=255)
   gtgGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGTGAc                                                                                                                            <  1:209425/140‑1 (MQ=255)
    tgGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGt                                                                                                                               <  2:96096/136‑1 (MQ=255)
    tgGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGt                                                                                                                               >  1:96096/1‑136 (MQ=255)
                 gCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGTGACCTCATGGCTCGTGc                                                                                                              >  2:154985/1‑140 (MQ=255)
                                               tGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTgc                                                                                <  1:283154/140‑1 (MQ=255)
                                               tGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTgc                                                                                <  1:253330/140‑1 (MQ=255)
                                               tGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTgc                                                                                >  2:253330/1‑140 (MQ=255)
                                               tGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTgc                                                                                >  2:283154/1‑140 (MQ=255)
                                                 cGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTgccg                                                                              >  2:172749/1‑140 (MQ=255)
                                                 cGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTgccg                                                                              >  2:309176/1‑140 (MQ=255)
                                                               gCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAgggggg                                                                >  1:257418/1‑138 (MQ=255)
                                                                            gTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCtggtgg                                                   >  2:51417/1‑140 (MQ=255)
                                                                                         ccGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGAt                                      <  1:309176/140‑1 (MQ=255)
                                                                                         ccGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGAt                                      <  1:172749/140‑1 (MQ=255)
                                                                                            tCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATggg                                   <  2:145385/140‑1 (MQ=255)
                                                                                                  aaGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTgcgc                             <  1:51417/140‑1 (MQ=255)
                                                                                                              tCCTGAGTTTTACCTGCGTCATCGCCCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTat                 <  2:231358/140‑1 (MQ=255)
                                                                                                                   aGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTatgatg             >  2:156979/1‑140 (MQ=255)
                                                                                                                              cGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGGGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAAt  >  1:71674/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                  
GTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAAT  >  NC_000913/469072‑469337

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