Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	208,932	A→C	H322Q (CAT→CAG)	h16_A0195 ←	Putative GTPase (G3E family)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	208,932	0	A	C	100.0%	196.3 / NA	50	H322Q (CAT→CAG)	h16_A0195	Putative GTPase (G3E family)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (24/26);  total (24/26)

GGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCATGGAACACCACTTTGCGGTCAACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAG  >  AM260479/208805‑209012                                                                                                                                |                                                                                 ggcAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGaa                                                                             >  2:851564/1‑133 (MQ=255)      ttCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAcaccac                                                                       >  2:23398/1‑134 (MQ=255)               cccTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTc                                                              <  1:851564/134‑1 (MQ=255)                    gAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTccac                                                           <  2:1036927/132‑4 (MQ=255)                    gAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                             <  2:895507/130‑2 (MQ=255)                    gAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                             >  1:895507/1‑129 (MQ=255)                       gATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGTGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGt                                                      >  2:1014527/1‑134 (MQ=255)                       gATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTccac                                                           <  2:428283/129‑4 (MQ=255)                       gATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTccac                                                           >  1:428283/1‑126 (MQ=255)                       gATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                             >  2:576308/1‑126 (MQ=255)                       gATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                             <  1:576308/127‑2 (MQ=255)                       gATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGCCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGt                                                      >  2:735376/1‑134 (MQ=255)                        aTGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGCCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTc                                                     <  1:735376/134‑1 (MQ=255)                                gcgGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAcacca                                                                        <  1:561189/107‑1 (MQ=255)                                gcgGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAAcacca                                                                        >  2:561189/1‑107 (MQ=255)                                gcgGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCatatat                                              <  2:461193/133‑1 (MQ=255)                                gcgGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCatatat                                              <  1:1014527/133‑1 (MQ=255)                                   gggCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTccac                                                           >  2:1005066/1‑114 (MQ=255)                                   gggCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTccac                                                           <  1:1005066/117‑4 (MQ=255)                                   gggCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCcaca                                                          >  1:256263/1‑118 (MQ=255)                                   gggCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCcaca                                                          <  2:256263/118‑1 (MQ=255)                                                  gATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                             <  1:277820/100‑2 (MQ=255)                                                  gATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                             >  2:277820/1‑99 (MQ=255)                                                  gATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCatatat                                              <  1:6212/115‑1 (MQ=255)                                                  gATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTAGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCatatat                                              >  2:6212/1‑115 (MQ=255)                                                       acacCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                             >  1:390736/1‑94 (MQ=255)                                                       acacCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                             <  2:390736/95‑2 (MQ=255)                                                          ccATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGtttt                   <  1:23398/134‑1 (MQ=255)                                                                 cgTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                        >  2:30248/1‑90 (MQ=255)                                                                 cgTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                        >  2:945437/1‑90 (MQ=255)                                                                 cgTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                        <  1:945437/90‑1 (MQ=255)                                                                 cgTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                        <  1:30248/90‑1 (MQ=255)                                                                 cgTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCatatat                                              >  2:1062512/1‑100 (MQ=255)                                                                 cgTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCatatat                                              <  1:1062512/100‑1 (MQ=255)                                                                       cGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTccac                                                           >  2:289477/1‑78 (MQ=255)                                                                       cGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTccac                                                           <  1:289477/81‑4 (MQ=255)                                                                       cGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACgg      >  2:885166/1‑134 (MQ=255)                                                                         gggTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTccac                                                           >  2:58177/1‑76 (MQ=255)                                                                         gggTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTccac                                                           <  1:58177/79‑4 (MQ=255)                                                                         gggTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                             <  2:1033728/77‑2 (MQ=255)                                                                         gggTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTcc                                                             >  1:1033728/1‑76 (MQ=255)                                                                         gggTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCatatat                                              >  1:886325/1‑92 (MQ=255)                                                                         gggTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCatatat                                              <  2:886325/92‑1 (MQ=255)                                                                           gTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAg  <  2:163550/134‑1 (MQ=255)                                                                                    tccCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                        >  1:882256/1‑71 (MQ=255)                                                                                    tccCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCACAcc                                                        <  2:882256/71‑1 (MQ=255)                                                                                            ggCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCacac                                                         >  1:524095/1‑62 (MQ=255)                                                                                            ggCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCCTGGAACACCACTTTGCGGTCCacac                                                         <  2:524095/62‑1 (MQ=255)                                                                                                                          gACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGttt                    <  1:646522/69‑1 (MQ=255)                                                                                                                          gACCCCCTGGAACACCACTTTGCGGTCCACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGttt                    >  2:646522/1‑69 (MQ=255)                                                                                                                                |                                                                                 GGCAGTTCTCCAGGCCCTTGAGGATGGCGTCGCGGGGCAGGTCCACGCCGATGAACACCATGCGCGTGCCCGGGGTCTCTTCCTCCCACTTGGCGCCGATGTCGCTGCCCATCAGCTGGTGGACCCCATGGAACACCACTTTGCGGTCAACACCGTCCATATATAGTACGCCCTTGTAGCGCAGCAGTTTTTCCCCGTATACGGCCAG  >  AM260479/208805‑209012

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