Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,862,976	G→A	G69D (GGT→GAT)	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,862,976	0	G	A	100.0%	230.2 / NA	70	G69D (GGT→GAT)	gapA	glyceraldehyde‑3‑phosphate dehydrogenase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/35);  total (35/35)

TTCTGACATCGAGATCGTTGCAATCAACGACCTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGTGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGT  >  NC_000913/1862842‑1863109                                                                                                                                       |                                                                                                                                      ttCTGACATCGAGATCGTTGCAATCAACGACCTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATaaa                                                                                                                                   <  2:1332574/139‑1 (MQ=255)       cATCGAGATCGTTGCAATCAACGACCTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATc                                                                                                                             <  1:796565/139‑1 (MQ=255)              aTCGTTGCAATCAACGACCTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTa                                                                                                                      <  1:1577438/139‑1 (MQ=255)                  ttGCAATCAACGACCTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGc                                                                                                                  >  1:1296204/1‑139 (MQ=255)                   tGCAATCAACGACCTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCt                                                                                                                 <  2:1077538/139‑1 (MQ=255)                          aaCGACCTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTAATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACgcg                                                                                                          >  2:1120662/1‑137 (MQ=255)                                cTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCgg                                                                                                    <  1:1147216/139‑1 (MQ=255)                                cTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCgg                                                                                                    <  2:1296204/139‑1 (MQ=255)                                   ttAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGc                                                                                                   <  2:304574/137‑1 (MQ=255)                                   ttAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGc                                                                                                   >  1:304574/1‑137 (MQ=255)                                             aTTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCgtgt                                                                                                                        <  2:629075/106‑1 (MQ=255)                                             aTTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCgtgt                                                                                                                        >  1:629075/1‑106 (MQ=255)                                             aTTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATg                                                                                       >  2:28284/1‑139 (MQ=255)                                                aCATGGCATACATGCTGAAATTTGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGTTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGaaacgcga                                                                                    <  2:962466/139‑6 (MQ=255)                                                   tGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGa                                                                                    <  2:954825/136‑1 (MQ=255)                                                   tGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGa                                                                                    >  1:954825/1‑136 (MQ=255)                                                       aTACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCgg                                                                                                    <  2:1338413/116‑1 (MQ=255)                                                       aTACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCgg                                                                                                    >  1:1338413/1‑116 (MQ=255)                                                       aTACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGtt                                                                             >  1:1236307/1‑139 (MQ=255)                                                           aTGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGgtg                                                                         >  1:372839/1‑139 (MQ=255)                                                                gAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGAc                                                                    >  1:1524016/1‑139 (MQ=255)                                                                         cTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCg                                                             >  1:1247720/1‑137 (MQ=255)                                                                         cTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCg                                                             <  2:1247720/137‑1 (MQ=255)                                                                           ccACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGa                                                         >  1:1551732/1‑139 (MQ=255)                                                                           ccACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGa                                                         >  2:229204/1‑139 (MQ=255)                                                                              cTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGc                                                      >  2:1258936/1‑139 (MQ=255)                                                                              cTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGc                                                      >  2:229806/1‑139 (MQ=255)                                                                              cTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGc                                                      >  1:1470032/1‑139 (MQ=255)                                                                                cACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCaa                                                    <  1:229806/139‑1 (MQ=255)                                                                                cACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCaa                                                    <  2:1570568/139‑1 (MQ=255)                                                                                 aCGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAAc                                                   >  1:1593630/1‑139 (MQ=255)                                                                                   ggCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTa                                                                                                                      <  1:1417574/70‑1 (MQ=255)                                                                                   ggCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTa                                                                                                                      >  2:1417574/1‑70 (MQ=255)                                                                                     ccGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                               <  1:11999/139‑1 (MQ=255)                                                                                     ccGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                               <  1:28284/139‑1 (MQ=255)                                                                                     ccGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                               <  2:1524016/139‑1 (MQ=255)                                                                                        tttCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTa                                                                                                 <  1:373207/86‑1 (MQ=255)                                                                                        tttCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTa                                                                                                 >  2:373207/1‑86 (MQ=255)                                                                                          tCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGtt                                          <  2:1130712/139‑1 (MQ=255)                                                                                                   ccGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                               >  2:1515528/1‑125 (MQ=255)                                                                                                   ccGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                               <  1:1515528/125‑1 (MQ=255)                                                                                                   ccGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCctgactga                                 >  1:1413145/1‑139 (MQ=255)                                                                                                        gAAGTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGaaa                            <  2:793833/139‑1 (MQ=255)                                                                                                           gTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                               >  2:228253/1‑117 (MQ=255)                                                                                                           gTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                               <  1:228253/117‑1 (MQ=255)                                                                                                           gTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTg                         <  2:372839/139‑1 (MQ=255)                                                                                                           gTGAAAGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTg                         <  1:1120662/139‑1 (MQ=255)                                                                                                                aGACGGTCATCTGATCGTTTACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGa                              <  2:659462/129‑1 (MQ=255)                                                                                                                aGACGGTCATCTGATCGTTTACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGa                              >  1:659462/1‑129 (MQ=255)                                                                                                                aGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                               >  2:317832/1‑112 (MQ=255)                                                                                                                aGACGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                               <  1:317832/112‑1 (MQ=255)                                                                                                                   cGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGtt                                                                 <  2:350811/91‑1 (MQ=255)                                                                                                                   cGGTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGtt                                                                 >  1:350811/1‑91 (MQ=255)                                                                                                                   cGGCCATCTGATCGTTAACGATAAAGAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTaaa                 <  1:1068245/139‑1 (MQ=255)                                                                                                                    ggTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCaa                                                    <  2:532812/103‑1 (MQ=255)                                                                                                                    ggTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCaa                                                    >  1:532812/1‑103 (MQ=255)                                                                                                                     gTCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCAGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAaca               >  1:449391/1‑139 (MQ=255)                                                                                                                      tCATCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAacac              >  2:1319745/1‑139 (MQ=255)                                                                                                                        aTCTGATCGTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACAt            <  2:1236307/139‑1 (MQ=255)                                                                                                                                gTTAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTg    >  2:1098773/1‑139 (MQ=255)                                                                                                                                  tAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGa                                                                                    <  2:311293/57‑1 (MQ=255)                                                                                                                                  tAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGa                                                                                    >  1:311293/1‑57 (MQ=255)                                                                                                                                  tAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                               >  2:637066/1‑94 (MQ=255)                                                                                                                                  tAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                               <  1:637066/94‑1 (MQ=255)                                                                                                                                  tAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGt  >  1:345702/1‑139 (MQ=255)                                                                                                                                  tAACGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGt  >  2:766982/1‑139 (MQ=255)                                                                                                                                   aaCGATAAAAAAATCCGTGTTACCGCTGAACgcga                                                                                                         >  2:1315663/1‑32 (MQ=255)                                                                                                                                   aaCGATAAAAAAATCCGTGTTACCGCTGAACgcga                                                                                                         <  1:1315663/35‑4 (MQ=255)                                                                                                                                    aCGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                               <  1:1198944/92‑1 (MQ=255)                                                                                                                                    aCGATAAAAAAATCCGTGTTACCGCTGAACGCGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGt                                               >  2:1198944/1‑92 (MQ=255)                                                                                                                                       |                                                                                                                                      TTCTGACATCGAGATCGTTGCAATCAACGACCTGTTAGACGCTGATTACATGGCATACATGCTGAAATATGACTCCACTCACGGCCGTTTCGACGGTACCGTTGAAGTGAAAGACGGTCATCTGATCGTTAACGGTAAAAAAATCCGTGTTACCGCTGAACGTGATCCGGCTAACCTGAAATGGGACGAAGTTGGTGTTGACGTTGTCGCTGAAGCAACTGGTCTGTTCCTGACTGACGAAACTGCTCGTAAACACATCACCGCTGGT  >  NC_000913/1862842‑1863109

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