Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,068,276 Δ1 bp intergenic (+22/+235) 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,2760G.100.0% 181.7 / NA 43intergenic (+22/+235)ymdF/rutGKGG family protein/pyrimidine permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (20/23);  total (20/23)

GTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCA  >  NC_000913/1068154‑1068410
                                                                                                                          |                                                                                                                                      
gTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTc                                                                                                                      >  1:893527/1‑139 (MQ=255)
           agcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGcc                                                                                                           >  1:497479/1‑139 (MQ=255)
             cGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCAt                                                                                                         >  2:568056/1‑139 (MQ=255)
               ggggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCa                                                                                                       <  2:497479/139‑1 (MQ=255)
                    aaTTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCAccgcc                                                                                                  <  2:162735/139‑1 (MQ=255)
                     aTTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCAc                                                                                                      >  2:744869/1‑134 (MQ=255)
                     aTTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCAc                                                                                                      <  1:744869/134‑1 (MQ=255)
                             aaTGACCCGCAGCGCGCATCTGGAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                         <  2:121407/139‑1 (MQ=255)
                             aaTGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                         <  2:680962/139‑1 (MQ=255)
                             aaTGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                         <  1:568056/139‑1 (MQ=255)
                                         cgcgcATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACt                                                                             >  1:732487/1‑139 (MQ=255)
                                               tCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                                       <  2:698072/139‑1 (MQ=255)
                                                        ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCaaa                                                                                                                >  1:588978/1‑89 (MQ=255)
                                                        ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCaaa                                                                                                                <  2:588978/89‑1 (MQ=255)
                                                        ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCAcc                                                                                                     <  1:790470/100‑1 (MQ=255)
                                                        ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCAcc                                                                                                     >  2:790470/1‑100 (MQ=255)
                                                          tAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtgatga                                                                        >  2:397318/1‑127 (MQ=255)
                                                          tAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtgatga                                                                        <  1:397318/127‑1 (MQ=255)
                                                                    ggtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCaaaaa                                                                                                              >  2:236882/1‑79 (MQ=255)
                                                                    ggtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCaaaaa                                                                                                              <  1:236882/79‑1 (MQ=255)
                                                                    ggtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgc                                                  >  2:511244/1‑139 (MQ=255)
                                                                     gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                 <  2:893527/139‑1 (MQ=255)
                                                                      tAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                 <  1:108404/138‑1 (MQ=255)
                                                                      tAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                 >  2:108404/1‑138 (MQ=255)
                                                                               tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTCGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTAt                                       >  2:817709/1‑139 (MQ=255)
                                                                               tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTCGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTAt                                       <  1:817709/139‑1 (MQ=255)
                                                                                   ggCAAAAGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg                                   <  2:732487/139‑1 (MQ=255)
                                                                                         aGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                                       <  1:505611/97‑1 (MQ=255)
                                                                                         aGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGa                                              <  2:207494/122‑1 (MQ=255)
                                                                                         aGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGa                                              >  1:207494/1‑122 (MQ=255)
                                                                                         aGCGACAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCAATCACTGCCGGTTACCCGGTGGTTGATACTGATGAc                                                                       >  2:505611/1‑97 (MQ=255)
                                                                                             aCAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACt                                                                             <  1:790590/87‑1 (MQ=255)
                                                                                             aCAACTAACCGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACt                                                                             >  2:790590/1‑87 (MQ=255)
                                                                                                    accGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCtt                                                           <  1:579469/97‑1 (MQ=255)
                                                                                                    accGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCtt                                                           >  2:579469/2‑98 (MQ=255)
                                                                                                    accGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCa                  <  1:511244/138‑1 (MQ=255)
                                                                                                    accGGGCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCa                  >  2:16709/2‑139 (MQ=255)
                                                                                                        ggCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                 >  1:312139/1‑104 (MQ=255)
                                                                                                        ggCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                 <  2:312139/104‑1 (MQ=255)
                                                                                                         gCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAAt                                                                   >  1:424395/1‑85 (MQ=255)
                                                                                                         gCTAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAAt                                                                   <  2:424395/85‑1 (MQ=255)
                                                                                                           tAATCAATGACGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATAtt           >  2:782870/1‑139 (MQ=255)
                                                                                                                    aCGAATCATTTTTGCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCa  <  1:782870/139‑1 (MQ=255)
                                                                                                                          |                                                                                                                                      
GTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCA  >  NC_000913/1068154‑1068410

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: