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%	212.2 / NA	59	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 . (29/30);  total (29/30)

GTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATC  >  NC_000913/1208400‑1208664                                                                                                                                       |                                                                                                                                   gTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGcc                                                                                                                               <  1:41669/139‑1 (MQ=255)    aGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGa                                                                                                                            >  1:833900/1‑139 (MQ=255)          cATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACAt                                                                                                                      >  2:1699806/1‑139 (MQ=255)          cATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACAt                                                                                                                      >  1:1555979/1‑139 (MQ=255)          cATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACAt                                                                                                                      >  1:1621308/1‑139 (MQ=255)          cATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACAt                                                                                                                      >  2:698346/1‑139 (MQ=255)                gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                <  2:701986/139‑1 (MQ=255)                  aGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGcgcg                                                                                                              <  1:235262/139‑1 (MQ=255)                     cTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTg                                                                                                           <  2:833900/139‑1 (MQ=255)                     cTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTg                                                                                                           <  1:698346/139‑1 (MQ=255)                        aaaGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTg                                                                                                        <  2:1555793/139‑1 (MQ=255)                         aagaagATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTgg                                                                                                       >  2:1802878/1‑139 (MQ=255)                           gaagaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTa                                                                                                     <  2:1621308/139‑1 (MQ=255)                           gaagaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTa                                                                                                     <  2:1555979/139‑1 (MQ=255)                              gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                                                                                                  >  1:443776/1‑139 (MQ=255)                               aTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTg                                                                                                 <  1:1802878/139‑1 (MQ=255)                                                   ttGTGGCTGGATTATCTAGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGAt                                                                                                                           >  1:854684/1‑93 (MQ=17)                                                   ttGTGGCTGGATTATCTAGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGAt                                                                                                                           <  2:854684/93‑1 (MQ=17)                                                       ggCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAAt                                                                         >  1:1796836/1‑139 (MQ=255)                                                            gATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                >  2:239201/1‑95 (MQ=255)                                                            gATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                <  1:239201/95‑1 (MQ=255)                                                                            tGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTa                                                                           >  1:1468339/1‑116 (MQ=255)                                                                            tGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTa                                                                           <  2:1468339/116‑1 (MQ=255)                                                                                    gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCa                                                                                        <  1:566850/95‑1 (MQ=255)                                                                                    gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCa                                                                                        >  2:566850/1‑95 (MQ=255)                                                                                    gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                             <  1:698664/106‑1 (MQ=255)                                                                                    gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                             >  2:698664/1‑106 (MQ=255)                                                                                    gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTct                                            <  2:1107203/139‑1 (MQ=255)                                                                                        aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGca                                                                   <  2:1684672/112‑1 (MQ=255)                                                                                        aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGca                                                                   >  1:1684672/1‑112 (MQ=255)                                                                                         tACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                <  1:508345/66‑1 (MQ=14)                                                                                         tACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                >  2:508345/1‑66 (MQ=14)                                                                                            aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                    <  1:1699806/139‑1 (MQ=255)                                                                                              ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATc                                                                                                                     <  2:148394/56‑1 (MQ=14)                                                                                              ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATc                                                                                                                     >  1:148394/1‑56 (MQ=14)                                                                                              ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                                                                                                  <  2:687691/75‑1 (MQ=11)                                                                                              ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                                                                                                  >  1:687691/1‑75 (MQ=11)                                                                                              ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa                                  >  2:80647/1‑139 (MQ=255)                                                                                              ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa                                  >  1:310455/1‑139 (MQ=255)                                                                                                            atTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCgg                      >  2:504187/1‑137 (MQ=255)                                                                                                            atTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCgg                      <  1:504187/137‑1 (MQ=255)                                                                                                              tGATTGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGGGTTAAGTCTGGTTGTCTCTGCCTGCGTCAGTTTACGTCCGGCCTg                  >  1:1181477/1‑139 (MQ=18)                                                                                                                aaTGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgc                                          >  2:241760/1‑113 (MQ=255)                                                                                                                aaTGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgc                                          <  1:241760/113‑1 (MQ=255)                                                                                                                   ggCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAAtt             <  1:80647/139‑1 (MQ=255)                                                                                                                    gcttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                             <  2:43952/71‑1 (MQ=11)                                                                                                                    gcttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                             >  1:43952/4‑74 (MQ=11)                                                                                                                     cttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGcgcg                                                                                                              <  2:740261/38‑1 (MQ=14)                                                                                                                     cttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGcgcg                                                                                                              >  1:740261/3‑40 (MQ=14)                                                                                                                     cttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGt                                                                                                      >  1:1005287/3‑48 (MQ=11)                                                                                                                     cttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGt                                                                                                      <  2:1005287/46‑1 (MQ=11)                                                                                                                     cttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTc           >  2:778157/3‑139 (MQ=255)                                                                                                                     cttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTc           >  2:793611/3‑139 (MQ=255)                                                                                                                        aCGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGt        <  1:793611/139‑1 (MQ=255)                                                                                                                           ccTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGa     <  2:1796836/139‑1 (MQ=255)                                                                                                                             tCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTACAATTTCAGTTGAAt   <  2:310455/139‑1 (MQ=255)                                                                                                                              ccGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTCACGTCCGGCCTGCAATTTCAGTTGAATc  >  2:1100267/1‑139 (MQ=255)                                                                                                                               cGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                                         >  2:1308779/1‑99 (MQ=255)                                                                                                                               cGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                                         <  1:1308779/99‑1 (MQ=255)                                                                                                                                       |                                                                                                                                   GTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATC  >  NC_000913/1208400‑1208664

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