Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,626,475	C→T	R302R (CGG→CGA)	dcp ←	peptidyl‑dipeptidase Dcp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,626,475	0	C	T	100.0%	99.0 / NA	30	R302R (CGG→CGA)	dcp	peptidyl‑dipeptidase Dcp
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/11);  total (19/11)

AACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCC‑AGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTTGGGTGGTATTCAG  >  NC_000913/1626236‑1626717                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   aaCTTACTACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCC‑AGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCCCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGtt                                                                                                                                                                                                                                         >  2:86032/1‑251 (MQ=255)                  cTGATTCGCGGTACAGAATACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCACGATCAAATTTCTCCCGCCGTACATGTTCCGCATAAAATGCCC‑AGTCCCACGGCTGCGCGCTAACCCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCCCGCACGTTGACGCGCCGCTGGACCAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGgtg                                                                                                                                                                                                                        >  1:510546/1‑251 (MQ=255)                   tGATTCGCGGTACAGAATACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCCCATAAAATGCCC‑AGTCCCACGGCTGCCCGCTAAACCCGCCCCGCTGCTTATCGATAACCGCCTGTCTGGAGGCTAATTCATCGCTTGCACGTTGACCCGCCGCGGGGACAAATTCTCGCATAAAGTTAAGTGCTGCTTCCGgtgt                                                                                                                                                                                                                       >  2:40373/1‑251 (MQ=255)                     aTTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCCTCCAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCC‑AGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGCGGCTAATTCCTCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCCGGGGttt                                                                                                                                                                                                                     >  1:582839/1‑251 (MQ=255)                                                                   aaaTATGGCTTGAGATGCGCCTCATCAAGATCAAATTTCTCCCGACGTACCTGTTCCGCATAAAATGCCC‑AGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCTCCGCTGGAACAATTTCTCGCATAAAGTTAAGGGCTGCTTCAGGTGTTTTTGCCCTCTGATCGGCGATTTTCCTTGCGGCATAATGAGGAAAAcc                                                                                                                                                                       >  2:114400/1‑251 (MQ=255)                                                                      tatGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCC‑AGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCa                                                                                                                                                                      >  1:320800/1‑249 (MQ=255)                                                                                         ccacaaaaaaaatcccccccccccgTACCTGTTCC‑CAAAAAATCCCCCGGCCCCCCGGCGGCCCGTTAACCCCGCCCTGCTTCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCa                                                                                                                                                 <  1:79443/229‑1 (MQ=255)                                                                                                            ccgccgTACCTGTTCCGCATAAAATGCCC‑AGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGtt                                                                                                                              >  1:38699/1‑251 (MQ=255)                                                                                                                   tcctctccccaataaaagcccc‑AGCCCCACGGCTGCGCGCGACACCCTCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGa                                                                                                                        <  2:419485/231‑1 (MQ=255)                                                                                                                        ttCCGAATAAAATGCCC‑AGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAgcgc                                                                                                                  >  2:483322/1‑251 (MQ=255)                                                                                                                                      ccc‑AGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGAt                                                                                                                                                                                                         >  1:381016/1‑150 (MQ=255)                                                                                                                                        c‑AGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGcatc                                                                                                  >  1:273032/1‑251 (MQ=255)                                                                                                                                          aGTCCCACGGATACGCGATAAACCCGCCCTGCTCCGTATCCATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCCCGTTGACGCGCCGCTGGAAAAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCCAGCGGCATAATGAGGCAAACCAAGAAGGGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGAGCGGGTATCATTGGcctcc                                                                                                 >  2:412885/1‑247 (MQ=255)                                                                                                                                            ccccACGGCTGGGCGCAAACCCCCCCCTGCTGCTTCTCGAAAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCTCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCAtt                                                                                               <  1:118041/250‑1 (MQ=255)                                                                                                                                                                 ccGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGcc                                                                          >  1:459760/1‑251 (MQ=255)                                                                                                                                                                     ccTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGATGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTACGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGAAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGc                                                                       >  2:125514/1‑250 (MQ=255)                                                                                                                                                                               aTCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGtt                                                            >  1:17056/1‑251 (MQ=255)                                                                                                                                                                               aTCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGt                                                             <  1:297251/250‑1 (MQ=255)                                                                                                                                                                                    tAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTc                                                        >  2:21116/1‑250 (MQ=255)                                                                                                                                                                                           cTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCg                                                >  2:84014/1‑251 (MQ=255)                                                                                                                                                                                            tGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGc                                               >  1:525308/1‑251 (MQ=255)                                                                                                                                                                                             gTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGGTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTGGCCTGTTGTGCACGGATCGCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATGTTTTTCCGCTCGCGTCCAGCCCGCAATAAACATTTTTTCACGCGTCg                                                <  2:232380/249‑1 (MQ=255)                                                                                                                                                                                              tATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCAc                                             >  2:165016/1‑251 (MQ=255)                                                                                                                                                                                                        aaTTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCAt                                   <  2:17056/251‑1 (MQ=255)                                                                                                                                                                                                            cATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGC‑TGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGAGCGCGCATTTCgg                              <  1:494953/251‑1 (MQ=255)                                                                                                                                                                                                             aTCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCgg                              <  1:21116/251‑1 (MQ=255)                                                                                                                                                                                                                   cgcACGTTGACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAg                          <  1:86032/249‑1 (MQ=255)                                                                                                                                                                                                                           gACGCGCCGCTGGAACAATTTCTCGCATCAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGgttgt                >  2:303533/1‑251 (MQ=255)                                                                                                                                                                                                                           gACGCGCCGCTGGAACAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGgttgt                <  2:446347/251‑1 (MQ=255)                                                                                                                                                                                                                                         aCAATTTCTCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGCATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTTGGGTGGTATTcag  <  1:483322/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   AACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCC‑AGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTTGGGTGGTATTCAG  >  NC_000913/1626236‑1626717

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