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_000913	1,068,285	0	T	.	100.0%	175.9 / NA	42	intergenic (+31/+226)	ymdF/rutG	KGG family protein/pyrimidine permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (21/21);  total (21/21)

CGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTG  >  NC_000913/1068167‑1068413                                                                                                                       |                                                                                                                                 cGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCAt                                                                                                            >  2:939371/1‑139 (MQ=255)      gTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCAccg                                                                                                       >  1:186386/1‑139 (MQ=255)         aTTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGcc                                                                                                              <  2:559719/129‑1 (MQ=255)         aTTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGcc                                                                                                              >  1:559719/1‑129 (MQ=255)         aTTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCAc                                                                                                         >  2:528495/1‑134 (MQ=255)         aTTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCAc                                                                                                         <  1:528495/134‑1 (MQ=255)            tCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGtt                                                                                                 <  1:808119/139‑1 (MQ=255)                 aaTGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                            <  2:186386/139‑1 (MQ=255)                     aCCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtgg                                                                                        >  1:511027/1‑139 (MQ=255)                     aCCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                            >  2:1003808/1‑135 (MQ=255)                     aCCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                            <  1:1003808/135‑1 (MQ=255)                           aGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTGTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATa                                                                                  <  1:428625/139‑1 (MQ=255)                                       aaGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACTACT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAAt                                                                      >  2:648417/1‑139 (MQ=255)                                                    aGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                            <  1:10342/104‑1 (MQ=255)                                                    aGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                            >  2:10342/1‑104 (MQ=255)                                                         gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                    >  2:692378/1‑139 (MQ=255)                                                          tAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcgc                                                   >  1:130644/1‑139 (MQ=255)                                                             gagCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtggt                                                                                       >  2:5591/1‑100 (MQ=255)                                                             gagCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtggt                                                                                       <  1:5591/100‑1 (MQ=255)                                                             gagCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGAt                                                <  2:130644/139‑1 (MQ=255)                                                              agCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATg                                               <  1:648417/139‑1 (MQ=255)                                                                 aGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCtt                                            >  2:64319/1‑139 (MQ=255)                                                                   tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                    <  1:781198/129‑1 (MQ=255)                                                                   tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                    >  2:781198/1‑129 (MQ=255)                                                                   tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTAt                                          >  1:705137/1‑139 (MQ=255)                                                                    cACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtgatga                                                                           <  1:195711/105‑1 (MQ=255)                                                                    cACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtgatga                                                                           >  2:195711/1‑105 (MQ=255)                                                                       ggCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg                                      <  2:695031/139‑1 (MQ=255)                                                                             aGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg                                      <  2:591639/133‑1 (MQ=255)                                                                             aGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGCTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg                                      >  1:591639/1‑133 (MQ=255)                                                                               cGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGCTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCag                              <  1:692378/139‑1 (MQ=39)                                                                                                       gACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGAGACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGc      <  2:705137/139‑1 (MQ=255)                                                                                                           gaatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACt   <  1:939371/135‑1 (MQ=255)                                                                                                            aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATa                                                                                  >  2:1014226/4‑59 (MQ=39)                                                                                                            aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATa                                                                                  <  1:1014226/56‑1 (MQ=39)                                                                                                            aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTg  >  1:661390/4‑139 (MQ=255)                                                                                                             atcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCa                                        >  1:665837/3‑100 (MQ=255)                                                                                                             atcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCa                     <  2:308478/117‑1 (MQ=255)                                                                                                             atcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCa                     >  1:308478/3‑119 (MQ=255)                                                                                                             ataaTTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCa                                        <  2:665837/97‑1 (MQ=255)                                                                                                              tcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGa                                                                                    <  2:617594/54‑1 (MQ=39)                                                                                                              tcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGa                                                                                    >  1:617594/2‑55 (MQ=39)                                                                                                                       |                                                                                                                                 CGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTG  >  NC_000913/1068167‑1068413

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