Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,431,882	T→C	V945A (GTG→GCG)	stfR →	Rac prophage, putative tail fiber protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,431,882	0	T	C	100.0%	77.5 / NA	23	V945A (GTG→GCG)	stfR	Rac prophage, putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/10);  total (13/10)

CCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCAT  >  NC_000913/1431760‑1431973                                                                                                                           |                                                                                            cccAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGCCAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGt                                                                             >  2:680467/1‑139 (MQ=255)        ggCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGa                                                                      <  1:970922/139‑1 (MQ=255)             aCCAATCCCGTGGCCATCAGGTACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGAtt                                                                 <  1:241762/139‑1 (MQ=255)                              cAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAg                                                >  1:797901/1‑139 (MQ=255)                                   aCCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATc                                                                                         <  1:797994/93‑1 (MQ=255)                                   aCCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATc                                                                                         >  2:797994/1‑93 (MQ=255)                                       ttCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGc                                                           >  1:736499/1‑119 (MQ=255)                                       ttCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGc                                                           <  2:736499/119‑1 (MQ=255)                                        tCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGGGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCGAGTGGGCGGGcc                                      >  1:326388/1‑139 (MQ=255)                                           gTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAgg                                                             >  1:559953/1‑113 (MQ=255)                                           gTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAgg                                                             <  2:559953/113‑1 (MQ=255)                                                    tGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAg                          >  1:281298/1‑139 (MQ=255)                                                           aTGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGtctc                            >  1:45477/1‑130 (MQ=255)                                                           aTGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGtctc                            <  2:45477/130‑1 (MQ=255)                                                                 gggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCa             >  1:371738/1‑139 (MQ=255)                                                                                 aaTCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGc    <  2:118069/132‑1 (MQ=255)                                                                                 aaTCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGc    >  1:118069/1‑132 (MQ=255)                                                                                            ccGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATg                                                                                <  1:873148/45‑1 (MQ=255)                                                                                            ccGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATg                                                                                >  2:873148/1‑45 (MQ=255)                                                                                                    tGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGTCGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCAt  <  1:667174/115‑1 (MQ=255)                                                                                                    tGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCAt  >  2:667174/1‑115 (MQ=255)                                                                                                                   gTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCAt  <  2:543076/100‑1 (MQ=255)                                                                                                                   gTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCAt  >  1:543076/1‑100 (MQ=255)                                                                                                                           |                                                                                            CCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCAT  >  NC_000913/1431760‑1431973

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