Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,863,186	C→T	E86K (GAA→AAA)	SAUSA300_RS14645 ←	phage tail protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,863,186	0	C	T	100.0%	70.7 / NA	26	E86K (GAA→AAA)	SAUSA300_RS14645	phage tail protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (1/25);  total (1/25)

TTTTTTAATATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTCGTACCAGTTAGTAATATT‑AAAAAAACGACCAATTCTTGTTTTTTCTCTAAATTTAACAGGATATAAGTATTCAGACAAAAATGTACATTTCGTCTT  >  NC_007793/2863091‑2863282                                                                                                |                                                                                                  ttttttAATATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTag                                                                                               <  1:1509395/100‑2 (MQ=255)  tttttAATATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTAcc                                                                                              <  1:991573/100‑1 (MQ=255)  tttttAATATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTAcc                                                                                              <  1:5215817/100‑1 (MQ=255)  tttttAATATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTAcc                                                                                              <  1:5153361/100‑1 (MQ=255)  tttttAATATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTAcc                                                                                              <  1:2819404/100‑1 (MQ=255)  tttttAATATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTAcc                                                                                              <  1:2749352/100‑1 (MQ=255)  tttttAATATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTAcc                                                                                              <  1:2403025/100‑1 (MQ=255)    tttAGTATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAg                                                                                            <  1:4988175/100‑1 (MQ=255)    tttAATATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAg                                                                                            <  1:575662/100‑1 (MQ=255)      tAATATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAGtt                                                                                          <  1:4723081/100‑1 (MQ=255)      tAATATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAGtt                                                                                          <  1:1003768/100‑1 (MQ=255)      tAATATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAGtt                                                                                          <  1:827451/100‑1 (MQ=255)      tAATATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAGtc                                                                                          <  1:972372/100‑2 (MQ=255)      tAATACATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAGtt                                                                                          <  1:3137528/100‑1 (MQ=255)          atataAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAGTTAGTa                                                                                      <  1:3489577/100‑1 (MQ=255)                aGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAGTTAGTAATATT‑a                                                                               <  1:4563479/100‑1 (MQ=255)                                       tttttATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAGTTAGTAATATT‑AAAAAAACGACCAATTCTTGtggt                                                        <  1:3908280/100‑4 (MQ=255)                                        ctttATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAGTTAGTAATATT‑AAAAAAACGACCAATTCTTGttttt                                                       >  1:2556560/2‑100 (MQ=255)                                                           gCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAGTTAGTAATATT‑AAAAAAACGACCAATTCTTGTTTTTTCTCTAAATTTAACAGGat                                    <  1:4338585/100‑1 (MQ=255)                                                           gCAATAACATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAGTTAGTAATATT‑AAAAAAACGACCAATTCTTGTTTTTTCTCTAAATTTAACAGGat                                    <  1:4454925/100‑1 (MQ=255)                                                             aaTAACATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAGTTAGTAATATT‑AAAAAAACGACCAATTCTTGTTTTTTCTCTAAATTTAACAGGatat                                  <  1:2481065/100‑1 (MQ=255)                                                                 aCATACTTCTCATCGGTAGAATGGATTTCTTTGTACCAGTTAGTAATATTAAAAAAAACGACCAATTCTTGTTTTTTCTCTAAATTTAACAGGATATAAg                               <  1:2939999/100‑1 (MQ=255)                                                                     acctctcATCGGTAGAATGAATTTCTTTGTGCCAGTTAGTAATATT‑AAAAAAACGACCAATTCTTGTTTTTTCTCTAAATTTAACAGGATATAAGTATTc                          <  1:485850/97‑1 (MQ=255)                                                                     aCTTCTCATCGGTAGAATGGATTTCTTTGTACCAGTTAGTAATATT‑AAAAAAACGACCAATTCTTGTTTTTTCTCTAAATTTAACAGGATATAAGTATTc                          <  1:542089/100‑1 (MQ=255)                                                                         ctcATCGGTAGAATGGATTTCTTTGTACCAGTTAGTAATATT‑AAAAAAACGACCAATTCTTGTTTTTTCTCTAAATTTAACAGGATATAAGTATTCAgcg                      <  1:2580812/100‑3 (MQ=255)                                                                                   gAATGGATTTCTTTGTACCAGTTAGTAATATT‑AAAAAAACGACCAATTCTTGTTTTTTCTCTAAATTTAACAGGATATAAGTATTCAGACAAAAATGTAc            <  1:4133527/100‑1 (MQ=255)                                                                                             ctttGTACCAGTTAGTAATATT‑AAAAAAACGACCAATTCTTGTTTTTTCTCTAAATTTAACAGGATATAAGTATTCAGACAAAAATGTACATTTCGTCtt  <  1:5538952/100‑1 (MQ=255)                                                                                                |                                                                                                  TTTTTTAATATATAAAGAATAATCAATTAAGTCTTTGTTTTTTATATTTTTAATAGTTGCAATAACATACTTCTCATCGGTAGAATGGATTTCTTCGTACCAGTTAGTAATATT‑AAAAAAACGACCAATTCTTGTTTTTTCTCTAAATTTAACAGGATATAAGTATTCAGACAAAAATGTACATTTCGTCTT  >  NC_007793/2863091‑2863282

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