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_000913	1,068,276	0	G	.	100.0%	212.8 / NA	51	intergenic (+22/+235)	ymdF/rutG	KGG family protein/pyrimidine permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (25/26);  total (25/26)

GACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCA  >  NC_000913/1068157‑1068410                                                                                                                        |                                                                                                                                       gACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTc                                                                                                                   >  2:944526/1‑139 (MQ=255)           cGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCAt                                                                                                         >  1:80614/1‑139 (MQ=255)           cGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCAt                                                                                                         >  2:486985/1‑139 (MQ=255)                           aaTGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                         <  2:80614/139‑1 (MQ=255)                           aaTGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                         <  1:944526/139‑1 (MQ=255)                           aaTGACCCGCAGCGCGCATCTGAAGCAGATAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                         <  2:124399/139‑1 (MQ=255)                               aCCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtgg                                                                                     >  1:1069430/1‑139 (MQ=255)                               aCCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                         >  2:263916/1‑135 (MQ=255)                               aCCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                         <  1:263916/135‑1 (MQ=255)                                cccGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGGCGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGt                                                                                               <  1:1208778/128‑1 (MQ=255)                                cccGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGGCGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGt                                                                                               >  2:1208778/1‑128 (MQ=255)                                cccGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtggt                                                                                    <  2:556740/139‑1 (MQ=255)                                         cgcATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                         <  1:229992/125‑1 (MQ=255)                                         cgcATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                         >  2:229992/1‑125 (MQ=255)                                         cgcATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtga                                                                           >  1:1067826/1‑139 (MQ=255)                                           cATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtgatg                                                                         >  1:503138/1‑139 (MQ=255)                                              cTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACa                                                                      >  1:1123045/1‑139 (MQ=255)                                                        tAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACATATGTAAGCt                                                            <  2:817645/139‑1 (MQ=255)                                                        tAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCt                                                            >  1:1006524/1‑139 (MQ=255)                                                                 tggtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATg                                                   >  2:694417/1‑139 (MQ=255)                                                                 tggtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATg                                                   >  2:231709/1‑139 (MQ=255)                                                                   gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                 <  2:875603/139‑1 (MQ=255)                                                                   gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                 <  2:1006524/139‑1 (MQ=255)                                                                   gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                 <  2:1073775/139‑1 (MQ=255)                                                                          cAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCt                                          >  2:80244/1‑139 (MQ=255)                                                                             tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAg                                                                                                                         <  1:531523/57‑1 (MQ=255)                                                                             tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAg                                                                                                                         >  2:531523/1‑57 (MQ=255)                                                                             tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCaaa                                                                                                                <  2:671565/66‑1 (MQ=255)                                                                             tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCaaa                                                                                                                >  1:671565/1‑66 (MQ=255)                                                                             tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTAt                                       >  2:234731/1‑139 (MQ=255)                                                                             tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTAt                                       >  1:1001145/1‑139 (MQ=255)                                                                                 ggCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg                                   <  1:80244/139‑1 (MQ=255)                                                                                 ggCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg                                   <  1:234731/139‑1 (MQ=255)                                                                                 ggCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg                                   <  1:231709/139‑1 (MQ=255)                                                                                 ggCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg                                   <  2:1001145/139‑1 (MQ=255)                                                                                  gCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCAt                                                                                                         >  2:1213178/1‑68 (MQ=255)                                                                                  gCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCAt                                                                                                         <  1:1213178/68‑1 (MQ=255)                                                                                      aaGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTAc                              <  1:486985/139‑1 (MQ=255)                                                                                       aGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTAcc                             >  1:255208/1‑139 (MQ=255)                                                                                       aGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTAcc                             >  1:1193545/1‑139 (MQ=255)                                                                                       aGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTAcc                             >  2:505235/1‑139 (MQ=255)                                                                                       aGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTAcc                             >  1:990809/1‑139 (MQ=255)                                                                                        gCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCa                            <  2:1193545/139‑1 (MQ=255)                                                                                           aCAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                         >  1:144901/1‑75 (MQ=255)                                                                                           aCAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                         <  2:144901/75‑1 (MQ=255)                                                                                             aaCTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCagaaga                       <  1:505235/139‑1 (MQ=255)                                                                                                  accGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCa                  <  2:592231/138‑1 (MQ=255)                                                                                                       gCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                         <  2:586262/63‑1 (MQ=255)                                                                                                       gCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                         >  1:586262/1‑63 (MQ=255)                                                                                                                  aCGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCa  <  2:1069430/139‑1 (MQ=255)                                                                                                                  aCGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCa  <  2:255208/139‑1 (MQ=255)                                                                                                                        |                                                                                                                                       GACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCA  >  NC_000913/1068157‑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: ‑