Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	536,499	C→G	H150D (CAT→GAT)	unknown →	probable glucanotransferase (endo alpha‑1,4 polygalactosaminidase related protein)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	536,499	0	C	G	100.0%	121.7 / NA	43	H150D (CAT→GAT)	unknown	probable glucanotransferase (endo alpha‑1,4 polygalactosaminidase related protein)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (20/23);  total (20/23)

GCTTCAGGGAAAACTGTGTATGGGTATATTAATATTGGAGCAATTGAGAAACACCGCCCGTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATCATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAGGGGATCTTTTTGTAAAAAAATGTATTAATCAAAATACCGTTTTAATTGATGGTATTAACCAAGAAAGTGTTTTTACCACTATTGATTTTAATAAAAAAAAATATGGAATCCAGCAGA  >  1/536265‑536703                                                                                                                                                                                                                                           |                                                                                                                                                                                                             gCTTCAGGGAAAACTGTGTATGGGTATATTAATATTGGAGCAATTGAGAAACACCGCCCGTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATga                                                                                                                                                                                              >  1:48697/1‑251 (MQ=255)  cTTCAGGGAAAACTGTGTATGGGTATATTAATATTGGAGCAATTGAGAAACACCGCCCGTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATgag                                                                                                                                                                                             >  1:795022/1‑251 (MQ=255)        ggAAAACTGTGTATGGGTATATTAATATTGGAGCAATTGAGAAACACCGCCCGTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGatat                                                                                                                                                                                         <  1:408610/249‑1 (MQ=255)              cTGTGTATGGGTATATTAATATTGGAGCAATTGAGAAACACCGCCCGTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAggg                                                                                                                                                                                 <  1:242664/251‑1 (MQ=255)                       ggTATATTAATATTGGAGCAATTGAGAAACACCGCCCGTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGat                                                                                                                                                                                           >  1:379565/1‑232 (MQ=255)                       ggTATATTAATATTGGAGCAATTGAGAAACACCGCCCGTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAAc                                                                                                                                                                         <  1:436216/250‑1 (MQ=255)                          atatTAATATTGGAGCAATTGAGAAACACCGCCCGTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGc                                                                                                                                                                                >  1:145177/1‑240 (MQ=255)                                 tatTGGAGCAATTGAGAAACACCGCCCGTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGAc                                                                                                                                                                <  1:848116/249‑1 (MQ=255)                                          aaTTGAGAAACACCGCCCGTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGat                                                                                                                                                     <  1:743079/251‑1 (MQ=255)                                              gagaAACACCGCCCGTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATAttt                                                                                                                                                                                       <  1:500088/213‑1 (MQ=255)                                              gagaAACACCGCCCGTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGat                                                                                                                                                     >  1:184795/1‑247 (MQ=255)                                                 aaaCACCGCCCGTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTataa                                                                                                                                              <  1:515856/251‑1 (MQ=255)                                                            gTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAggg                                                                                                                                                                                 >  1:153183/1‑205 (MQ=255)                                                            gTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTataat                                                                                                                                             >  1:517684/1‑241 (MQ=255)                                                                        ttttCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACttatt                                                                                                                                 >  1:228588/1‑241 (MQ=255)                                                                                       gCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTataat                                                                                                                                             <  1:822263/214‑1 (MQ=255)                                                                                            gTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAAttt                                                                                                                                                                    >  1:585326/1‑186 (MQ=255)                                                                                                cTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAgg                                                                                                                      <  1:166282/228‑1 (MQ=255)                                                                                                             ccGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAGGGGATCTTTTTGTAAAAAAATGTATTAATCAAAATAc                                                                                  <  1:130351/251‑1 (MQ=255)                                                                                                               ggATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATgg                                                                                                                         >  1:522653/1‑210 (MQ=255)                                                                                                                gATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGAtattat                                                                                                                                                <  1:881238/186‑1 (MQ=255)                                                                                                                     gCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTc                                                                                                                                                                                      <  1:293666/143‑1 (MQ=255)                                                                                                                        aTGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAGGGGATCTTTTTGTAAAAAAATGTATTAATCAAAATACCGTTTTAAtt                                                                        <  1:596878/250‑1 (MQ=255)                                                                                                                          gggTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTc                                                                                                                                         <  1:128980/183‑1 (MQ=255)                                                                                                                                        tctcCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGGGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACttattatt                                                                                                                              >  1:459352/1‑180 (MQ=255)                                                                                                                                         ctcCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTataat                                                                                                                                             <  1:73123/164‑1 (MQ=255)                                                                                                                                           ccGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTc                                                                                                                                         >  1:755491/1‑166 (MQ=255)                                                                                                                                                 tGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTataa                                                                                                                                              <  1:824087/155‑1 (MQ=255)                                                                                                                                                        gTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTc                                                                                                                                         <  1:103323/153‑1 (MQ=255)                                                                                                                                                                           aGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAGGGGATCTTTTTGTAAAAAAATGTATTAATCAAAATACCGTTTTAATTGATGGTATTAACCAAGAAAGTGTTTTTAcc                                          <  1:324900/229‑1 (MQ=255)                                                                                                                                                                                 aaaGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAGGGGATCTTTTTGTAAAAAAATGTATTAATCAAAATACCGTTTTAATTGATGGTATTAACCAAGAAAGTGTTTTTACCACTAt                                     <  1:66624/228‑1 (MQ=255)                                                                                                                                                                                     tATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAGGGGATCTTTTTGTAAAAAAATGTATTAATCAAAATACCGTTTTAATTGATGGTATTAACCAAGAAAGTGTTTTTACCACTATTGATTTTAAT‑AAAAAAAATATGGAAt          <  1:127710/250‑1 (MQ=255)                                                                                                                                                                                           gATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGat                                                                                                                                                     >  1:472050/1‑106 (MQ=255)                                                                                                                                                                                           gATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAGGGGATCTTTTTGTAAAAAAATGTATTAATCAAAATACCGTTTTAATTGATGGTATTAACCAAGAAAGTGTTTTTACCACTATTGATTTTAAT‑AAAAAAAATATGGAATCcagcag   >  1:314439/1‑251 (MQ=255)                                                                                                                                                                                           gATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAGGGGATCTTTTTGTAAAAAAATGTATTAATCAAAATACCGTTTTAATTGATGGTATTAACCAAGAAAGTGTTTTTACCACTATTGATTTTAAT‑‑AAAAAAATATGGAAACCAACAGa  >  1:606918/1‑251 (MQ=255)                                                                                                                                                                                                ggggTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAGGGGATCTTTTTGTAAAAAAATGTAt                                                                                              >  1:786031/1‑156 (MQ=255)                                                                                                                                                                                                 gggTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTataat                                                                                                                                             >  1:861151/1‑108 (MQ=255)                                                                                                                                                                                                      ggATGGTTTTTTTCTTGATAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAGGGGATCTTTTTGTAAAAAAATGTATTAATCAAAATACCGTTTTAATTGATGGTATTAACCAAGAAAGTGTTTTTAc                                           >  1:272385/1‑201 (MQ=255)                                                                                                                                                                                                                        tAATGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAGGGGATCTTTTTGTAAAAAAATGTATTAATCAAAATACCGTTTTAATTGATGGTATTAACCAAGAAAgtg                                                  <  1:557292/176‑1 (MQ=255)                                                                                                                                                                                                                          aTGCAGATGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTaa                                                                                                                                                                          <  1:803807/54‑1 (MQ=255)                                                                                                                                                                                                                                 tGTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAGGGGATCTTTTTGTAAAAAAATGTATTAATCAAAATACCGTTTTAATTGATGGTATTAACCAAGAAAGTGTTTTTACCACTATTGATTTTAAT‑AAAAAAAATATGGAATCcagcag   >  1:656712/1‑213 (MQ=255)                                                                                                                                                                                                                                  gTTTATTATGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAGGGGATCTTTTTGTAAAAAAATGTATTAATCAAAAt                                                                                    <  1:567330/132‑1 (MQ=255)                                                                                                                                                                                                                                    ttattaTGATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGata                                                                                                                                                    >  1:180044/1‑66 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                             GCTTCAGGGAAAACTGTGTATGGGTATATTAATATTGGAGCAATTGAGAAACACCGCCCGTATTATAGTCATTTTCAGTCAGTCATGCTTGGTATCTATGAAACTTGGCCGGATGAGCGATGGGTTGATGTATCCTCTCCGGTATGGCAGCGTTTTATCGTTAATGAGCTAGGAAAAAAGTATGCTGATATGGGGTTGGATGGTTTTTTTCTTGATAATGCAGATGTTTATTATCATTATCAAACAAATGAGATATTTCAAGGGCTTGTAACAATTTTGACTGGGCTTCGATATTATAATATTCCACTTATTATTAATGGAGGGGATCTTTTTGTAAAAAAATGTATTAATCAAAATACCGTTTTAATTGATGGTATTAACCAAGAAAGTGTTTTTACCACTATTGATTTTAATAAAAAAAAATATGGAATCCAGCAGA  >  1/536265‑536703

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