Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	85	2,871	G→T	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	85	2,871	0	G	T	100.0%	249.2 / NA	72	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (31/41);  total (31/41)

TTTTTCCACGAACAGCAGTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCAGAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGCACATCAACCTCAAA  >  85/2773‑2968                                                                                                   |                                                                                                  tttttCCACGAACAGCAGTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATaa                                                                                                 <  2:805663/101‑1 (MQ=255)   tttCCACGAACAGCAGTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAAt                                                                                               <  1:779528/101‑1 (MQ=255)   tttCCACGAACAGCAGTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAAt                                                                                               <  1:1395965/101‑1 (MQ=255)      ccACGAACAGCAGTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACt                                                                                            <  1:1136949/101‑1 (MQ=255)      ccACGAACAGCAGTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACt                                                                                            <  1:1542054/101‑1 (MQ=255)      ccACGAACAGCAGTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACt                                                                                            <  1:1198708/101‑1 (MQ=255)       cACGAACAGCAGTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTg                                                                                           <  1:1336761/101‑1 (MQ=255)        aCGAACAGCAGTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTgg                                                                                          <  2:1252793/101‑1 (MQ=255)          gAACAGCAGTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTc                                                                                        <  2:1515234/101‑1 (MQ=255)           aaCAGCAGTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCa                                                                                       <  2:285858/101‑1 (MQ=255)             cagcagTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATc                                                                                     <  2:1077135/101‑1 (MQ=255)              agcagTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATcc                                                                                    <  1:564948/101‑1 (MQ=255)              agcagTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATcc                                                                                    >  2:1701903/1‑101 (MQ=255)                  gTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGccc                                                                                >  2:1109698/1‑101 (MQ=255)                   tCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCt                                                                               >  1:1337618/1‑101 (MQ=255)                    cACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCtt                                                                              >  1:326184/1‑101 (MQ=255)                    cACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCtt                                                                              <  2:1119052/101‑1 (MQ=255)                     aCGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTc                                                                             <  1:717891/101‑1 (MQ=255)                     aCGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTc                                                                             >  1:1335344/1‑101 (MQ=255)                      cGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCt                                                                            >  2:968146/1‑101 (MQ=255)                         aGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGc                                                                         >  1:1224611/1‑101 (MQ=255)                           ccGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCaa                                                                       >  2:173610/1‑101 (MQ=255)                             ggCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATg                                                                     <  1:95087/101‑1 (MQ=255)                               cACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGAc                                                                   <  1:651494/101‑1 (MQ=255)                               cACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGAc                                                                   <  2:975337/101‑1 (MQ=255)                               cACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGAc                                                                   <  1:1109585/101‑1 (MQ=255)                               cACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGAc                                                                   <  2:579096/101‑1 (MQ=255)                                       aGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTa                                                           >  2:112389/1‑101 (MQ=255)                                        gTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAg                                                          >  1:380653/1‑101 (MQ=255)                                         ttgattgaTTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGa                                                         >  1:764317/1‑101 (MQ=255)                                         ttgattgaTTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGa                                                         >  1:1549179/1‑101 (MQ=255)                                           gattgattATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATc                                                       >  2:89038/1‑101 (MQ=255)                                           gattgattATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATc                                                       <  1:288843/101‑1 (MQ=255)                                            attgattATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCa                                                      >  1:1164039/1‑101 (MQ=255)                                            attgattATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCa                                                      <  1:373346/101‑1 (MQ=255)                                            attgattATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCa                                                      <  2:68256/101‑1 (MQ=255)                                                  tattTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCAc                                                <  1:348409/101‑1 (MQ=255)                                                   attTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCAcg                                               <  1:1289753/101‑1 (MQ=255)                                                   attTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCAcg                                               <  2:1287956/101‑1 (MQ=255)                                                             tatGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCAcgtaatccgtaa                                     >  1:889158/1‑101 (MQ=255)                                                              atGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCAcgtaatccgtaat                                    <  2:195087/101‑1 (MQ=255)                                                              atGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCAcgtaatccgtaat                                    <  1:4078/101‑1 (MQ=255)                                                              atGATTACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCAcgtaatccgtaat                                    >  1:438265/1‑101 (MQ=255)                                                                  ttACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGt                                <  1:1392516/101‑1 (MQ=255)                                                                  ttACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGt                                <  1:350718/101‑1 (MQ=255)                                                                  ttACGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGt                                <  1:1413690/101‑1 (MQ=255)                                                                    aCGGGAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGa                              <  1:785688/101‑1 (MQ=255)                                                                       ggAGACGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAAc                           >  2:370511/1‑101 (MQ=255)                                                                        gagaCGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACt                          >  2:1396320/1‑101 (MQ=255)                                                                        gagaCGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACt                          <  1:1387696/101‑1 (MQ=255)                                                                           aCGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGt                       <  2:1692525/101‑1 (MQ=255)                                                                            cGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGtt                      <  1:1517852/101‑1 (MQ=255)                                                                            cGCCGTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGtt                      <  1:1291372/101‑1 (MQ=255)                                                                               tgTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAAc                   >  2:1069226/2‑101 (MQ=255)                                                                               tgTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAAc                   >  2:98256/2‑101 (MQ=255)                                                                                gTAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACt                  <  1:1701909/101‑1 (MQ=255)                                                                                 tAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTg                 >  2:36663/1‑101 (MQ=255)                                                                                 tAGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTg                 >  2:1223602/1‑101 (MQ=255)                                                                                  aGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGc                >  1:1007652/1‑101 (MQ=255)                                                                                  aGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGc                >  2:1210808/1‑101 (MQ=255)                                                                                  aGACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGc                >  2:1241190/1‑101 (MQ=255)                                                                                   gACGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGca               <  2:584976/101‑1 (MQ=255)                                                                                    aCGTCAGTCACTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGcac              <  2:1549173/101‑1 (MQ=255)                                                                                      gtcagtcaCTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGCACAt            >  1:863424/1‑101 (MQ=255)                                                                                         agtcaCTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGCACATCaa         >  2:1111810/1‑101 (MQ=255)                                                                                          gtcaCTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGGAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGCACATCAAc        >  2:558251/1‑101 (MQ=255)                                                                                            caCTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGCACATCAACCt      >  2:910520/1‑101 (MQ=255)                                                                                            caCTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGCACATCAACCt      <  2:1567738/101‑1 (MQ=255)                                                                                              cTTCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGCACATCAACCTCa    >  1:832600/1‑101 (MQ=255)                                                                                               ttCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGCACATCAACCTCaa   >  2:1441804/1‑101 (MQ=255)                                                                                                tCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGCACATCAACCTCaaa  <  2:861172/101‑1 (MQ=255)                                                                                                tCATAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGCACATCAACCTCaaa  <  1:940444/101‑1 (MQ=255)                                                                                                   |                                                                                                  TTTTTCCACGAACAGCAGTCACGTAGCCGGCACGACTAAGTTGATTGATTATTTTGACCATATGATTACGGGAGACGCCGTAGACGTCAGTCACTTCAGAAATACTGGTCATCCGCCCTTCTGGCAATGACGCCATGTAGATCAGCGCACGTAATCCGTAATCAGTGAAACTCGTTAACTGCACATCAACCTCAAA  >  85/2773‑2968

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