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%	132.2 / NA	37	V945A (GTG→GCG)	stfR	Rac prophage, putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/20);  total (17/20)

ACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAA  >  NC_000913/1431748‑1432014                                                                                                                                       |                                                                                                                                     accaGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATc                                                                                                                                  <  2:597327/139‑1 (MQ=255)          cTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATg                                                                                                                         <  2:378181/139‑1 (MQ=255)                       gcACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGACAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGa                                                                                                            >  2:710095/1‑139 (MQ=255)                         aCCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGACAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGAtt                                                                                                          <  1:710095/139‑1 (MQ=255)                            aaTCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAg                                                                                                       >  2:880097/1‑139 (MQ=255)                                        aTCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGcc                                                                                           <  1:880097/139‑1 (MQ=255)                                                  gTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCggg                                                                                 <  2:934237/139‑1 (MQ=255)                                                   ttCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCggg                                                                                 <  1:310217/138‑1 (MQ=255)                                                   ttCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCggg                                                                                 >  2:310217/1‑138 (MQ=255)                                                         cTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTAtt                                                                          <  2:1190143/139‑1 (MQ=255)                                                                      gATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAg                                                             >  1:233339/1‑139 (MQ=255)                                                                      gATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAg                                                             >  2:233852/1‑139 (MQ=255)                                                                        tGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGa                                                           <  1:692145/139‑1 (MQ=255)                                                                        tGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGa                                                           <  2:766258/139‑1 (MQ=255)                                                                        tGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGa                                                           <  1:1025026/139‑1 (MQ=255)                                                                            ggggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGc                                                       >  1:388196/1‑139 (MQ=255)                                                                             gggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTgg                                                                                                                    <  2:889796/77‑1 (MQ=255)                                                                             gggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCa                                                      >  2:508873/1‑139 (MQ=255)                                                                             gggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCACTGATATGCGTgg                                                                                                                    >  1:889796/1‑77 (MQ=255)                                                                                 aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATcc                                                                                                                                 >  1:653643/1‑60 (MQ=255)                                                                                 aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATcc                                                                                                                                 <  2:653643/60‑1 (MQ=255)                                                                                 aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGt                                                                                     <  2:757083/104‑1 (MQ=255)                                                                                 aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGt                                                                                     >  1:757083/1‑104 (MQ=255)                                                                                             aaTCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCa                                            <  2:1075848/133‑1 (MQ=255)                                                                                             aaTCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCa                                            >  1:1075848/1‑133 (MQ=255)                                                                                                      aCCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcgccag                             >  1:202291/1‑139 (MQ=255)                                                                                                       cccGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcgccagc                            <  2:233339/139‑1 (MQ=255)                                                                                                             aCTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATcc                      >  2:1317165/1‑139 (MQ=255)                                                                                                                         ttATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGtctc                                                                     <  2:11614/80‑1 (MQ=255)                                                                                                                         ttATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGtctc                                                                     >  1:11614/1‑80 (MQ=255)                                                                                                                                tCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTa                                                   >  2:203481/1‑91 (MQ=255)                                                                                                                                tCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTa                                                   <  1:203481/91‑1 (MQ=255)                                                                                                                                 cAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCggg                                                                                 <  2:397216/60‑1 (MQ=255)                                                                                                                                 cAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCggg                                                                                 >  1:397216/1‑60 (MQ=255)                                                                                                                                 cAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGaa  <  2:388196/139‑1 (MQ=255)                                                                                                                                  aGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcg                                 >  1:85053/1‑107 (MQ=255)                                                                                                                                  aGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcg                                 <  2:85053/107‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     ACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAA  >  NC_000913/1431748‑1432014

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