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%	167.0 / NA	53	*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 (28/25);  total (28/25)

AAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAGCTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACAC  >  NC_000913/1419040‑1419275                                                                                                                        |                                                                                                                     aaGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCa                                                                                                   >  2:1533756/1‑139 (MQ=255)   gggCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAAc                                                                                                 >  2:612019/1‑139 (MQ=255)              aaCATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcgcgc                                                                                      >  1:571289/1‑139 (MQ=255)                cATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCAc                                                                                    <  1:120671/139‑1 (MQ=255)                 aTCAGAGTGATAGCGCAAACTGTTCTTATTATCACGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCAcc                                                                                   <  1:592665/139‑1 (MQ=255)                     gagTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGtt                                                                               >  1:1505321/1‑139 (MQ=255)                           tAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGttt                                                                              >  1:1224772/1‑134 (MQ=255)                           tAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGttt                                                                              <  2:1224772/134‑1 (MQ=255)                              cgcAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTa                                                                      >  2:307453/1‑139 (MQ=255)                              cgcAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTa                                                                      >  2:133310/1‑139 (MQ=255)                              cgcAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTa                                                                      >  1:814911/1‑139 (MQ=255)                                cAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                         <  2:530713/118‑1 (MQ=255)                                cAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                         >  1:530713/1‑118 (MQ=255)                                cAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCa                                                                                     >  2:400132/1‑122 (MQ=255)                                cAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCa                                                                                     <  1:400132/122‑1 (MQ=255)                                  aaCTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATg                                                                  >  2:163216/1‑139 (MQ=255)                                    cTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGtt                                                                <  1:133310/139‑1 (MQ=255)                                       ttCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTg                                                             >  1:404165/1‑139 (MQ=255)                                        tCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGc                                                            <  2:1505321/139‑1 (MQ=255)                                           tattatCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                         <  1:307453/139‑1 (MQ=255)                                               atCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGt                                                     >  1:855141/1‑139 (MQ=255)                                                tctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                    <  1:163216/139‑1 (MQ=255)                                                tctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                    <  2:571289/139‑1 (MQ=255)                                                 ctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGttt                                                                              >  1:412272/1‑112 (MQ=255)                                                 ctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGttt                                                                              <  2:412272/112‑1 (MQ=255)                                                 ctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                    >  2:741644/1‑138 (MQ=255)                                                 ctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                    <  1:741644/138‑1 (MQ=255)                                                    gTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACg                                                >  1:935056/1‑139 (MQ=255)                                                      cACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                    >  1:1135354/1‑133 (MQ=255)                                                      cACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                    <  2:1135354/133‑1 (MQ=255)                                                      cACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTc                                              >  1:1492092/1‑139 (MQ=255)                                                           aaaaaTACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGc                                         <  1:612019/139‑1 (MQ=255)                                                            aaaaTACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCa                                        <  2:814911/139‑1 (MQ=255)                                                                   tCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACg                                 >  2:1103055/1‑139 (MQ=255)                                                                       aTTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCg                                                                                            <  1:1530111/76‑1 (MQ=255)                                                                       aTTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCg                                                                                            >  2:1530111/1‑76 (MQ=255)                                                                              gAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTgg                      <  1:1533756/139‑1 (MQ=255)                                                                                   tATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTggg                 >  2:1454156/1‑139 (MQ=255)                                                                                          gAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                         <  1:1336749/60‑1 (MQ=255)                                                                                          gAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                         >  2:1336749/1‑60 (MQ=255)                                                                                          gAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGt                                      >  2:839158/1‑111 (MQ=255)                                                                                          gAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGt                                      <  1:839158/111‑1 (MQ=255)                                                                                             cagcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACAcgccg       <  1:1454156/139‑1 (MQ=255)                                                                                             cagcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACAcgccg       <  2:855141/139‑1 (MQ=255)                                                                                               gcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                         >  2:1299608/1‑55 (MQ=255)                                                                                               gcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                         <  1:1299608/55‑1 (MQ=255)                                                                                               gcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCtt                                                          <  1:717808/86‑1 (MQ=255)                                                                                               gcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCtt                                                          >  2:717808/1‑86 (MQ=255)                                                                                               gcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGca     >  1:184098/1‑139 (MQ=255)                                                                                                cagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACAcgccg       <  2:338049/136‑1 (MQ=255)                                                                                                cagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACAcgccg       >  1:338049/1‑136 (MQ=255)                                                                                                             tttAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCacac  <  1:743060/128‑1 (MQ=255)                                                                                                             tttAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCacac  >  2:743060/1‑128 (MQ=255)                                                                                                                        |                                                                                                                     AAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAGCTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACAC  >  NC_000913/1419040‑1419275

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