Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	4	36,757	G→C	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	4	36,757	0	G	C	100.0%	156.2 / NA	44	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (24/20);  total (24/20)

GGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAAC  >  4/36664‑36853                                                                                              |                                                                                                 gggTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCa                                                                                           >  1:1379619/1‑101 (MQ=255) gggTCTGATCGGGATCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGGAATCGCTGGCa                                                                                           >  1:529757/1‑101 (MQ=255)       gATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTgcgc                                                                                     >  2:330146/1‑101 (MQ=255)        aTCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGcgcg                                                                                    >  2:1540528/1‑101 (MQ=255)         tCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCgg                                                                                   <  1:1083992/101‑1 (MQ=255)               tCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAg                                                                             <  2:703928/101‑1 (MQ=255)               tCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAg                                                                             <  1:404242/101‑1 (MQ=255)                cTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGa                                                                            >  2:566844/1‑101 (MQ=255)                 tACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGaa                                                                           >  1:1052363/1‑101 (MQ=255)                  aCCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAAc                                                                          <  1:386173/101‑1 (MQ=255)                     aGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCGCTGCGCGGTTCCAGAACAat                                                                       <  2:865113/101‑1 (MQ=255)                     aGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAat                                                                       <  2:620293/101‑1 (MQ=255)                       cGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAatat                                                                     <  2:95361/101‑1 (MQ=255)                        gTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAAtata                                                                    <  2:397381/101‑1 (MQ=255)                        gTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAAtata                                                                    >  1:238095/1‑101 (MQ=255)                        gTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAAtata                                                                    >  1:1598824/1‑101 (MQ=255)                         tCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAg                                                                   <  2:1247546/101‑1 (MQ=255)                                    tttGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCaa                                                        >  1:1110406/1‑101 (MQ=255)                                        ccgccCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAacg                                                    <  2:140612/101‑1 (MQ=255)                                            ccGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACgg                                                >  2:1144635/1‑101 (MQ=255)                                             cGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACggg                                               <  1:433963/101‑1 (MQ=255)                                               gCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGgcg                                             <  2:66988/101‑1 (MQ=255)                                                 gTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCa                                           >  1:1079638/1‑101 (MQ=255)                                                 gTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCa                                           >  1:946850/1‑101 (MQ=255)                                                 gTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCa                                           >  2:1090812/1‑101 (MQ=255)                                                           caccaGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGgttgt                                 >  1:221672/1‑101 (MQ=255)                                                           caccaGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGgttgt                                 <  1:203738/101‑1 (MQ=255)                                                                ggCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGa                            <  2:1383399/101‑1 (MQ=255)                                                                  cGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATc                          >  1:1422254/1‑101 (MQ=255)                                                                  cGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATc                          <  2:492821/101‑1 (MQ=255)                                                                   gggATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATcg                         <  1:1087914/101‑1 (MQ=255)                                                                    ggATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATcgc                        >  1:1205503/1‑101 (MQ=255)                                                                    ggATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATcgc                        <  2:118011/101‑1 (MQ=255)                                                                     gATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATcgcg                       <  1:1289029/101‑1 (MQ=255)                                                                      aTTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGcga                      >  1:1645584/1‑101 (MQ=255)                                                                      aTTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGcga                      >  1:375070/1‑101 (MQ=255)                                                                          ttGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGaa                  >  1:1653202/1‑101 (MQ=255)                                                                           tGATCGCTTAACAGGTCATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCTCCGAGCAAAACGACGGGCTCACTGTGGTTGTTCAGATCGCGACGAAt                 <  2:435794/101‑1 (MQ=255)                                                                             aTCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCg               >  1:996185/1‑101 (MQ=255)                                                                             aTCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCg               >  1:502666/1‑101 (MQ=255)                                                                                gCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCa            <  1:956332/101‑1 (MQ=255)                                                                                gCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCa            >  2:783542/1‑101 (MQ=255)                                                                                gCTTAACAGGGAATCGCTGGCACTGAGCGGTTCCAGAACAATATAGCCCCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCa            >  2:1282551/1‑101 (MQ=255)                                                                                   tAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCAtct         >  2:1585311/1‑101 (MQ=255)                                                                                          tAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAAc  >  1:670538/1‑101 (MQ=255)                                                                                          tAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAAc  <  1:992471/101‑1 (MQ=255)                                                                                              |                                                                                                 GGGTCTGATCGGGTTCTACCAGCGTCGCAGCCAGTTTTGCCGCCCGGCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGACGGGCGCACTGTGGTTGTTCAGATCGCGACGAATCGTCATCTCCGAAAC  >  4/36664‑36853

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