Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,648,666	C→A	A1208S (GCG→TCG)	yfhM ←	alpha2‑macroglobulin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,648,666	0	C	A	100.0%	133.3 / NA	37	A1208S (GCG→TCG)	yfhM	alpha2‑macroglobulin
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (18/19);  total (18/19)

CATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATAC  >  NC_000913/2648589‑2648742                                                                              |                                                                             cATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACt                                                                            <  2:473858/80‑1 (MQ=255) cATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACt                                                                            <  1:473859/80‑1 (MQ=255)  aTTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGc                                                                          <  1:473860/81‑1 (MQ=255)      gcagcaAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTCTGATGCCCAACGACTGCAGtt                                                                      >  1:473898/1‑81 (MQ=255)        agcaAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTgg                                                                    >  2:473911/1‑81 (MQ=255)            aaCGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTg                                                                >  1:473895/1‑81 (MQ=255)             aCGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCGCTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGc                                                               >  1:473962/1‑81 (MQ=255)               gggAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGt                                                             >  1:473910/1‑81 (MQ=255)               gggAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGt                                                             >  2:473893/1‑81 (MQ=255)                ggAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGtt                                                            <  1:473861/81‑1 (MQ=255)                  aaaTGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTgg                                                          <  1:473862/81‑1 (MQ=255)                   aaTGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGt                                                         <  2:473863/81‑1 (MQ=255)                          atatCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGt                                                  >  1:473900/1‑81 (MQ=255)                           tatCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGt                                                  >  1:473904/1‑80 (MQ=255)                           tatCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTg                                                 <  1:473864/81‑1 (MQ=255)                            atCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGa                                                >  1:473905/1‑81 (MQ=255)                            atCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGa                                                >  1:473897/1‑81 (MQ=255)                             tCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGa                                                <  1:473865/80‑1 (MQ=255)                               gACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACg                                              >  1:473902/1‑80 (MQ=255)                                  cGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGaaa                                          >  1:473899/1‑81 (MQ=255)                                   gATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGaaa                                          >  2:473907/1‑80 (MQ=255)                                   gATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAAc                                         <  2:473866/81‑1 (MQ=255)                                        aCGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGcc                                    <  1:473868/81‑1 (MQ=255)                                        aCGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGcc                                    <  2:473867/81‑1 (MQ=255)                                         cGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCg                                   >  1:473915/1‑81 (MQ=255)                                            cGTTTCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTg                                >  2:473914/1‑81 (MQ=255)                                              tttCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGc                              <  2:473869/81‑1 (MQ=255)                                               ttCTCATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCg                             >  2:473909/1‑81 (MQ=255)                                                tctcATCACTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCg                             <  2:473870/80‑1 (MQ=255)                                                       aCTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTc                     <  2:473872/81‑1 (MQ=255)                                                       aCTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTc                     <  2:473873/81‑1 (MQ=255)                                                       aCTGTCGCCTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTc                     <  1:473871/81‑1 (MQ=255)                                                              ccTTTGATGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGTAAACAGGCCGCTGGCGGTTTGCTCAAGACAg              <  2:473874/81‑1 (MQ=255)                                                                    aTGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTa         <  1:473875/80‑1 (MQ=255)                                                                     tGCCCAACGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACg       >  2:473908/1‑81 (MQ=255)                                                                          aaCGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATAc  >  1:473917/1‑81 (MQ=255)                                                                          aaCGACTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATAc  <  1:473876/81‑1 (MQ=255)                                                                              |                                                                             CATTTGCAGCAAACGGGAAATGCCGATATCGACCGATGCACGGCGTTTCTCATCACTGTCGCCTTTGATGCCCAACGCCTGCAGTTGGGCTGCGTTGGTATAAAGTGACGGAAACAGGCCGCTGGCGGTTTGCTCAAGACAGCCGTACGGATAC  >  NC_000913/2648589‑2648742

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