Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,472,123	A→G	Y89H (TAC→CAC)	yfdM ←	CPS‑53 (KpLE1) prophage; putative methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,472,123	0	A	G	100.0%	143.2 / NA	41	Y89H (TAC→CAC)	yfdM	CPS‑53 (KpLE1) prophage; putative methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/23);  total (18/23)

AGCTGGTCATGCCGCCCCTTCTCCAGGCTGGCACCGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTACGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGCAGTGTAATACGCGGCACATTTGGCGTTATCACCATCAGTAAACAGATCCAGAACAA  >  NC_000913/2471989‑2472257                                                                                                                                       |                                                                                                                                       agcTGGTCATGCCGCCCCTTCTCCAGGCTGGCACCGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCgggg                                                                                                                                   <  2:302515/140‑1 (MQ=255)   cTGGTCATGCCGCCCCTTCTCCAGGCTGGCACCGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGa                                                                                                                                 <  1:472052/140‑1 (MQ=255)                cccTTCTCCAGGCTGGCACCGGAGGCCTCCACAACGCTGCAAAGCTCTTCCGGCAACATGTGCAAGTAGTCATCATCCAGATCGGTGGCTGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGgcag                                                                                                                    <  2:451635/140‑1 (MQ=17)                cccTTCTCCAGGCTGGCACCGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGgcag                                                                                                                    >  1:186429/1‑140 (MQ=255)                 ccTTCTCCAGGCTGGCACCGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGgcagc                                                                                                                   >  2:498041/1‑140 (MQ=255)                          ggCTGGCACCGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAg                                                                                                          <  2:45741/140‑1 (MQ=255)                           gCTGGCACCGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGc                                                                                                         >  1:140790/1‑140 (MQ=255)                           gCTGGCACCGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGc                                                                                                         <  2:12080/140‑1 (MQ=255)                               gCACCGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCg                                                                                                     >  1:122360/1‑140 (MQ=255)                                cACCGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAg                                                                                                          >  2:64675/1‑134 (MQ=255)                                cACCGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAg                                                                                                          <  1:64675/134‑1 (MQ=255)                                cACCGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGc                                                                                                    >  1:351481/1‑140 (MQ=255)                                  ccGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCaa                                                                                                  >  2:557697/1‑140 (MQ=255)                                   cGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAg                                                                                                 <  1:126779/140‑1 (MQ=255)                                   cGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAg                                                                                                 <  1:147574/140‑1 (MQ=255)                                   cGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAg                                                                                                 <  1:466836/140‑1 (MQ=255)                                          cTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCt                                                                                          >  2:86097/1‑140 (MQ=255)                                                  cGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTc                                                                                  >  2:299479/1‑140 (MQ=255)                                                    cTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCACAAGACGTTCTGACCAGTCAt                                                                                >  1:61183/1‑140 (MQ=255)                                                     tGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATg                                                                               >  1:564773/1‑140 (MQ=255)                                                               ctcctGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGcgcg                                                                     <  1:557697/140‑1 (MQ=255)                                                                 cctGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGtt                                                                   <  1:523907/140‑1 (MQ=255)                                                                 cctGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGtt                                                                   <  2:351481/140‑1 (MQ=255)                                                                         aTGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCg                                                           <  2:564773/140‑1 (MQ=255)                                                                         aTGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCg                                                           <  2:186429/140‑1 (MQ=255)                                                                         aTGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCACAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCg                                                           <  2:61183/140‑1 (MQ=255)                                                                          tGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGc                                                          <  2:71857/140‑1 (MQ=255)                                                                           gTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGcgcg                                                                     <  1:458298/128‑1 (MQ=255)                                                                           gTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGcgcg                                                                     >  2:458298/1‑128 (MQ=255)                                                                               aaGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGCAgtgt                                                     >  1:286545/1‑140 (MQ=255)                                                                                    gtCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGCAGTGTAATAc                                                <  1:299479/140‑1 (MQ=255)                                                                                     tcatcaTCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGCa                                                         >  2:78096/1‑130 (MQ=255)                                                                                     tcatcaTCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGCa                                                         <  1:78096/130‑1 (MQ=255)                                                                                     tcatcaTCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGCAGTGTAATAcg                                               >  2:22842/1‑140 (MQ=255)                                                                                                   ggTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGCAGTGTAATACGCGGCACATTTGGCg                                 <  2:122360/140‑1 (MQ=255)                                                                                                     tGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGCAGTGTAATACGCGGCACATTTGGCGtt                               >  2:404181/1‑140 (MQ=255)                                                                                                     tGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGCAGTGTAATACGCGGCACATTTGGCGtt                               <  1:498041/140‑1 (MQ=255)                                                                                                     tGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGCAGTGTAATACGCGGCACATTTGGCGtt                               >  1:587766/1‑140 (MQ=255)                                                                                                            cGTGCCGGTCTGGGGGTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGCAGTGTAATACGCGGCACATTTGGCGTTATCACCa                        <  1:86097/140‑1 (MQ=255)                                                                                                                        ggggTTACCACGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGCAGTGTAATACGCGGCACATTTGGCGTTATCACCATCAGTAAACAGa            >  2:478939/1‑140 (MQ=255)                                                                                                                                  cGGGTGCGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGCAGTGTAATACGCGGCACATTTGGCGTTATCACCATCAGTAAACAGATCCAGAACaa  <  2:587766/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       AGCTGGTCATGCCGCCCCTTCTCCAGGCTGGCACCGGAGGCCTCCACAACGCTGCAAAGCTCCTCCTGCAACATGTCAAAGTAGTCATCATCCAGATCGGTGGCAGGCGTGCCGGTCTGGGGGTTACCACGGGTACGGGGGATTACCAAAGGCAGCACCTTTAAGCTCCGCAAGACGTTCTGACCAGTCATGCGCCAGCGCGTTGTCTTCCGCAGTGTAATACGCGGCACATTTGGCGTTATCACCATCAGTAAACAGATCCAGAACAA  >  NC_000913/2471989‑2472257

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