Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	573,777	(G)5→4	coding (48/141 nt)	ylcG →	DLP12 prophage; uncharacterized protein YlcG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	573,773	0	G	.	100.0%	76.0 / NA	19	coding (44/141 nt)	ylcG	DLP12 prophage; uncharacterized protein YlcG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (7/12);  total (7/12)

GCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGACCGATTACCGAATACTGAAGAATTACGCCAAAACCCTGACAGGAGCAGGAGTATGAAGTCAGAGA  >  NC_000913/573644‑573879                                                                                                                                  |                                                                                                           gCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTgggg                                                                                                         >  1:590744/1‑133 (MQ=255)       gTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTg                                                                                                  >  1:444778/1‑133 (MQ=255)                 aCCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATc                                                                                        >  2:397285/1‑133 (MQ=255)                        tgtgAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCACCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTccc                                                                                 >  2:63573/1‑133 (MQ=255)                            aaGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGt                                                                             <  1:391065/133‑1 (MQ=255)                             aGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGtt                                                                            <  2:590744/133‑1 (MQ=255)                             aGATGCCTAGTACCACAGATGGCAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGtt                                                                            <  1:286370/133‑1 (MQ=255)                                        aCCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCACCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGAcc                                                                 <  1:63573/133‑1 (MQ=255)                                          cAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGACCGa                                                               <  1:397285/133‑1 (MQ=255)                                          cAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTCAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGACCGa                                                               <  2:74988/133‑1 (MQ=255)                                                                              ggggAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGACCGATTACCGAATACTGAAGAATTACGCCAAAACCCTGAc                           <  2:69310/133‑1 (MQ=255)                                                                              ggggAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGACCGATTACCGAATACTGAAGAATTACGCCAAAACCCTGAc                           <  1:74435/133‑1 (MQ=255)                                                                                   atgaatgaTAGTTGAGTTTAATATAGCAGAACTTCTTCGCCGCCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGACCGATTACCGAATACTGAAGAATTACGCCAAAACCCTGAcaggag                      <  2:286604/133‑1 (MQ=255)                                                                                           tgTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGACCGATTACCGAATACTGAAGAATTACGCCAAAACCCTGAcaggagca                    <  2:469/127‑1 (MQ=255)                                                                                           tgTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGACCGATTACCGAATACTGAAGAATTACGCCAAAACCCTGAcaggagca                    >  1:469/1‑127 (MQ=255)                                                                                             tttGAGTTTAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGACCGATTACCGAATACTGAAGAATTACGCCAAAACCCTGACAGGAGCAGGAGTATg            <  2:76089/133‑1 (MQ=255)                                                                                                gAGTTTAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGACCGATTACCGAATACTGAAGAATTACGCCAAAACCCTGAcaggagca                    <  2:77516/122‑1 (MQ=255)                                                                                                gAGTTTAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGACCGATTACCGAATACTGAAGAATTACGCCAAAACCCTGAcaggagca                    >  1:77516/1‑122 (MQ=255)                                                                                                gAGTTTAATATGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGACCGATTACCGAATACTGAAGAATTACGCCAAAACCCTGACAGGAGCAGGAGTATGAAg         >  2:21154/1‑133 (MQ=255)                                                                                                       atatGGCAGAACTTCTTCGCCACCGCT‑GGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGACCGATTACCGAATACTGAAGAATTACGCCAAAACTCTGACAGTAGCGGGTGTATGAAGTCAgaga  >  1:183331/1‑133 (MQ=255)                                                                                                                                  |                                                                                                           GCTCAGGTCGTTGATTACCGCGTTGTGAAGATGCCTGTTACCAAAGGTGGGAGGCTGGAACTGACCATCACCGAAATGGGGAATGAATGATGTTTGAGTTTAATATGGCAGAACTTCTTCGCCACCGCTGGGGGCGTCTGCGCTTATATCGTTTCCCCGGTTCTGTTTTGACCGATTACCGAATACTGAAGAATTACGCCAAAACCCTGACAGGAGCAGGAGTATGAAGTCAGAGA  >  NC_000913/573644‑573879

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