Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,702,757	C→T	F175F (TTC→TTT)	add →	adenosine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,702,757	0	C	T	100.0%	99.7 / NA	28	F175F (TTC→TTT)	add	adenosine deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/13);  total (15/13)

GAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTC  >  NC_000913/1702518‑1702975                                                                                                                                                                                                                                                |                                                                                                                                                                                                                           gAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCt                                                                                                                                                                                                                 >  1:245638/1‑251 (MQ=255)    atGCGTTTTTCACCAGGCTACATGGCACTGGCACATCAGCTGCCTGTTGCTGGTGTTGTCGAAGCGGTGCTCGATGGCGTACGGGAAGGTTGCCGCACCTTTGGTGTGCAGCCTAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGACGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGCGACCAGCTTACCACACTTGATTTAGCCGGTGTTGAACTTGGTTTTCCgggtctcctttt                                                                                                                                                                                                              >  2:230549/2‑242 (MQ=255)     tGCGTTTTTCACCAGGCTACATGGCACTGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTc                                                                                                                                                                                                             >  2:206433/1‑251 (MQ=255)             tCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCCGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTtctc                                                                                                                                                                                                      >  2:123895/1‑250 (MQ=255)             tCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGCTTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTtctc                                                                                                                                                                                                      >  2:130283/1‑250 (MQ=255)                 cAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTc                                                                                                                                                                                                 <  1:114090/251‑1 (MQ=255)                         atggcaatggcaCATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCTCCGTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACcgcgc                                                                                                                                                                                         <  2:239932/251‑1 (MQ=255)                                                   aGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCg                                                                                                                                                               >  1:289291/1‑251 (MQ=255)                                                      gggTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTcc                                                                                                                                                            <  1:68975/251‑1 (MQ=255)                                                                gAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCg                                                                                                                                                   >  1:202140/1‑250 (MQ=255)                                                                     ggTGAGCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTg                                                                                                                                             <  2:320182/251‑1 (MQ=255)                                                                        gATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTgcc                                                                                                                                           <  2:95738/250‑1 (MQ=255)                                                                                 cGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCAGTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGaaa                                                                                                                                 <  1:123895/251‑1 (MQ=255)                                                                                                                cagacGAAGCTTATCGGCATTATTAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCGCTTGATTTAGCCGGTGCTGAACTTGGTTTTCCGGAAGGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGCTGCGTGCGGGCATATTTCCGTCCATGCAGGCGACGCTGCCGGGCCGGGAATCATCGGGCAGGCGATTCGttatcttggtgc                                                                                                  >  2:252253/5‑240 (MQ=255)                                                                                                                      aaGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCg                                                                                            >  1:207927/1‑251 (MQ=255)                                                                                                                                 cATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATg                                                                                 <  1:27984/251‑1 (MQ=255)                                                                                                                                          ccGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCATTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGCGCGTATTGGACATGGCGTaaa                                                                          >  1:372383/1‑249 (MQ=255)                                                                                                                                           cGGACCTTCGGCGAAGCCGCCTGGCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTTCCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTTCCGTCCATGCAGGCGAAGCTGCCGGGCCTGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGc                                                                       <  2:381408/251‑1 (MQ=255)                                                                                                                                            ggACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGCTTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGcc                                                                      >  2:401864/1‑251 (MQ=255)                                                                                                                                                                      cAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGAttt                                            <  1:176547/251‑1 (MQ=255)                                                                                                                                                                        agagCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGCCATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTtc                                          >  2:26359/1‑251 (MQ=255)                                                                                                                                                                                  gCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCgggc                                 >  2:348409/1‑247 (MQ=255)                                                                                                                                                                                                  gTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCACCAAATTGGTATTGa                 >  1:397429/1‑250 (MQ=255)                                                                                                                                                                                                  gTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGAACTGGGTTCGGAGCTTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGaa                <  1:39559/251‑1 (MQ=255)                                                                                                                                                                                                        aGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGt          <  1:206433/251‑1 (MQ=255)                                                                                                                                                                                                        aGATTACCGCACTTGATTTAGCAGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGt          >  2:46347/1‑251 (MQ=255)                                                                                                                                                                                                            tACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTg       <  1:401864/250‑1 (MQ=255)                                                                                                                                                                                                                gCACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTGCTCACGTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACGGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTc  <  2:372383/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                           GAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCGCACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTC  >  NC_000913/1702518‑1702975

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