Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,208,534	Δ1 bp	coding (403/414 nt);coding (586/603 nt)	tfaP → tfaE ←	e14 prophage, uncharacterized protein e14 prophage, putative tail fiber assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,208,534	0	C	.	100.0%	147.6 / NA	37	coding (403/414 nt) coding (586/603 nt)	tfaP tfaE	e14 prophage, uncharacterized protein e14 prophage, putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (21/16);  total (21/16)

AGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCA  >  NC_000913/1208403‑1208665                                                                                                                                    |                                                                                                                                    aGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGa                                                                                                                             >  1:391536/1‑139 (MQ=255) aGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGa                                                                                                                             >  1:680549/1‑139 (MQ=255)             gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                 <  1:259/139‑1 (MQ=255)             gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                 <  2:546452/139‑1 (MQ=255)             gAAGACTTAAAGAAGATGAGCTAAAAAATTATAATTTGTGGCTGGATTATCTGGAAGCACTGAAACTGGTTGATACATCCAGTGCGCCAAATATTGAATGGCTTACCCCTCCGGCAGTT‑AGGCCAGATGACATCCgggg                                                                                                                 <  1:734207/139‑3 (MQ=18)                    tAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCt                                                                                                          >  2:234219/1‑139 (MQ=255)                      aagaagATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTg                                                                                                            >  2:631911/1‑135 (MQ=255)                      aagaagATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTg                                                                                                            <  1:631911/135‑1 (MQ=255)                           gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGTAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                 <  2:163176/125‑1 (MQ=255)                           gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGTAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                 >  1:163176/1‑125 (MQ=255)                           gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                                                                                                   >  1:482244/1‑139 (MQ=255)                                 cTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGcc                                                                                             >  1:767994/1‑139 (MQ=255)                                          tataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGcgcg                                                                                                               <  2:379595/112‑1 (MQ=255)                                          tataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGcgcg                                                                                                               >  1:379595/1‑112 (MQ=255)                                          tataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcgcg                                                                                    >  2:295867/1‑139 (MQ=255)                                           ataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                                   <  2:242359/139‑1 (MQ=255)                                           ataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                                   <  1:234219/139‑1 (MQ=255)                                                             aTCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAg                                                                 >  2:256002/1‑139 (MQ=255)                                                                   aaGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATc                                                                         >  2:17979/1‑125 (MQ=255)                                                                   aaGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATc                                                                         <  1:17979/125‑1 (MQ=255)                                                                         tGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                                     <  2:158182/139‑1 (MQ=255)                                                                         tGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                                     <  1:267948/139‑1 (MQ=255)                                                                            aaCTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGc                                                                     >  1:548835/1‑120 (MQ=255)                                                                            aaCTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGc                                                                     <  2:548835/120‑1 (MQ=255)                                                                            aaCTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGgttg                                                  >  2:152424/1‑139 (MQ=255)                                                                                ggTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGc                                                                     <  1:59046/116‑1 (MQ=255)                                                                                ggTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGc                                                                     >  2:59046/1‑116 (MQ=255)                                                                                           ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa                                   >  1:42937/1‑139 (MQ=255)                                                                                               tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGttt                               >  2:274280/1‑139 (MQ=255)                                                                                                   ccAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                 <  1:7808/53‑1 (MQ=14)                                                                                                   ccAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                 >  2:7808/1‑53 (MQ=14)                                                                                                   ccAGATATTGAATGACTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGt                           >  2:395716/1‑139 (MQ=255)                                                                                                               tGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTTACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAAt               <  2:482244/139‑1 (MQ=255)                                                                                                                  cttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCAACGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTc            >  2:475852/3‑139 (MQ=255)                                                                                                                          tCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTGTACGTCCGGCCTGCAATTTCAGTTGAAt    <  2:495157/139‑1 (MQ=255)                                                                                                                           ccGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATc   >  1:616445/1‑139 (MQ=255)                                                                                                                            cGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCa  >  2:164663/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                    AGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCA  >  NC_000913/1208403‑1208665

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: ‑