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

GTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCAC  >  NC_000913/4424427‑4424674                                                                                                                     |                                                                                                                                    gTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTc                                                                                                              >  1:662953/1‑139 (MQ=255)     aaTTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGc                                                                                                          >  1:1058057/1‑139 (MQ=255)        tCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATc                                                                                                       <  2:287510/139‑1 (MQ=255)            cAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGaa                                                                                                   >  1:930393/1‑139 (MQ=255)                 gTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                              <  1:301211/139‑1 (MQ=255)                    tttCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGTCGTTACTCCCCGGTTCCCGCATCAGAATTACGTgg                                                                                           >  1:476892/1‑139 (MQ=255)                                gCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTgcg                                                                               <  2:930393/139‑1 (MQ=255)                                           aGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATc                                                                    <  1:1145530/139‑1 (MQ=255)                                           aGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATc                                                                    <  1:33383/139‑1 (MQ=255)                                                   aagcaaaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttctt                                                            <  2:578199/139‑1 (MQ=255)                                                        aaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCAttt                                                       <  2:1058057/139‑1 (MQ=255)                                                            aTGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgcg                                                                        <  1:1213536/118‑1 (MQ=255)                                                            aTGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgcg                                                                        >  2:1213536/1‑118 (MQ=255)                                                              ggCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCAt                                                         >  1:636649/1‑131 (MQ=255)                                                              ggCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCAt                                                         <  2:636649/131‑1 (MQ=255)                                                                  aTGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGAc                                             >  1:545142/1‑139 (MQ=255)                                                                              ttATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCttt                                 >  2:827362/1‑139 (MQ=255)                                                                                       gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCt                                        >  1:377448/5‑123 (MQ=255)                                                                                       gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCt                                        <  2:377448/119‑1 (MQ=255)                                                                                       gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGa                        >  2:244566/5‑139 (MQ=255)                                                                                            aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAAt                                                                                                  <  1:456933/60‑1 (MQ=255)                                                                                            aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAAt                                                                                                  >  2:456933/1‑60 (MQ=255)                                                                                            aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                              >  2:855111/1‑64 (MQ=255)                                                                                            aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                              <  1:855111/64‑1 (MQ=255)                                                                                            aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGgtagt                                                                                       <  1:62225/71‑1 (MQ=255)                                                                                            aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGgtagt                                                                                       >  2:62225/1‑71 (MQ=255)                                                                                            aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTa                   >  2:474889/1‑139 (MQ=255)                                                                                            aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGAAGTACGTTGCGCCTTGCGCATCTGTTTTTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTa                   >  1:39430/1‑139 (MQ=255)                                                                                               gATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtga                <  1:827362/139‑1 (MQ=255)                                                                                                aTAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtgat               <  1:474889/139‑1 (MQ=255)                                                                                                  aaaTATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGc                          >  1:912022/1‑126 (MQ=255)                                                                                                  aaaTATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGc                          <  2:912022/126‑1 (MQ=255)                                                                                                  aaaTATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtgattg             >  1:521635/1‑139 (MQ=255)                                                                                                     tatCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGata          <  1:244566/139‑1 (MQ=255)                                                                                                        ccGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATaa       >  1:568815/1‑139 (MQ=255)                                                                                                          gcgcTGGCGTG‑ATATTGTCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt     <  2:476892/139‑1 (MQ=255)                                                                                                             cTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCAc  >  1:912817/1‑139 (MQ=255)                                                                                                                     |                                                                                                                                    GTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAATCAC  >  NC_000913/4424427‑4424674

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