Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,068,285 Δ1 bp intergenic (+31/+226) ymdF → / ← rutG KGG family protein/pyrimidine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,068,2850T.100.0% 188.0 / NA 43intergenic (+31/+226)ymdF/rutGKGG family protein/pyrimidine permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (18/25);  total (18/25)

GTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACT  >  NC_000913/1068154‑1068412
                                                                                                                                   |                                                                                                                               
gTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTc                                                                                                                        <  2:499889/139‑1 (MQ=255)
 tGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCg                                                                                                                       >  2:301128/1‑139 (MQ=255)
             cGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCAt                                                                                                           >  1:503402/1‑139 (MQ=255)
                           aggaTGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGTTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                                             <  2:701015/136‑1 (MQ=255)
                             aaTGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGGGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           <  2:665543/139‑1 (MQ=255)
                             aaTGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           <  2:503402/139‑1 (MQ=255)
                             aaTGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           <  1:301128/139‑1 (MQ=255)
                                      cAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGAt                                                                                  >  1:568090/1‑139 (MQ=255)
                                             cATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           <  1:209089/123‑1 (MQ=255)
                                             cATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           >  2:209089/1‑123 (MQ=255)
                                                        ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGc                                                                                                                          >  1:194058/1‑81 (MQ=255)
                                                        ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGc                                                                                                                          <  2:194058/81‑1 (MQ=255)
                                                        ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAg                                                                >  1:777980/1‑139 (MQ=255)
                                                          tAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                                         <  1:771439/128‑1 (MQ=255)
                                                          tAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                                         >  2:771439/1‑128 (MQ=255)
                                                          tAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCt                                                              >  1:794001/1‑139 (MQ=255)
                                                             aaaaGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGc                                                           <  2:663246/139‑1 (MQ=255)
                                                                          agCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATg                                              <  1:532460/139‑1 (MQ=255)
                                                                               tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           >  2:629671/1‑89 (MQ=255)
                                                                               tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           <  1:629671/89‑1 (MQ=255)
                                                                               tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           <  1:458599/89‑1 (MQ=255)
                                                                               tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           >  2:458599/1‑89 (MQ=255)
                                                                               tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTAt                                         >  1:469005/1‑139 (MQ=255)
                                                                                cACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATc                                        <  2:777980/139‑1 (MQ=255)
                                                                                        aaGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTAc                                <  2:631458/139‑1 (MQ=255)
                                                                                        aaGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTAc                                <  2:568090/139‑1 (MQ=255)
                                                                                           cGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                   <  2:231080/117‑1 (MQ=255)
                                                                                           cGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                   >  1:231080/1‑117 (MQ=255)
                                                                                               aaCTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCagaaga                         <  2:365320/139‑1 (MQ=255)
                                                                                                    accGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTTAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                                         <  2:659535/85‑1 (MQ=255)
                                                                                                    accGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCAc                                                                                                        >  2:696759/2‑55 (MQ=255)
                                                                                                    accGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCAc                                                                                                        <  1:696759/54‑1 (MQ=255)
                                                                                                    accGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                                         >  1:659535/2‑86 (MQ=255)
                                                                                                    accGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCa                    <  1:623440/138‑1 (MQ=255)
                                                                                                                  tGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGcc                                   >  2:872242/1‑110 (MQ=255)
                                                                                                                  tGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGcc                                   <  1:872242/110‑1 (MQ=255)
                                                                                                                   gACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGt                                                                   >  2:110983/1‑77 (MQ=255)
                                                                                                                   gACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGt                                                                   <  1:110983/77‑1 (MQ=255)
                                                                                                                   gACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGc     >  1:489693/1‑139 (MQ=255)
                                                                                                                    aCGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCa    <  1:606521/139‑1 (MQ=255)
                                                                                                                       gaatcATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACt  <  2:489693/135‑1 (MQ=255)
                                                                                                                         atcATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCa                                       >  2:319144/3‑100 (MQ=255)
                                                                                                                         atcATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCa                                       <  1:319144/98‑1 (MQ=255)
                                                                                                                                   |                                                                                                                               
GTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACT  >  NC_000913/1068154‑1068412

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