Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,211,652	C→T	intergenic (+75/+28)	icdC → / ← iraM	pseudogene, isocitrate dehydrogenase C‑terminal gene fragment, idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion/RpoS stabilzer during Mg starvation, anti‑RssB factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,211,652	0	C	T	100.0%	59.7 / NA	20	intergenic (+75/+28)	icdC/iraM	pseudogene, isocitrate dehydrogenase C‑terminal gene fragment, idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion/RpoS stabilzer during Mg starvation, anti‑RssB factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/13);  total (7/13)

CGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTGCTT  >  NC_000913/1211523‑1211776                                                                                                                                  |                                                                                                                             cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGc                                                                                                                     <  2:733396/139‑1 (MQ=255) cgcTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGc                                                                                                                     >  1:797488/1‑139 (MQ=255)            tgAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATg                                                                                                          <  2:797488/139‑1 (MQ=255)               aaTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAg                                                                                                       <  2:706195/139‑1 (MQ=255)                         gTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTAAAATATTAAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACtt                                                                                             <  2:588242/139‑1 (MQ=255)                             ggTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCACATTACTTGAGc                                                                                         >  2:5225/1‑139 (MQ=255)                                  aGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCata                                                                                    <  1:691466/139‑1 (MQ=255)                                           cGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCatat                                                                           <  1:67305/139‑1 (MQ=255)                                                        tctctcCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCaa                                                              <  2:640020/139‑1 (MQ=255)                                                             ccATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACtatt                                                         <  2:301813/139‑1 (MQ=255)                                                                                   ggCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAg                                                                                                                      >  1:184266/1‑56 (MQ=255)                                                                                   ggCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAg                                                                                                                      <  2:184266/56‑1 (MQ=255)                                                                                   ggCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTAtt                                   >  2:628308/1‑139 (MQ=255)                                                                                   ggCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTAtt                                   <  1:628308/139‑1 (MQ=255)                                                                                   ggCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTAtt                                   >  1:382345/1‑139 (MQ=255)                                                                                    gCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCaa                                                              <  1:776194/111‑1 (MQ=255)                                                                                    gCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCaa                                                              >  2:776194/1‑111 (MQ=255)                                                                                              cGCCCGTTGTTTTATTTATGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTcc                        >  1:273356/1‑139 (MQ=255)                                                                                                              tatGTGGATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGa        <  2:382345/139‑1 (MQ=255)                                                                                                                    gATATTATTAATAGTATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTGCtt  <  1:211373/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                             CGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGTTTTATTTATGTGGATATTATTAATAGCATATCGAGCATATTTATATGAAGCCCATTACTTGAGCCCATATGGGCATATTTTTATAATGCAACTATTATGTAAACATTTATTTGTTATTTTGCTTTCTCCTGGAGGACACTCTTGACTGCTT  >  NC_000913/1211523‑1211776

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