Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	2	58,599	G→A	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	2	58,599	0	G	A	100.0%	223.1 / NA	63	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/29);  total (34/29)

ACCACCAGCTCTACGGAACATCGCTGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAGCCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGTGGGTTCGCTGCTTT  >  2/58505‑58694                                                                                               |                                                                                                accaccAGCTCTACGGAACATCGCTGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTg                                                                                           >  1:144381/1‑101 (MQ=255)     ccaGCTCTACGGAATATCGCTGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTg                                                                                       <  2:1132682/101‑1 (MQ=255)           ctACGGAACATCGCTGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAg                                                                                 <  1:49474/101‑1 (MQ=255)            tACGGAACATCGCTGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGt                                                                                <  1:1070840/101‑1 (MQ=255)                gAACATCGCTGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACgg                                                                            <  1:1193960/101‑1 (MQ=255)                 aaCATCGCTGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGa                                                                           <  1:860155/101‑1 (MQ=255)                  aCATCGCTGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGAc                                                                          >  2:1230180/1‑101 (MQ=255)                     tCGCTGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAAt                                                                       <  2:875822/101‑1 (MQ=255)                      cGCTGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAAtt                                                                      >  1:417552/1‑101 (MQ=255)                       gCTGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAAttt                                                                     >  2:808367/1‑101 (MQ=255)                        cTGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAACGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATttg                                                                    >  1:1566735/1‑101 (MQ=255)                         tGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATttgt                                                                   >  1:689537/1‑101 (MQ=255)                         tGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATttgt                                                                   >  2:281343/1‑101 (MQ=255)                          gTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATttgtt                                                                  >  2:1056459/1‑101 (MQ=255)                              aaTTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAg                                                              <  1:31493/101‑1 (MQ=255)                              aaTTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAg                                                              >  1:1197809/1‑101 (MQ=255)                               aTTATCGATTTACCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAgg                                                             >  2:1615609/1‑101 (MQ=255)                                  aTCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCAt                                                          <  1:1625894/101‑1 (MQ=255)                                  aTCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCAt                                                          <  1:68104/101‑1 (MQ=255)                                    cGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGt                                                        <  1:451720/101‑1 (MQ=255)                                     gATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTc                                                       >  1:616761/1‑101 (MQ=255)                                          aGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATaca                                                  >  2:473244/1‑101 (MQ=255)                                          aGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATaca                                                  >  2:724483/1‑101 (MQ=255)                                           gCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATacac                                                 <  2:1229864/101‑1 (MQ=255)                                           gCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTACAGGCATGTCATacac                                                 <  1:742328/101‑1 (MQ=255)                                             cAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTa                                               >  2:1618315/1‑101 (MQ=255)                                              aGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTAt                                              <  1:118922/101‑1 (MQ=255)                                              aGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTAt                                              >  1:1672357/1‑101 (MQ=255)                                               gTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATg                                             >  1:1369986/1‑101 (MQ=255)                                               gTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATg                                             >  2:1554841/1‑101 (MQ=255)                                               gTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATg                                             <  2:1717267/101‑1 (MQ=255)                                                ttATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGt                                            >  1:1041946/1‑101 (MQ=255)                                                  ataCCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCg                                          >  2:1061669/1‑101 (MQ=255)                                                  ataCCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCg                                          <  2:253316/101‑1 (MQ=255)                                                    aCCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGAt                                        <  2:1643212/101‑1 (MQ=255)                                                       ttATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCaa                                     >  2:1238696/1‑101 (MQ=255)                                                       ttATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCaa                                     <  1:1459932/101‑1 (MQ=255)                                                          tatGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAact                                  >  2:596260/1‑101 (MQ=255)                                                           atGGCTATATCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAacta                                 >  1:970124/1‑101 (MQ=255)                                                                tataTCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAg                            >  2:165362/1‑101 (MQ=255)                                                                 atatCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGt                           <  2:1076637/101‑1 (MQ=255)                                                                  tatCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTa                          >  1:1674444/1‑101 (MQ=255)                                                                  tatCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTa                          <  1:625636/101‑1 (MQ=255)                                                                  tatCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTa                          >  2:1409434/1‑101 (MQ=255)                                                                   atCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTAt                         <  1:1505527/101‑1 (MQ=255)                                                                   atCTCAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTAt                         >  2:1151768/1‑101 (MQ=255)                                                                     ctcAATGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCt                       <  1:59302/101‑1 (MQ=255)                                                                         aTGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTacac                   >  2:1462860/1‑101 (MQ=255)                                                                         aTGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTacac                   <  1:1144059/101‑1 (MQ=255)                                                                          tGGCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACAcc                  <  1:139993/101‑1 (MQ=255)                                                                            gCAACGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGt                <  1:1230185/101‑1 (MQ=255)                                                                              aaCGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGTgg              >  2:446499/1‑101 (MQ=255)                                                                              aaCGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGTgg              >  1:1560363/1‑101 (MQ=255)                                                                                cGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGTGGGt            <  2:797446/101‑1 (MQ=255)                                                                                cGTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGTGGGt            <  1:248963/101‑1 (MQ=255)                                                                                 gTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGTGGGtt           >  2:1646053/1‑101 (MQ=255)                                                                                 gTTTTATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGTGGGtt           >  1:185061/1‑101 (MQ=255)                                                                                  ttttATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGTGGGTTc          <  2:1672351/101‑1 (MQ=255)                                                                                   tttATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGTGGGTTCg         >  1:1639629/1‑101 (MQ=255)                                                                                    ttATCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGTGGGTTCgc        >  1:1353750/1‑101 (MQ=255)                                                                                      aTCTGGTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGTGGGTTCgctg      >  2:1265427/1‑101 (MQ=255)                                                                                          ggTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGTGGGTTCGCTGCttt  <  1:1618321/101‑1 (MQ=255)                                                                                          ggTTAACCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGTGGGTTCGCTGCttt  <  1:1615615/101‑1 (MQ=255)                                                                                               |                                                                                                ACCACCAGCTCTACGGAACATCGCTGTTCAATTATCGATTTAGCCAGTTATACCTTATATGGCTATATCTCAATGGCAACGTTTTATCTGGTTAGCCAGTGGCTGACACAGTACGGACAATTTGTTGCAGGCATGTCATACACTATGTCGATCAAACTACTCAGTATCTACACCGTGGGTTCGCTGCTTT  >  2/58505‑58694

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