Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,428,030	C→G	A233P (GCC→CCC)	insH‑5 ←	Rac prophage; IS5 transposase and trans‑activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,428,030	0	C	G	100.0%	95.5 / NA	30	A233P (GCC→CCC)	insH‑5	Rac prophage; IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (12/18);  total (12/18)

TTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTAC  >  NC_000913/1427954‑1428107                                                                             |                                                                              ttCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGc                                                                           <  2:245474/81‑1 (MQ=255)      gcgcGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGctcctc                                                                      >  2:245518/1‑81 (MQ=255)      gcgcGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGctcctc                                                                      >  2:245516/1‑81 (MQ=255)      gcgcGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGctcctc                                                                      >  2:245507/1‑81 (MQ=255)      gcgcGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGctcct                                                                       <  1:245475/80‑1 (MQ=255)        gcgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTcgc                                                                    <  1:245476/81‑1 (MQ=255)        gcgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTcgc                                                                    >  2:245510/1‑81 (MQ=255)         cgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTcgcg                                                                   >  2:245520/1‑81 (MQ=255)         cgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTcgcg                                                                   >  1:245508/1‑81 (MQ=255)          ggATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCgcgc                                                                  >  2:245511/1‑81 (MQ=255)            aTTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCtg                                                                <  1:245477/81‑1 (MQ=255)             ttCTGCTTCAGGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCtgt                                                               <  1:245478/81‑1 (MQ=255)                  cTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGgcgc                                                          >  2:245522/1‑81 (MQ=255)                      aaGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTcc                                                       <  1:245479/80‑1 (MQ=255)                          tttttACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGt                                                  <  2:245480/81‑1 (MQ=255)                            tttACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTa                                                 <  2:245481/80‑1 (MQ=255)                                ccTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCg                                             <  2:245482/80‑1 (MQ=255)                                    gCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATc                                        >  2:245521/1‑81 (MQ=255)                                      cGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCgg                                      <  2:245483/81‑1 (MQ=255)                                        ggACGCTAGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCt                                    <  2:245484/81‑1 (MQ=255)                                                 gCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAAt                           >  1:245519/1‑81 (MQ=255)                                                  cGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAAtt                          <  1:245485/81‑1 (MQ=18)                                                    aTCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGc                        <  1:245486/81‑1 (MQ=18)                                                     tCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGct                       <  2:245487/81‑1 (MQ=18)                                                     tCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGct                       <  1:245488/81‑1 (MQ=18)                                                       aGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGctc                      <  1:245489/80‑1 (MQ=21)                                                              ccACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCAt              >  2:245517/1‑81 (MQ=18)                                                                 caTCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGaa           <  2:245490/81‑1 (MQ=37)                                                                       ccTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGAt     <  1:245491/81‑1 (MQ=255)                                                                          cGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTAc  >  1:245528/1‑81 (MQ=255)                                                                             |                                                                              TTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTAC  >  NC_000913/1427954‑1428107

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑