Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,651,972	(T)7→8	coding (316/4962 nt)	yfhM ←	alpha2‑macroglobulin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,651,965	1	.	T	100.0%	61.0 / NA	18	S108N (AGC→AAC)	yfhM	alpha2‑macroglobulin
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (8/10);  total (8/10)

GTTTCGGTTCGAGGTGGCGTAAACGTAACTCTTTAAGATTATCTGACAGCTCCCAGGCACCATCCACTTTGCCGC‑TTTTTTTATCGACGACATGAATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCAC  >  NC_000913/2651891‑2652038                                                                            |                                                                          gTTTCGGTTCGAGGTGGCGTAAACGTAACTCTTTAAGATTATCTGACAGCTCCCAGGCACCATCCACTTTGCCGC‑tttttt                                                                     >  2:474632/1‑81 (MQ=255)      ggTTCGAGGTGGCGTAAACGTAACTCTTTAAGATTATCTGACAGCTCCCAGGCACCATCCACTTTGCCGC‑ttttttttat                                                                 <  1:474601/80‑4 (MQ=255)        ttCGAGGTGGCGTAAACGTAACTCTTTAAGATTATCTGACAGCTCCCAGGCACCATCCACTTTGCCGCTTTTTTTTATcga                                                               >  2:474637/1‑81 (MQ=255)          cGAGGTGGCGTAAACGTAACTCTTTAAGATTATCTGACAGCTCCCAGGCACCATCCACTTTGCCGCTTTTTTTTATcgac                                                              >  2:474626/1‑80 (MQ=255)            aGGTGGCGTAAACGTAACTCTTTAAGATTATCTGACAGCTCCCAGGCACCATCCACTTTGCCGCTTTTTTTTATCgacgac                                                           <  1:474603/81‑1 (MQ=255)                  cGTAAACGTAACTCTTTAAGATTATCTGACAGCTCCCAGGCACCATCCACTTTGCCGCTTTTTTTTATCGACGACATGaa                                                      >  2:474640/1‑80 (MQ=255)                  cGTAAACGTAACTCTTTAAGATTATCTGACAGCTCCCAGGCACCATCCACTTTGCCGCTTTTTTTTATCGACGACATGAAt                                                     <  1:474604/81‑1 (MQ=255)                   gTAAACGTAACTCTTTAAGATTATCTGACAGCTCCCAGGCACCATCCACTTTGCCGCTTTTTTTTATCGACGACATGAATa                                                    <  1:474605/81‑1 (MQ=255)                       aCGTAACTCTTTAAGATTATCTGACAGCTCCCAGGCACCATCCACTTTGCCGCTTTTTTTTATCGACGACATGAATAAcgc                                                <  1:474606/81‑1 (MQ=255)                          tAACTCTTTAAGATTATCTGACAGCTCCCAGGCACCATCCACTTTGCCGCTTTTTTTTATCGACGACATGAATAACGCGt                                              >  1:474647/1‑80 (MQ=255)                                  tAAGATTATCTGACAGCTCCCAGGCACCATCCACTTTGCCGCTTTTTTTTATCGACGACATGAATAACGCGTGAGAAATc                                      <  2:474607/80‑1 (MQ=255)                                      aTTATCTGACAGCTCCCAGGCACCATCCACTTTGCCGCTTTTTTTTATCGACGACATGAATAACGCGTGAGAAATCCTGAt                                 >  1:474635/1‑81 (MQ=255)                                              aCAGCTCCCAGGCACCATCCACTTTGCCGCTTTTTTTTATCGACGACATGAATAACGCGTGAGAAATCCTGATCCGGGTc                          <  2:474608/80‑1 (MQ=255)                                              aCAGCTCCCAGGCACCATCCACTTTGCCGCTTTTTTTTATCGACGACATGAATAACGCGTGAGAAATCCTGATCCGGGTc                          <  1:474609/80‑1 (MQ=255)                                                aGCTCCCAGGAACCATCCACTTTGCCGCTTTTTTTTATCGACGACATGAATAACGCGTGAGAAATCCTGATCCGGGTCgag                       <  2:474610/81‑1 (MQ=255)                                                 gCTCCCAGGCACCATCCACTTTGCCGCTTTTTTTTATCGACGACATGAATAACGCGTGAGAAATCCTGATCCGGGTCga                        >  2:474643/1‑79 (MQ=255)                                                       aGGCACCATCCACTTTGCCGCTTTTTTTTATCGACGACATGAATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATg                >  2:474642/1‑81 (MQ=255)                                                           aCCATCCACTTTGCCGCTTTTTTTTATCGACGACATGAATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGgaga            <  2:474611/81‑1 (MQ=255)                                                                      tGCCGCTTTTTTTTATCGACGACATGAATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCAc  <  2:474612/80‑1 (MQ=255)                                                                       gCCGCTTTTTTTTATCGACGACATGAATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCAc  >  1:474645/1‑79 (MQ=255)                                                                            |                                                                          GTTTCGGTTCGAGGTGGCGTAAACGTAACTCTTTAAGATTATCTGACAGCTCCCAGGCACCATCCACTTTGCCGC‑TTTTTTTATCGACGACATGAATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCAC  >  NC_000913/2651891‑2652038

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