Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,635,293	G→A	T185T (ACC→ACT)	ydfN ←	Qin prophage; putative side tail fiber assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,635,293	0	G	A	100.0%	33.9 / NA	14	T185T (ACC→ACT)	ydfN	Qin prophage; putative side tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/6);  total (8/6)

ATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATG  >  NC_000913/1635164‑1635374                                                                                                                                  |                                                                                  ataACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGTGCCCCGGTGTTATTCGTGGATTTAGtg                                                                                <  2:537508/133‑2 (MQ=255)          aGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGTGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCa                                                                       >  1:68425/1‑133 (MQ=255)          aGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTCTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGTGCCCCGGTGTTATACGTGGAGTTAGTGACGTAATCa                                                                       >  1:175179/1‑133 (MQ=255)            ggAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGTGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCaaa                                                                     >  2:64050/1‑133 (MQ=255)                aaTGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGTGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAacgac                                                                 >  2:619440/1‑133 (MQ=255)                 aTGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGTGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAacgacg                                                                <  2:534236/133‑1 (MQ=255)                        tgtTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGTGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGt                                                         >  2:508828/1‑133 (MQ=255)                                                     ttGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGTGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATgcgc                            >  2:106159/1‑133 (MQ=17)                                                      tgtgTTGGTGCTCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGTGCCTCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCTAAATCTGTACTGGATgcgct                           <  2:104447/133‑1 (MQ=17)                                                       gtgtTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGTGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATgcgctg                          >  1:89694/1‑133 (MQ=17)                                                          ttGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGTGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATgcgctggcg                       <  1:38077/133‑1 (MQ=17)                                                              tGCGCACCCGCACTATTTCCGGACCCGCTATTACTATGGGTATGTGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCtgt                   <  1:195910/133‑1 (MQ=11)                                                                     cGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGTGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGtgtg            <  1:619440/133‑1 (MQ=14)                                                                               tGCGGTCCCGCTAATACTATGGGTATGTGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATg  >  2:606828/1‑133 (MQ=11)                                                                                                                                  |                                                                                  ATAACCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATG  >  NC_000913/1635164‑1635374

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