Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	155	9,022	T→G	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	155	9,022	0	T	G	100.0%	220.9 / NA	62	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (25/37);  total (25/37)

GTTGCGCTTTTGGTGCGGCTTAGCCGTGTGTTTTCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACTTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGAAATGCAAGGTACCGTTCTTGAA  >  155/8924‑9119                                                                                                   |                                                                                                  gTTGCGCTTTTGGTGCGGCTTAGCCGTGTGTTTTCGGAGTAATGTGCCGAACCTGTTTGTTGAGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTc                                                                                                 >  2:383747/1‑101 (MQ=255)     cgcTTTTGGTGCGGCTTAGCCGTGTGTTTTCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCa                                                                                             >  1:1282961/1‑101 (MQ=255)       cTTTTGGTGCGGCTTAGCCGTGTGTTTTCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCAtt                                                                                           <  1:759348/101‑1 (MQ=255)       cTTTTGGTGCGGCTTAGCCGTGTGTTTTCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCAtt                                                                                           >  1:458375/1‑101 (MQ=255)          ttGGTGCGGCTTAGCCGTGTGTTTTCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATc                                                                                        <  2:510667/101‑1 (MQ=255)           tGGTGCGGCTTAGCCGTGTGTTTTCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCt                                                                                       <  2:186881/101‑1 (MQ=255)            ggTGCGGCTTAGCCGTGTGTTTTCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCtt                                                                                      <  1:1658928/101‑1 (MQ=255)                     tAGCCGTGTGTTTTCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCa                                                                             <  2:1660677/101‑1 (MQ=255)                        ccGTGTGTTTTCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAAt                                                                          >  1:1226199/1‑101 (MQ=255)                         cGTGTGTTTTCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAAtt                                                                         >  1:346207/1‑101 (MQ=255)                          gtgtgtTTTCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTa                                                                        >  1:1062299/1‑101 (MQ=255)                           tgtgtTTTCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTAc                                                                       <  2:1324342/101‑1 (MQ=255)                                 ttCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGt                                                                 <  2:1304123/101‑1 (MQ=255)                                 ttCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGt                                                                 <  1:1093741/101‑1 (MQ=255)                                   cGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGa                                                               >  2:113319/1‑101 (MQ=255)                                   cGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGa                                                               >  1:1086244/1‑101 (MQ=255)                                    ggAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGAt                                                              <  1:1234883/101‑1 (MQ=255)                                     gAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATa                                                             <  2:68220/101‑1 (MQ=255)                                      aGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATAc                                                            <  1:627123/101‑1 (MQ=255)                                        tAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATAccc                                                          >  2:467818/1‑101 (MQ=255)                                         aaTGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATAcccc                                                         >  2:582195/1‑101 (MQ=255)                                         aaTGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATAcccc                                                         <  1:51961/101‑1 (MQ=255)                                          aTGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCa                                                        >  1:1621066/1‑101 (MQ=255)                                          aTGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCa                                                        >  2:834313/1‑101 (MQ=255)                                            gtgCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCaga                                                      >  2:1326150/1‑101 (MQ=255)                                            gtgCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTCATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCaga                                                      >  1:1413514/1‑101 (MQ=255)                                              gacGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAgg                                                    >  1:1439343/3‑101 (MQ=255)                                              gCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAgg                                                    >  1:1439344/1‑101 (MQ=255)                                              gCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAgg                                                    <  1:1515048/101‑1 (MQ=255)                                                cGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGAt                                                  <  2:624761/101‑1 (MQ=255)                                                cGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGAt                                                  >  2:166845/1‑101 (MQ=255)                                                cGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGAt                                                  <  1:77405/101‑1 (MQ=255)                                                  aaCCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTa                                                <  1:372739/101‑1 (MQ=255)                                                  aaCCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTa                                                <  2:415833/101‑1 (MQ=255)                                                     cTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGAt                                             <  1:298850/101‑1 (MQ=255)                                                     cTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGAt                                             <  2:1226194/101‑1 (MQ=255)                                                      tgtttgttGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGCGGATTAGATg                                            >  1:296251/1‑101 (MQ=255)                                                         ttgttgCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGcc                                         <  2:721250/101‑1 (MQ=255)                                                            ttgCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCaaa                                      >  2:1661034/1‑101 (MQ=255)                                                             tgCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAaag                                     <  1:1288464/101‑1 (MQ=255)                                                               cGATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAaagaa                                   <  2:998378/101‑1 (MQ=255)                                                                gATTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAaagaag                                  <  1:1679055/101‑1 (MQ=255)                                                                 aTTTAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAcgaaga                                 <  1:1423482/101‑1 (MQ=255)                                                                  tttAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGAc                                <  1:278394/101‑1 (MQ=255)                                                                  tttAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGAc                                <  1:1149890/101‑1 (MQ=255)                                                                   ttAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACa                               <  2:1086239/101‑1 (MQ=255)                                                                    tAGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACaa                              <  1:467823/101‑1 (MQ=255)                                                                     aGCGCGCAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAat                             <  2:789612/101‑1 (MQ=255)                                                                       cgcgcAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAatat                           >  1:1301469/1‑101 (MQ=255)                                                                         cgcAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTg                         >  2:1542695/1‑101 (MQ=255)                                                                           cAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGaa                       <  1:1075630/101‑1 (MQ=255)                                                                           cAAATCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGaa                       <  2:641943/101‑1 (MQ=255)                                                                            aaaTCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGaaa                      <  1:792969/101‑1 (MQ=255)                                                                            aaaTCTTTACTTATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGaaa                      <  2:971561/101‑1 (MQ=255)                                                                                    ctttATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGAAATGCAAGGt              <  1:1131630/99‑1 (MQ=255)                                                                                      ttATTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGAAATGCAAGGTag            >  2:1131625/1‑100 (MQ=255)                                                                                       tATTTACAGAACGTTGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGAAATGCAAGGTAcc           >  2:252244/1‑101 (MQ=255)                                                                                        aTTTACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGAAATGCAAGGTACCg          >  1:562911/1‑101 (MQ=255)                                                                                           tACAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGAAATGCAAGGTACCGTTc       <  1:1331391/101‑1 (MQ=255)                                                                                             cAGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGAAATGCAAGGTACCGTTCtt     <  2:1397730/101‑1 (MQ=255)                                                                                              aGAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGAAATGCAAGGTACCGTTCTTg    <  2:765920/101‑1 (MQ=255)                                                                                               gAACGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGAAATGCAAGGTACCGTTCTTGa   >  1:1071024/1‑101 (MQ=255)                                                                                                aaCGTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGAAATGCAAGGTACCGTTCTTGaa  >  2:1142925/1‑101 (MQ=255)                                                                                                   |                                                                                                  GTTGCGCTTTTGGTGCGGCTTAGCCGTGTGTTTTCGGAGTAATGTGCCGAACCTGTTTGTTGCGATTTAGCGCGCAAATCTTTACTTATTTACAGAACTTCGGCATTATCTTGCCGGTTCAAATTACGGTAGTGATACCCCAGAGGATTAGATGGCCAAAGAAGACAATATTGAAATGCAAGGTACCGTTCTTGAA  >  155/8924‑9119

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