Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,791,609	Δ1 bp	coding (412/714 nt)	ydiV ←	anti‑FlhDC factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,791,606	0	G	.	100.0%	62.5 / NA	16	coding (415/714 nt)	ydiV	anti‑FlhDC factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (7/9);  total (7/9)

CAGCATGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGATACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCT  >  NC_000913/1791529‑1791674                                                                              |                                                                     cAGCATGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAggga                                                                   <  1:323580/81‑2 (MQ=255) cAGCATGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAggga                                                                   <  2:323579/81‑2 (MQ=255)       gACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAA‑GGGATACACt                                                            >  1:323605/1‑81 (MQ=255)         cgcgGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAA‑GGGATACACTTg                                                          <  1:323581/81‑1 (MQ=255)          gcgGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAA‑GGGATACACTTGc                                                         <  2:323582/81‑1 (MQ=255)              gTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAA‑GGGATACACTTGCGATa                                                     <  1:323583/81‑1 (MQ=255)                 aaaaGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAA‑GGGATACACTTGCGATAAg                                                   <  2:323584/80‑1 (MQ=255)                   aaGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAA‑GGGATACACTTGCGATAAGGa                                                 >  1:323613/1‑80 (MQ=255)                   aaGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAA‑GGGATACACTTGCGATAAGGAt                                                <  2:323585/81‑1 (MQ=255)                    aGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAA‑GGGATACACTTGCGATAAGGAt                                                >  1:323604/1‑80 (MQ=255)                     gACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAA‑GGGATACACTTGCGATAAGGATaa                                              <  1:323586/81‑1 (MQ=255)                        cATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAA‑GGGATACACTTGCGATAAGGATAAAAg                                           <  2:323587/81‑1 (MQ=255)                                       tCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAA‑GGGATACACTTGCGATAAGGATAAAAGACCGCGATTATCttt                            >  2:323615/1‑81 (MQ=255)                                                      cTGCTCCCAAATTACCTAACACAA‑GGGATACACTTTCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATg             >  2:323617/1‑81 (MQ=255)                                                         cTCCCAAATTACCTAACACAA‑GGGATACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGc            <  2:323588/79‑1 (MQ=255)                                                              aaaTTACCTAACACAA‑GGGATACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAAt     >  1:323609/1‑81 (MQ=255)                                                                 ttACCTAACACAA‑GGGATACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATtct  >  1:323612/1‑81 (MQ=255)                                                                 ttACCTAACACAA‑GGGATACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATtct  <  2:323589/81‑1 (MQ=255)                                                                              |                                                                     CAGCATGACGCGGGTGAAAAGACCATCAAAAACCGCTTTCATTGTGCTATTGCCTGCTCCCAAATTACCTAACACAAGGGGATACACTTGCGATAAGGATAAAAGACCGCGATTATCTTTTCCTTCATTGAGATGCGGGTAATTCT  >  NC_000913/1791529‑1791674

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