Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,419,159	G→A	*163* (TAG→TAA)  S51S (AGC→AGT)	sieB → ydaF ←	phage superinfection exclusion protein, Rac prophage uncharacterized protein, Rac prophage

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,419,159	0	G	A	100.0%	133.6 / NA	43	*163* (TAG→TAA)  S51S (AGC→AGT)	sieB ydaF	phage superinfection exclusion protein, Rac prophage uncharacterized protein, Rac prophage
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/19);  total (24/19)

ATTAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAGCTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGT  >  NC_000913/1419030‑1419289                                                                                                                                  |                                                                                                                                   aTTAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACtgtg                                                                                                                           >  2:1279667/1‑139 (MQ=255)   tAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACtgtgtg                                                                                                                         <  1:450611/139‑1 (MQ=255)                             cAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAg                                                                                               >  1:1051046/1‑139 (MQ=255)                                agTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGttt                                                                                            >  1:639742/1‑139 (MQ=255)                                agTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGttt                                                                                            >  2:34235/1‑139 (MQ=255)                                agTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGttt                                                                                            >  1:18929/1‑139 (MQ=255)                                          cAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTaaa                                                                                  >  2:417215/1‑139 (MQ=255)                                                   cTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCt                                                                         <  1:1279667/139‑1 (MQ=255)                                                     tattatCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                                       <  2:1647013/139‑1 (MQ=255)                                                     tattatCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                                       <  2:639742/139‑1 (MQ=255)                                                     tattatCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                                       <  2:1294863/139‑1 (MQ=255)                                                     tattatCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                                       <  1:466385/139‑1 (MQ=255)                                                          tctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                                  <  2:1051046/139‑1 (MQ=255)                                                          tctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                                  <  1:417215/139‑1 (MQ=255)                                                           ctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCaa                                                                                                                <  1:1262554/92‑1 (MQ=255)                                                           ctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCaa                                                                                                                >  2:1262554/1‑92 (MQ=255)                                                           ctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCt                                                                         <  2:1024668/131‑1 (MQ=255)                                                           ctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCt                                                                         >  1:1024668/1‑131 (MQ=255)                                                            tcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGtata                                                                <  2:18929/139‑1 (MQ=255)                                                             cGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGt                                                                               <  1:109517/123‑1 (MQ=255)                                                             cGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGt                                                                               >  2:109517/1‑123 (MQ=255)                                                               tCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGt                                                             >  1:357165/1‑139 (MQ=255)                                                                               ccATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTACATGTTTTGCTTTTGGTATACGTCATGTCAGTGAACGtt                                             >  2:396958/1‑139 (MQ=255)                                                                               ccATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGtt                                             >  1:867400/1‑139 (MQ=255)                                                                                     tATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATcc                                       <  1:34235/139‑1 (MQ=255)                                                                                                    gAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACAcg                        >  1:14526/1‑139 (MQ=255)                                                                                                        agcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                                       <  2:468412/56‑1 (MQ=255)                                                                                                        agcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                                       >  1:468412/1‑56 (MQ=255)                                                                                                         gcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                                       <  2:1053545/55‑1 (MQ=255)                                                                                                         gcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGca                   <  2:357165/139‑1 (MQ=255)                                                                                                         gcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGca                   >  2:314810/1‑139 (MQ=255)                                                                                                         gcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTACTTCTCCAACAATCGgcg                                                                                                       >  1:1053545/1‑55 (MQ=255)                                                                                                          cagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                                       >  2:158338/1‑86 (MQ=255)                                                                                                          cagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                                       <  1:158338/86‑1 (MQ=255)                                                                                                          cagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGcac                  >  2:261490/1‑139 (MQ=255)                                                                                                                  gACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGa          >  1:103544/1‑139 (MQ=255)                                                                                                                    cGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                                       >  2:208638/1‑76 (MQ=255)                                                                                                                    cGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                                       <  1:208638/76‑1 (MQ=255)                                                                                                                       tttAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCa                                                                                                   >  2:1643940/1‑45 (MQ=255)                                                                                                                       tttAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCa                                                                                                   <  1:1643940/45‑1 (MQ=255)                                                                                                                        ttAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTT‑GGGAACACGCCGCACACCAGGGATTCgttg   >  2:871994/1‑139 (MQ=255)                                                                                                                          aaTCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCgttgt  >  1:889352/1‑139 (MQ=255)                                                                                                                          aaTCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCgttgt  >  2:644811/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                   ATTAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAGCTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGT  >  NC_000913/1419030‑1419289

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