Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	991,794	C→T	I391I (ATC→ATT)	pepN →	aminopeptidase N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	991,794	0	C	T	100.0%	111.9 / NA	36	I391I (ATC→ATT)	pepN	aminopeptidase N
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/15);  total (21/15)

ACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTC  >  NC_000913/991660‑991929                                                                                                                                       |                                                                                                                                        aCCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTcaca                                                                                                                                     <  2:452564/139‑1 (MQ=255) aCCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTcaca                                                                                                                                     >  1:339327/1‑139 (MQ=255)     tGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCt                                                                                                                                 >  1:661532/1‑139 (MQ=255)      gcgcgcATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTg                                                                                                                                >  1:435496/1‑139 (MQ=255)                     tGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAAc                                                                                                                 >  2:337630/1‑139 (MQ=255)                             aCGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGaa                                                                                                         >  1:577240/1‑139 (MQ=255)                                    cccGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATg                                                                                                  <  2:435496/139‑1 (MQ=255)                                                   gATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAg                                                                                   <  1:337630/139‑1 (MQ=255)                                                     tCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCg                                                                                 >  1:103021/1‑139 (MQ=255)                                                      ccGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGt                                                                                <  2:339327/139‑1 (MQ=255)                                                         cccGGATATGGTCATTGAGATGAACAACTTCTACACCCTGGCCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCat                                                                             <  2:561585/139‑1 (MQ=255)                                                         cccGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCat                                                                             >  2:653840/1‑139 (MQ=255)                                                                         gagaTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGa                                                             <  2:208693/139‑1 (MQ=255)                                                                              gAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCtgt                                                        >  2:106436/1‑139 (MQ=255)                                                                                acaacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTga                                                      >  2:244841/1‑139 (MQ=255)                                                                                      tCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACtt                                                >  1:553068/1‑139 (MQ=255)                                                                                            cccTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCa                                          <  1:653840/139‑1 (MQ=255)                                                                                               tGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGc                                       >  1:543629/1‑139 (MQ=255)                                                                                               tGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGc                                       >  1:680594/1‑139 (MQ=255)                                                                                               tGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGc                                       >  1:97930/1‑139 (MQ=255)                                                                                               tGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGc                                       >  2:755319/1‑139 (MQ=255)                                                                                                 aCCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATg                                                                                                  <  1:545501/78‑1 (MQ=255)                                                                                                 aCCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATg                                                                                                  >  2:545501/1‑78 (MQ=255)                                                                                                 aCCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCa                                                                 >  1:290253/1‑111 (MQ=255)                                                                                                 aCCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCa                                                                 <  2:290253/111‑1 (MQ=255)                                                                                                 aCCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGa                                     >  1:200031/1‑139 (MQ=255)                                                                                                       tACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGc                                                                                  >  2:524585/1‑88 (MQ=255)                                                                                                       tACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGc                                                                                  <  1:524585/88‑1 (MQ=255)                                                                                                       tACGAGAAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGACGATGGacg                               >  2:454439/1‑137 (MQ=255)                                                                                                            gaAGGGCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCg                          <  2:103021/139‑1 (MQ=255)                                                                                                               gggCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGa                                                                                                          >  2:204290/1‑56 (MQ=255)                                                                                                               gggCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGa                                                                                                          <  1:204290/56‑1 (MQ=255)                                                                                                                ggCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGt                                                                                <  1:608231/81‑1 (MQ=255)                                                                                                                ggCGCGGAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGt                                                                                >  2:608231/1‑81 (MQ=255)                                                                                                                      gAAGTGATTCGCATGATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCg                <  1:275275/139‑1 (MQ=255)                                                                                                                                    gATTCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTc  <  2:661532/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        ACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTC  >  NC_000913/991660‑991929

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