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

GTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAT  >  NC_000913/4424427‑4424671                                                                                                                     |                                                                                                                                  gTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTc                                                                                                            >  1:1342399/1‑139 (MQ=255) gTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTc                                                                                                            >  1:257753/1‑139 (MQ=255)            cAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGaa                                                                                                 >  1:104287/1‑139 (MQ=255)            cAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGaa                                                                                                 >  2:338116/1‑139 (MQ=255)               ccGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGTGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTa                                                                                              <  2:1324134/139‑1 (MQ=255)                 gTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                            <  2:104287/139‑1 (MQ=255)                        aaaaaaaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑Ggtagt                                                                                    <  2:1420933/139‑1 (MQ=255)                                gCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTgcg                                                                            <  2:209428/139‑1 (MQ=255)                                       gCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTgcg                                                                     <  2:1381514/139‑1 (MQ=255)                                               ggCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGtt                                                             >  2:1511780/1‑139 (MQ=255)                                                   aagcaaaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTttctt                                                         <  2:1342399/139‑1 (MQ=255)                                                              ggCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGcc                                              >  1:368342/1‑139 (MQ=255)                                                              ggCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGcc                                              >  1:831287/1‑139 (MQ=255)                                                                  aTGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGAc                                          >  1:217120/1‑139 (MQ=255)                                                                   tGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGaa                                                                                                 <  1:1104050/84‑1 (MQ=255)                                                                   tGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGaa                                                                                                 >  2:1104050/1‑84 (MQ=255)                                                                   tGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACg                                         >  1:684691/1‑139 (MQ=255)                                                                     tttCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTa                                                                                              >  2:94887/1‑85 (MQ=255)                                                                     tttCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTa                                                                                              <  1:94887/85‑1 (MQ=255)                                                                            ggTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCt                                >  2:636243/1‑139 (MQ=255)                                                                                     tGGTGTTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGc                       <  2:257753/139‑1 (MQ=255)                                                                                       gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTccc                                                                                                                  >  1:1652766/5‑47 (MQ=37)                                                                                       gtgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCg                      >  2:1598884/5‑139 (MQ=255)                                                                                       gtgtTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTccc                                                                                                                  <  2:1652766/47‑1 (MQ=38)                                                                                       gtgtTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTttctt                                                         >  1:1376585/1‑103 (MQ=255)                                                                                       gtgtTAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTttctt                                                         <  2:1376585/103‑1 (MQ=255)                                                                                        tgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTcc                                                                                                           <  2:480377/50‑1 (MQ=38)                                                                                        tgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTcc                                                                                                           >  1:480377/4‑53 (MQ=38)                                                                                        tgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTgcg                                                                     <  2:154024/87‑1 (MQ=255)                                                                                        tgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTgcg                                                                     >  1:154024/4‑90 (MQ=255)                                                                                        tgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCAttt                                                    <  2:1218884/104‑1 (MQ=255)                                                                                        tgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCAttt                                                    >  1:1218884/4‑107 (MQ=255)                                                                                        tgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGc                       <  2:952534/133‑1 (MQ=255)                                                                                        tgttAGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGc                       >  1:952534/4‑136 (MQ=255)                                                                                            aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAAtt                                                                                               >  1:390442/1‑61 (MQ=255)                                                                                            aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAAtt                                                                                               <  2:390442/61‑1 (MQ=255)                                                                                            aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTa                >  2:909995/1‑139 (MQ=255)                                                                                            aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTa                >  1:381267/1‑139 (MQ=255)                                                                                            aGCGATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTa                >  2:321584/1‑139 (MQ=255)                                                                                               gATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtga             <  2:217120/139‑1 (MQ=255)                                                                                               gATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtga             <  2:81652/139‑1 (MQ=255)                                                                                               gATAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtga             <  1:1511780/139‑1 (MQ=255)                                                                                                 tAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtgat            >  2:1316113/1‑138 (MQ=255)                                                                                                 tAAATATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtgat            <  1:1316113/138‑1 (MQ=255)                                                                                                  aaaTATCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtgattg          >  2:116311/1‑139 (MQ=255)                                                                                                       tCCGCGCTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGAtata     >  1:1225296/1‑139 (MQ=255)                                                                                                          gcgcTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAt  <  1:338116/139‑1 (MQ=255)                                                                                                           cgcTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                            >  1:1623095/1‑49 (MQ=39)                                                                                                           cgcTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                            <  2:1623095/49‑1 (MQ=39)                                                                                                           cgcTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCAt                                                                  <  2:770674/74‑1 (MQ=255)                                                                                                           cgcTGGCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCAt                                                                  >  1:770674/1‑74 (MQ=255)                                                                                                               ggCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑Ggtag                                                                                     <  2:84651/51‑1 (MQ=255)                                                                                                               ggCGTG‑ATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑Ggtag                                                                                     >  1:84651/1‑51 (MQ=255)                                                                                                                     |                                                                                                                                  GTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGT‑GGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGATTGATATAAAT  >  NC_000913/4424427‑4424671

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