Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,948,728	T→C	M343V (ATG→GTG)	tap ←	methyl‑accepting protein IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,948,728	0	T	C	100.0%	186.6 / NA	48	M343V (ATG→GTG)	tap	methyl‑accepting protein IV
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/23);  total (25/23)

GCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGT  >  REL606/1948631‑1948824                                                                                                  |                                                                                                 gCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGt                                                                                               >  2:1310067/1‑101 (MQ=255)   cAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTAc                                                                                             <  1:1017759/101‑1 (MQ=255)    aGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACt                                                                                            <  2:79391/101‑1 (MQ=255)    aGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACt                                                                                            <  1:541723/101‑1 (MQ=255)    aGAATATTGGTCTGGAAAGAAATTCCGTCGATAGCGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGGGTGAGTCACGGTACt                                                                                            >  2:513614/1‑101 (MQ=37)       atatTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGCGGCGATCTCCTGCATGGGGTGAGTCACGGTACTGAc                                                                                         >  2:950247/1‑101 (MQ=255)        tatTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGAcc                                                                                        >  2:192277/1‑101 (MQ=255)         atTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCt                                                                                       >  1:415736/1‑101 (MQ=255)         atTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGT‑GCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTg                                                                                      <  2:1496676/101‑1 (MQ=255)                tGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACAcc                                                                                <  1:1547155/101‑1 (MQ=255)                tGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACAcc                                                                                >  1:893798/1‑101 (MQ=255)                tGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACAcc                                                                                <  2:1587759/101‑1 (MQ=255)                tGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACAcc                                                                                <  1:45456/101‑1 (MQ=255)                tGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACAcc                                                                                <  1:192277/101‑1 (MQ=255)                        aaTTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCt                                                                        >  2:396727/1‑101 (MQ=255)                         aTTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTg                                                                       >  2:985230/1‑101 (MQ=255)                         aTTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTg                                                                       >  1:1606561/1‑101 (MQ=255)                         aTTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTg                                                                       >  1:804721/1‑101 (MQ=255)                            ccGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTgcgc                                                                    <  1:1566359/101‑1 (MQ=255)                               tCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGt                                                                 >  1:416767/1‑101 (MQ=255)                                   tAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTc                                                             >  1:1387550/1‑101 (MQ=255)                                     aCGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTcgc                                                           >  1:468314/1‑101 (MQ=255)                                       gCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCgg                                                         >  1:1349185/1‑101 (MQ=255)                                       gCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCgg                                                         <  1:1577632/101‑1 (MQ=255)                                       gCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCgg                                                         <  2:468309/101‑1 (MQ=255)                                       gCTGATAATGTCGCCAATTTTTGGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGCCCTGAACACCGCCCGCCTGCGCCGTTGTCGCgg                                                         <  2:1281090/101‑1 (MQ=255)                                         tGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCa                                                       >  2:452761/1‑101 (MQ=255)                                             aaTGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCAtttt                                                   >  1:981485/1‑101 (MQ=255)                                               tGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCAtttttt                                                 <  1:224922/101‑1 (MQ=255)                                                  cGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGcc                                              <  1:1452186/101‑1 (MQ=255)                                                    ccAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAg                                            <  1:1684218/101‑1 (MQ=255)                                                       aTTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTc                                         >  2:341644/1‑101 (MQ=255)                                                        ttttttGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTcc                                        <  1:1485937/101‑1 (MQ=255)                                                                          cGATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTa                      <  1:498697/101‑1 (MQ=255)                                                                           gATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTAt                     >  1:957888/1‑101 (MQ=255)                                                                           gATCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTAt                     >  1:421256/1‑101 (MQ=255)                                                                            ctctcCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATc                    <  2:727722/100‑1 (MQ=255)                                                                            aTCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATc                    >  1:1463724/1‑101 (MQ=255)                                                                            aTCTCCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATc                    >  1:1550127/1‑101 (MQ=255)                                                                             tctcCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCg                   <  2:1318245/101‑1 (MQ=255)                                                                              ctcCTGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCgg                  <  2:981480/101‑1 (MQ=255)                                                                                 cTGCATGGTGTGAGTCACGGTACTGCCCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGt               <  1:985235/101‑1 (MQ=255)                                                                                  tGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGtt              >  2:133977/1‑101 (MQ=255)                                                                                  tGCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGtt              >  1:1586617/1‑101 (MQ=255)                                                                                   gCATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGttt             >  2:1566503/1‑101 (MQ=255)                                                                                    cATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGtttt            <  2:375901/101‑1 (MQ=255)                                                                                    cATGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGtttt            <  1:318676/101‑1 (MQ=255)                                                                                     aTGGTGTGAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTg           >  2:1434903/1‑101 (MQ=255)                                                                                          gtgAGTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTa      >  1:975033/1‑101 (MQ=255)                                                                                              gTCACGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGt  >  2:1431147/1‑101 (MQ=255)                                                                                                  |                                                                                                 GCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGT  >  REL606/1948631‑1948824

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