Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,714,326 (T)9→8 intergenic (+55/‑51) dtpA → / → gstA dipeptide/tripeptide:H(+) symporter DtpA/glutathione S‑transferase GstA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,714,3180T.95.7% 84.5 / ‑1.2 23intergenic (+47/‑59)dtpA/gstAdipeptide/tripeptide:H(+) symporter DtpA/glutathione S‑transferase GstA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (7/15);  total (7/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

TGATGCTGCTGACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGTGCCTGCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTA  >  NC_000913/1714189‑1714452
                                                                                                                                 |                                                                                                                                      
tGATGCTGCTGACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCttttttttgtt                                                                                                                              <  1:331036/140‑4 (MQ=255)
   tgctgcTGACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAg                                                                                                                          <  1:1063093/140‑1 (MQ=255)
   tgctgcTGACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGAAACCTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAg                                                                                                                          <  2:511827/140‑1 (MQ=255)
      tgctgACCGCGCCGAAACTGCCCCGCATGACGCAGGATGACGCTGCAGACAAAGGGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCaa                                                                                                                       <  1:336485/140‑1 (MQ=255)
       gctgACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAAc                                                                                                                      >  1:571671/1‑140 (MQ=255)
        ctgACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACg                                                                                                                     <  1:696257/140‑1 (MQ=255)
         tgACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACGc                                                                                                                    >  2:744869/1‑140 (MQ=255)
               cgacGAAACTGACCCGCATGCCGCAGGATGACCCGGAAGACAAAGCGGCGAACGCCGCCGTAGGGTAAAATTCAGGGAAACTTTTTTCCAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACGCACTAcc                                                                                                              <  2:346233/137‑1 (MQ=255)
                gcCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTAcc                                                                                                              <  1:1000789/140‑1 (MQ=255)
                            ccGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACca                                                                                                 <  2:28018/140‑1 (MQ=255)
                              gCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGAAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACcaca                                                                                               <  2:784148/140‑1 (MQ=255)
                                  gACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACa                                                                                           >  2:157046/1‑140 (MQ=255)
                                   aCGCAGGATGACGCTGCAGAAAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAACCTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACaa                                                                                          <  2:718957/140‑1 (MQ=255)
                                                                              cGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGTGCCTGct                                               >  1:976374/1‑140 (MQ=255)
                                                                                    aTTCAGGGAACCTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGTGCCTGCTCTCTCg                                         <  1:744869/140‑1 (MQ=255)
                                                                                              aCTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGTGCCTGCTCTCTCGCTTCCCatat                               <  2:976374/140‑1 (MQ=255)
                                                                                                tctTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGTGCCTGCTCTCTCGCTTCCCATATCa                             <  1:833336/140‑1 (MQ=255)
                                                                                                  ttttACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGTGCCTGCTCTCTCGCTTCCCATATCAcc                           <  1:705876/140‑1 (MQ=255)
                                                                                                          gCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGTGCCTGCTCTCTCGCTTc                                     >  2:628921/1‑122 (MQ=255)
                                                                                                          gCCGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGTGCCTGCTCTCTCGCTTc                                     <  1:628921/122‑1 (MQ=255)
                                                                                                            cGCTAACTTTTCGTTAGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGTGCCTGCTCTCTCGCTTCCCATATCACCCTGCGTgaga                 >  1:928633/1‑140 (MQ=255)
                                                                                                                   ttttCGTTAGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGTGCCTGCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGaaa          <  1:157046/140‑1 (MQ=255)
                                                                                                                        gTTAGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGTGCCTGCTCTCTCGCTTCCCATATCACCCTGCGTAAGAGCGGAAAGGAtt     <  1:712842/140‑1 (MQ=255)
                                                                                                                           aGCGGCTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGTGCCTGCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTa  >  1:354857/1‑140 (MQ=255)
                                                                                                                                 |                                                                                                                                      
TGATGCTGCTGACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGTGCCTGCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTA  >  NC_000913/1714189‑1714452

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