Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,863,004	T→C	R78R (CGT→CGC)	gapA →	glyceraldehyde‑3‑phosphate dehydrogenase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,863,004	0	T	C	100.0%	247.6 / NA	68	R78R (CGT→CGC)	gapA	glyceraldehyde‑3‑phosphate dehydrogenase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/35);  total (33/35)

GACCTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGTGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTGGTTATGACTGGTC  >  NC_000913/1862870‑1863134                                                                                                                                       |                                                                                                                                   gACCTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATc                                                                                                                                <  1:1775558/139‑1 (MQ=255)      gTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGc                                                                                                                            >  2:1545602/1‑138 (MQ=255)      gTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGc                                                                                                                            <  1:1545602/138‑1 (MQ=255)                 aTTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGaaa                                                                                                                  <  1:848549/137‑1 (MQ=255)                 aTTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGaaa                                                                                                                  >  2:848549/1‑137 (MQ=255)                 aTTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATg                                                                                                                >  1:946663/1‑139 (MQ=255)                 aTTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATg                                                                                                                >  2:555102/1‑139 (MQ=255)                     cATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGAc                                                                                                            <  2:1295835/139‑1 (MQ=255)                      aTGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACg                                                                                                           >  1:974409/1‑139 (MQ=255)                       tGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGa                                                                                                          <  2:74169/139‑1 (MQ=255)                         gCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAg                                                                                                        <  2:75608/139‑1 (MQ=255)                              cATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGgt                                                                                                   >  2:1180986/1‑139 (MQ=255)                                    gAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGAc                                                                                             >  2:953604/1‑139 (MQ=255)                                     aaaTATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACg                                                                                            >  2:464840/1‑139 (MQ=255)                                          tGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTc                                                                                       <  1:953604/139‑1 (MQ=255)                                           gACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAg                                                                                                        >  2:881864/1‑121 (MQ=255)                                           gACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAg                                                                                                        <  1:881864/121‑1 (MQ=255)                                           gACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCg                                                                                      >  1:1687791/1‑139 (MQ=255)                                             cTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCt                                                                                    >  1:1415837/1‑139 (MQ=255)                                               ccACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGa                                                                                  <  1:1400359/139‑1 (MQ=255)                                               ccACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGa                                                                                  <  1:417847/139‑1 (MQ=255)                                               ccACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGa                                                                                  >  1:787182/1‑139 (MQ=255)                                                  cTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGa                                                                                                             <  1:702932/109‑1 (MQ=255)                                                  cTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGa                                                                                                             >  2:702932/1‑109 (MQ=255)                                                  cTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGc                                                                               >  1:1612947/1‑139 (MQ=255)                                                         ccGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                                                        <  1:555102/139‑1 (MQ=255)                                                            tttCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGAc                                                                                                            <  2:404881/100‑1 (MQ=255)                                                            tttCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGAc                                                                                                            >  1:404881/1‑100 (MQ=255)                                                            tttCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTg                                                                     >  1:1348956/1‑139 (MQ=255)                                                            tttCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTg                                                                     >  1:163355/1‑139 (MQ=255)                                                                     tACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCctgact                                                            >  2:1250366/1‑139 (MQ=255)                                                                            gAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGaaa                                                     <  2:946663/139‑1 (MQ=255)                                                                               gTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGtt                                                                                          >  2:1807546/1‑99 (MQ=255)                                                                               gTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGtt                                                                                          <  1:1807546/99‑1 (MQ=255)                                                                               gTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTg                                                  >  1:1342597/1‑139 (MQ=255)                                                                                  aaaGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                                                        <  1:1665382/114‑1 (MQ=255)                                                                                  aaaGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                                                        >  2:1665382/1‑114 (MQ=255)                                                                                       cGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGaa                                                                                                                   <  2:656764/66‑1 (MQ=255)                                                                                       cGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGaa                                                                                                                   >  1:656764/1‑66 (MQ=255)                                                                                       cGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTaaa                                          <  1:464840/139‑1 (MQ=255)                                                                                        ggTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAac                                         <  2:1687791/139‑1 (MQ=255)                                                                                         gTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAaca                                        <  2:1415837/139‑1 (MQ=255)                                                                                         gTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAaca                                        >  2:36897/1‑139 (MQ=255)                                                                                          tCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGAC‑GACGAAACTGCTCGTAAAcaca                                      >  1:870235/1‑139 (MQ=255)                                                                                               tGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCAc                                  >  1:1489750/1‑139 (MQ=255)                                                                                                gATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCAcc                                 >  2:452639/1‑139 (MQ=255)                                                                                                 aTCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCg                                >  2:1749433/1‑139 (MQ=255)                                                                                                  tCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATggg                                                                                                              <  1:650852/60‑1 (MQ=255)                                                                                                  tCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATggg                                                                                                              >  2:650852/1‑60 (MQ=255)                                                                                                    gTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCgg                                                                                                                             <  1:1137644/43‑1 (MQ=255)                                                                                                    gTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCgg                                                                                                                             >  2:1137644/1‑43 (MQ=255)                                                                                                      tAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTg                                                  >  2:594407/1‑116 (MQ=255)                                                                                                      tAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTg                                                  <  1:594407/116‑1 (MQ=255)                                                                                                           gTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCgaa                      <  1:1250366/139‑1 (MQ=255)                                                                                                            tAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCgaag                     >  2:1789320/1‑139 (MQ=255)                                                                                                               aaaaaTCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGaaa                  <  2:163355/139‑1 (MQ=255)                                                                                                               aaaaaTCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGaaa                  <  2:1612947/139‑1 (MQ=255)                                                                                                               aaaaaTCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGaaa                  <  2:1342597/139‑1 (MQ=255)                                                                                                                     ccGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTGGtt            <  2:183061/139‑1 (MQ=255)                                                                                                                        tgtTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGt                                                                                           >  1:80784/1‑57 (MQ=255)                                                                                                                        tgtTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGt                                                                                           <  2:80784/57‑1 (MQ=255)                                                                                                                         gtTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTGGTTATGa        <  1:1789320/139‑1 (MQ=255)                                                                                                                              cGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTGGTTATGACTGGt   <  1:1180986/139‑1 (MQ=255)                                                                                                                              cGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTGGTTATGACTGGt   <  2:1489750/139‑1 (MQ=255)                                                                                                                               gCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTGGTTATGACTGGTc  <  1:1749433/139‑1 (MQ=255)                                                                                                                               gCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTGGTTATGACTGGTc  <  2:974409/139‑1 (MQ=255)                                                                                                                                  gAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGt                                                                                           <  2:1518037/47‑1 (MQ=255)                                                                                                                                  gAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGt                                                                                           >  1:1518037/1‑47 (MQ=255)                                                                                                                                       |                                                                                                                                   GACCTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGTGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTGGTTATGACTGGTC  >  NC_000913/1862870‑1863134

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