Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,704,009	G→A	H19Y (CAT→TAT)	rnc ←	RNase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,704,009	0	G	A	100.0%	104.1 / NA	31	H19Y (CAT→TAT)	rnc	RNase III
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/14);  total (17/14)

GAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATGATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCG  >  NC_000913/2703880‑2704141                                                                                                                                  |                                                                                                                                     gAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAgt                                                                                                                             >  2:653861/1‑139 (MQ=255)  aGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGTAACAGTTCCTGATAATTAAAAgtg                                                                                                                            >  1:853952/1‑139 (MQ=255)   gggAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAgtgt                                                                                                                           >  1:428687/1‑139 (MQ=255)    ggAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTa                                                                                                                          >  2:344129/1‑139 (MQ=255)             gATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGGTGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCtt                                                                                                                 <  2:208213/139‑1 (MQ=255)                  aGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGtt                                                                                                            >  1:902509/1‑139 (MQ=255)                      cATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAg                                                                                                        <  2:902509/139‑1 (MQ=255)                               aaaCGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGattaat                                                                                               <  2:503832/138‑1 (MQ=255)                                             aTAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCa                                                                                 >  2:1084156/1‑139 (MQ=255)                                               agagTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGTAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATg                                                                               <  2:853952/139‑1 (MQ=255)                                                      ccTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGCTGGCGTTCATGCGATAcc                                                                        >  2:686558/1‑139 (MQ=255)                                                      ccTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATAcc                                                                        <  1:80818/139‑1 (MQ=255)                                                      ccTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATAcc                                                                        >  1:694385/1‑139 (MQ=255)                                                      ccTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATAcc                                                                        >  2:707845/1‑139 (MQ=255)                                                      ccTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATAcc                                                                        >  2:740669/1‑139 (MQ=255)                                                      ccTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATAcc                                                                        >  2:52615/1‑139 (MQ=255)                                                                                   aCTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGcac                                           <  1:1084156/139‑1 (MQ=255)                                                                                      gctgGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACg                                        >  2:359065/1‑139 (MQ=255)                                                                                        tgGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGCTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGaa                                      >  1:309473/1‑139 (MQ=255)                                                                                         ggCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGa                                                                                          <  2:62161/86‑1 (MQ=255)                                                                                         ggCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGa                                                                                          >  1:62161/1‑86 (MQ=255)                                                                                         ggCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGaaa                                     <  1:52615/139‑1 (MQ=255)                                                                                         ggCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGaaa                                     <  1:344129/139‑1 (MQ=255)                                                                                          gCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGAAAc                                    >  2:920971/1‑139 (MQ=255)                                                                                              tACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGAAACAGAc                                <  2:501836/139‑1 (MQ=255)                                                                                              tACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGAAACAGAc                                <  2:428821/139‑1 (MQ=255)                                                                                               aCGATGAGTTAATGCCTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGAAACAGAcc                               <  2:476917/139‑1 (MQ=255)                                                                                                             ccTGCTGCAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGtt                 <  1:920971/139‑1 (MQ=255)                                                                                                                 ctgcAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGAAACAg                                  >  2:490873/1‑118 (MQ=255)                                                                                                                 ctgcAACAGTTCCTGATAATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGAAACAg                                  <  1:490873/118‑1 (MQ=255)                                                                                                                            ccTGATAACTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCg  >  2:1005685/1‑139 (MQ=255)                                                                                                                               gataatTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGAAACAGACCTg                             <  2:951176/109‑1 (MQ=255)                                                                                                                                  |                                                                                                                                     GAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGCACTACGATGAGTTAATGCCTGCTGCAACAGTTCCTGATGATTAAAAGTGTAGCCCAGCTTCCGTTGAAGCCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCG  >  NC_000913/2703880‑2704141

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