Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,717,565	A→G	T86A (ACC→GCC)	ghrB →	glyoxylate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,717,565	0	A	G	100.0%	116.5 / NA	34	T86A (ACC→GCC)	ghrB	glyoxylate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/17);  total (17/17)

TGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCA  >  NC_000913/3717438‑3717690                                                                                                                                |                                                                                                                              tGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTgccgcc                                                                                                                          >  2:79771/1‑133 (MQ=255)   ccGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCg                                                                                                                        <  2:63398/133‑1 (MQ=255)        gCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGaaaaa                                                                                                                   <  2:427627/133‑1 (MQ=255)           gAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTc                                                                                                                <  1:438603/133‑1 (MQ=255)               gTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTctgct                                                                                                            <  2:329752/133‑1 (MQ=255)                  tACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGAt                                                                                                         >  1:314432/1‑133 (MQ=255)                    cTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGc                                                                                                       >  2:361514/1‑133 (MQ=255)                       ggTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGcaca                                                                                                    <  1:76663/133‑1 (MQ=255)                       ggTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGcaca                                                                                                    >  1:510977/1‑133 (MQ=255)                          tCAAACGACAATGTATATGCCGCATTGCTGTAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGc                                                                                                 <  1:216261/133‑1 (MQ=255)                                 agaATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGt                                                                                          >  2:77966/1‑133 (MQ=255)                                 agaATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGt                                                                                          >  2:33354/1‑133 (MQ=255)                                   aaTGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTAt                                                                                        >  2:16535/1‑133 (MQ=255)                                               gCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGTCAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCg                                                                            >  1:449387/1‑133 (MQ=255)                                               gCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGTTGCACACACCATCCGTATTAACAGAAACCg                                                                            >  1:145109/1‑133 (MQ=255)                                                cATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGt                                                                           <  1:87834/133‑1 (MQ=255)                                                cATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGt                                                                           <  2:50596/133‑1 (MQ=255)                                                cATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGt                                                                           >  1:281041/1‑133 (MQ=255)                                                 aTTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTc                                                                          >  1:383152/1‑133 (MQ=255)                                                 aTTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTc                                                                          >  2:380698/1‑133 (MQ=255)                                                  ttGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCg                                                                         >  1:500244/1‑133 (MQ=255)                                                      tGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGa                                                                     <  2:383152/133‑1 (MQ=255)                                                                    aaaCTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGgcgc                                                       <  1:79771/133‑1 (MQ=255)                                                                    aaaCTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGgcgc                                                       <  2:510977/133‑1 (MQ=255)                                                                      aCTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTg                                                     <  2:260130/133‑1 (MQ=255)                                                                                 caTCAACGATCTCAGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGTACACGCCAACCGTATTAACGGAAACCGTTGCCTATAAGCTGATGTAGCTGCTGTTGTCTAc                                          >  2:171872/1‑133 (MQ=255)                                                                                         aTCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCtcgtc                                  >  1:261652/1‑133 (MQ=255)                                                                                           ctcCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCgg                                >  2:363264/1‑133 (MQ=255)                                                                                            tcCGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCggg                               <  2:314432/133‑1 (MQ=255)                                                                                              cGTCGGCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGtt                             <  1:380698/133‑1 (MQ=255)                                                                                                tCGGCTATTACAATTATGATGTTGATGCGCATGCCGCACGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTtgt                           <  1:227246/133‑1 (MQ=255)                                                                                                   gCTATGACAATTTTGATGTCGATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGa                        <  2:281041/133‑1 (MQ=255)                                                                                                                       gATGCGCTTGCCGCCCGAAAATTTCTGCTGATGCACACGGCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAg    >  2:297346/1‑133 (MQ=255)                                                                                                                       gATGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAg    >  2:457827/1‑133 (MQ=255)                                                                                                                         tGCGCTTGCCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCa  <  1:77966/133‑1 (MQ=255)                                                                                                                                |                                                                                                                              TGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCA  >  NC_000913/3717438‑3717690

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