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%	232.8 / NA	57	coding (249/276 nt)	yjfY	YhcN family protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (30/27);  total (30/27)

CGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCAC  >  NC_000913/4424414‑4424674                                                                                                                                  |                                                                                                                                    cGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGt                                                                                                                           <  1:734695/139‑1 (MQ=255)    gagaTAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTAc                                                                                                                        <  1:1244993/139‑1 (MQ=255)        tAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTccc                                                                                                                    <  1:109107/139‑1 (MQ=255)         aaCTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTccc                                                                                                                    <  2:1555023/138‑1 (MQ=255)         aaCTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTccc                                                                                                                    >  1:1555023/1‑138 (MQ=255)              gTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTc                                                                                                              >  1:33084/1‑139 (MQ=255)               tCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTcc                                                                                                             <  2:978128/139‑1 (MQ=255)                 tAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCg                                                                                                           >  1:623912/1‑139 (MQ=255)                 tAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCg                                                                                                           >  2:798292/1‑139 (MQ=255)                       cGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAg                                                                                                     <  1:798292/139‑1 (MQ=255)                           gCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAAtt                                                                                                 <  1:492037/139‑1 (MQ=255)                            ccGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTa                                                                                                <  2:33084/139‑1 (MQ=255)                              gTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                              <  1:1484607/139‑1 (MQ=255)                                      aaaaaaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGgtagta                                                                                      >  2:1268848/1‑139 (MQ=255)                                             gCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTgcg                                                                               <  2:1350583/139‑1 (MQ=255)                                                 ccTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTcc                                                                                                             <  1:841222/105‑1 (MQ=255)                                                 ccTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTcc                                                                                                             >  2:841222/1‑105 (MQ=255)                                                             gcaaagcaaaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGttt                                                               >  2:538794/1‑139 (MQ=255)                                                                             cAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTg                                               <  2:360877/139‑1 (MQ=255)                                                                             cAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTg                                               >  1:1093077/1‑139 (MQ=255)                                                                                       ttGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttct                                                             >  2:1144870/1‑115 (MQ=255)                                                                                       ttGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttct                                                             <  1:1144870/115‑1 (MQ=255)                                                                                           ttATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCttt                                 <  2:1093077/139‑1 (MQ=255)                                                                                           ttATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCttt                                 <  2:1532523/139‑1 (MQ=255)                                                                                            tataCCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTg                                                                          >  1:46647/1‑97 (MQ=255)                                                                                            tataCCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTg                                                                          <  2:46647/97‑1 (MQ=255)                                                                                             ataCCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCa                               >  2:242524/1‑139 (MQ=255)                                                                                             ataCCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCa                               >  2:467493/1‑139 (MQ=255)                                                                                             ataCCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCa                               <  1:1215164/139‑1 (MQ=255)                                                                                                   ggTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa                                                     <  2:461640/111‑1 (MQ=255)                                                                                                   ggTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa                                                     >  1:461640/1‑111 (MQ=255)                                                                                                    gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGaa                                                                                                   <  2:1150541/60‑1 (MQ=255)                                                                                                    gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGaa                                                                                                   >  1:1150541/5‑64 (MQ=255)                                                                                                    gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGGATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTTAGTGGCGa                        >  2:956469/5‑139 (MQ=255)                                                                                                    gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCt                                                                   >  1:887584/5‑96 (MQ=255)                                                                                                    gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGAc                                             <  1:247261/114‑1 (MQ=255)                                                                                                    gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGc                          >  1:552980/5‑137 (MQ=255)                                                                                                    gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGa                        >  1:881514/5‑139 (MQ=255)                                                                                                    gtgtTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCt                                                                   <  2:887584/96‑1 (MQ=255)                                                                                                    gtgtTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGAc                                             >  2:247261/1‑118 (MQ=255)                                                                                                    gtgtTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGc                          <  2:552980/137‑1 (MQ=255)                                                                                                     tgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGaa                       >  1:843539/4‑139 (MQ=255)                                                                                                      gttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                              >  1:46432/3‑67 (MQ=255)                                                                                                      gttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                              <  2:46432/65‑1 (MQ=255)                                                                                                         aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTa                   >  2:156406/1‑139 (MQ=255)                                                                                                         aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTa                   >  1:358304/1‑139 (MQ=255)                                                                                                         aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTa                   >  2:1239612/1‑139 (MQ=255)                                                                                                         aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTAGCGAAATTa                   >  1:1105480/1‑139 (MQ=255)                                                                                                                      cgcgCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttctt                                                            >  2:367506/1‑85 (MQ=255)                                                                                                                      cgcgCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttctt                                                            <  1:367506/85‑1 (MQ=255)                                                                                                                       gcgcTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt     <  2:1105480/139‑1 (MQ=255)                                                                                                                        cgcTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCttt                                 <  1:559171/110‑1 (MQ=255)                                                                                                                        cgcTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCttt                                 >  2:559171/1‑110 (MQ=255)                                                                                                                        cgcTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATc    >  2:1117872/1‑139 (MQ=255)                                                                                                                          cTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttctt                                                            >  1:1593210/1‑81 (MQ=255)                                                                                                                          cTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttctt                                                            <  2:1593210/81‑1 (MQ=255)                                                                                                                          cTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCAc  >  1:255650/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                    CGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCAC  >  NC_000913/4424414‑4424674

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