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%	224.2 / NA	63	R78R (CGT→CGC)	gapA	glyceraldehyde‑3‑phosphate dehydrogenase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/33);  total (30/33)

CTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGTGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTGGTTATGACTGGT  >  NC_000913/1862873‑1863133                                                                                                                                    |                                                                                                                                  cTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCTACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCgg                                                                                                                            <  1:745227/139‑1 (MQ=255) cTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCTACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCgg                                                                                                                            <  1:747408/139‑1 (MQ=255)  tGTTAGACACTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGc                                                                                                                           <  2:273257/139‑1 (MQ=255)      aGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAAc                                                                                                                       <  2:521318/139‑1 (MQ=255)          gCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGa                                                                                                                   >  2:528551/1‑139 (MQ=255)          gCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGa                                                                                                                   >  2:1205115/1‑139 (MQ=255)            tGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTCCCGCTGAACGCGATCCGGCTAACCTGaaa                                                                                                                 <  1:703528/139‑1 (MQ=255)                     ggCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGaa                                                                                                        >  2:558282/1‑139 (MQ=255)                        aTACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGtt                                                                                                     >  2:1113006/1‑139 (MQ=255)                          aCATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTgg                                                                                                   >  2:1070010/1‑139 (MQ=255)                            aTGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGgtg                                                                                                 >  2:181121/1‑139 (MQ=255)                              gCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGtt                                                                                               >  1:507971/1‑139 (MQ=255)                               cTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCgg                                                                                                                            <  1:496723/109‑1 (MQ=255)                               cTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCgg                                                                                                                            >  2:496723/1‑109 (MQ=255)                               cTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTg                                                                                              >  2:658930/1‑139 (MQ=255)                                tGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCt                                                                                                                          >  1:277932/1‑110 (MQ=255)                                tGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCt                                                                                                                          <  2:277932/110‑1 (MQ=255)                                   aaTATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGt                                                                                          <  1:27412/139‑1 (MQ=255)                                            ccACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGa                                                                                 >  1:863738/1‑139 (MQ=255)                                              aCTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAg                                                                               <  1:1155076/139‑1 (MQ=255)                                              aCTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAg                                                                               >  1:1061162/1‑139 (MQ=255)                                               cTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGc                                                                              <  2:372988/139‑1 (MQ=255)                                                     gCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGa                                                                                                                                 <  2:1003115/82‑1 (MQ=255)                                                      ccGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                                                       <  1:782916/139‑1 (MQ=255)                                                         tttCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCgg                                                                                                                            >  2:912830/1‑83 (MQ=255)                                                         tttCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCgg                                                                                                                            <  1:912830/83‑1 (MQ=255)                                                         tttCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGa                                                                                                            <  1:271182/99‑1 (MQ=255)                                                         tttCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGa                                                                                                            >  2:271182/1‑99 (MQ=255)                                                         tttCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTg                                                                    >  1:1157048/1‑139 (MQ=255)                                                           tCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGtt                                                                  <  1:1113006/139‑1 (MQ=255)                                                               cGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCctg                                                              <  2:863738/139‑1 (MQ=255)                                                                ggTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCctga                                                             <  2:507971/139‑1 (MQ=255)                                                                    ccGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCctgactga                                                         >  1:685934/1‑139 (MQ=255)                                                                          aaGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAAc                                                   <  1:558282/139‑1 (MQ=255)                                                                              gAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCt                                               >  2:1170351/1‑139 (MQ=255)                                                                                 aGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCaa                                                                            <  1:1158890/107‑1 (MQ=255)                                                                                 aGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCaa                                                                            >  2:1158890/1‑107 (MQ=255)                                                                                     ggTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGTCTGACGAAACTGCTCGTAAac                                        >  1:1113800/1‑139 (MQ=255)                                                                                             gATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCAcc                                >  2:451759/1‑139 (MQ=255)                                                                                              aTCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCg                               <  1:1070010/139‑1 (MQ=255)                                                                                                  ttAACGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTgg                           <  1:630487/139‑1 (MQ=255)                                                                                                     aCGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCgg                                                                                                                            <  2:1060478/39‑1 (MQ=255)                                                                                                     aCGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCgg                                                                                                                            >  1:1060478/1‑39 (MQ=255)                                                                                                     aCGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGc                                                                                                                           >  1:518230/1‑40 (MQ=255)                                                                                                     aCGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGc                                                                                                                           <  2:518230/40‑1 (MQ=255)                                                                                                     aCGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGa                                                                                                            >  2:758079/1‑55 (MQ=255)                                                                                                     aCGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGa                                                                                                            <  1:758079/55‑1 (MQ=255)                                                                                                     aCGGTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGc                        >  2:665241/1‑139 (MQ=255)                                                                                                        gTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGaa                                                                                                        <  2:163974/56‑1 (MQ=255)                                                                                                        gTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGaa                                                                                                        >  1:163974/1‑56 (MQ=255)                                                                                                        gTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCgaa                     <  1:665241/139‑1 (MQ=255)                                                                                                        gTAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCgaa                     >  1:620184/1‑139 (MQ=255)                                                                                                            aaaaaTCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGa                                                                                 <  1:519642/75‑1 (MQ=255)                                                                                                            aaaaaTCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGaaa                 <  1:439155/139‑1 (MQ=255)                                                                                                            aaaaaTCCGTGTTACCGCTGAACGCGATCCGGCTAACCCGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGa                                                                                 >  2:519642/1‑75 (MQ=255)                                                                                                               aaTCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTg              <  2:685934/139‑1 (MQ=255)                                                                                                                 tCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTGGt            >  1:276636/1‑139 (MQ=255)                                                                                                                   cGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTGGTTa          >  1:519279/1‑139 (MQ=255)                                                                                                                          ccGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCAAAGAAAGTGGTTATGACTgg   <  2:1157048/139‑1 (MQ=255)                                                                                                                           cGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTGGTTATGACTGGt  <  1:1205115/139‑1 (MQ=255)                                                                                                                           cGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTGGTTATGACTGGt  <  1:528551/139‑1 (MQ=255)                                                                                                                             cTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCctgac                                                            <  2:509140/79‑1 (MQ=255)                                                                                                                             cTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCctgac                                                            >  1:509140/1‑79 (MQ=255)                                                                                                                                    |                                                                                                                                  CTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGTGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGTGCGAAGAAAGTGGTTATGACTGGT  >  NC_000913/1862873‑1863133

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