Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,424,543	Δ1 bp	coding (249/276 nt)	yjfY ←	YhcN family protein, periplasmic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,424,543	0	C	.	100.0%	158.0 / NA	38	coding (249/276 nt)	yjfY	YhcN family protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (20/18);  total (20/18)

TCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAT  >  NC_000913/4424410‑4424671                                                                                                                                      |                                                                                                                                 tcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTg                                                                                                                            <  1:164354/139‑1 (MQ=255)        gagaTAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTAc                                                                                                                     <  2:723931/139‑1 (MQ=255)                             cAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGt                                                                                                             <  2:562646/126‑1 (MQ=255)                             cAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGt                                                                                                             >  1:562646/1‑126 (MQ=255)                             cAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGaa                                                                                                >  2:137623/1‑139 (MQ=255)                                ccGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTa                                                                                             <  2:562305/139‑1 (MQ=255)                                  gTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                           <  1:137623/139‑1 (MQ=255)                                  gTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                           <  2:234586/139‑1 (MQ=255)                                  gTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                           <  2:140565/139‑1 (MQ=255)                                                         cGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgcgc                                                                    <  2:618929/139‑1 (MQ=255)                                                                   aaagcaaaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttct                                                          >  1:188055/1‑139 (MQ=255)                                                                     agcaaaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttcttc                                                        <  2:188055/139‑1 (MQ=255)                                                                          aCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTa                                                                                             >  2:754194/1‑97 (MQ=255)                                                                          aCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTa                                                                                             <  1:754194/97‑1 (MQ=255)                                                                             aTGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATg                                                >  1:58412/1‑139 (MQ=255)                                                                              tGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGc                                               <  1:513182/139‑1 (MQ=255)                                                                               ggCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgcg                                                                     >  1:199376/1‑116 (MQ=255)                                                                               ggCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgcg                                                                     <  2:199376/116‑1 (MQ=255)                                                                                          gTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTc                                   >  2:636515/1‑139 (MQ=255)                                                                                                tataCCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTc                             >  2:103555/1‑139 (MQ=255)                                                                                                  taCCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAg                           <  1:636515/139‑1 (MQ=255)                                                                                                        gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGa                     >  1:701020/5‑139 (MQ=255)                                                                                                        gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGa                     >  2:139139/5‑139 (MQ=255)                                                                                                        gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGa                     >  2:208645/5‑139 (MQ=255)                                                                                                        gtgtTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGa                     >  2:648204/1‑139 (MQ=255)                                                                                                         tgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                           <  2:388496/65‑1 (MQ=255)                                                                                                         tgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                           >  1:388496/4‑68 (MQ=255)                                                                                                         tgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGaa                    >  1:55987/4‑139 (MQ=255)                                                                                                         tgttAGCGATAAATATCCGCGCTGGCGCG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGaa                    >  1:140671/4‑139 (MQ=255)                                                                                                             aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGc                       <  2:430205/132‑1 (MQ=255)                                                                                                             aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAAtt                 <  1:654782/138‑1 (MQ=255)                                                                                                             aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAAtt                 >  2:654782/1‑138 (MQ=255)                                                                                                             aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTa                >  2:571997/1‑139 (MQ=255)                                                                                                             aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGCGCTTCGCTTTCAGTGGc                       >  1:430205/1‑132 (MQ=255)                                                                                                                 aTAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtgat            >  1:573488/1‑139 (MQ=255)                                                                                                                           gcgcTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt  <  1:203469/139‑1 (MQ=255)                                                                                                                             gcTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGgtagt                                                                                    <  2:718004/55‑1 (MQ=255)                                                                                                                             gcTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGgtagt                                                                                    >  1:718004/1‑55 (MQ=255)                                                                                                                                      |                                                                                                                                 TCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAT  >  NC_000913/4424410‑4424671

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