Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,728,154	T→C	G106G (GGT→GGC)	gloA →	glyoxalase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,728,154	0	T	C	100.0%	218.1 / NA	60	G106G (GGT→GGC)	gloA	glyoxalase I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/28);  total (32/28)

GGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGTACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGCCGGGCGTGAACTCATCGCGCCCGCAT  >  NC_000913/1728018‑1728278                                                                                                                                         |                                                                                                                             ggggCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGcact                                                                                                                           >  2:369940/1‑140 (MQ=255)   ggCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGcactac                                                                                                                         >  2:139630/1‑140 (MQ=255)    gCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACg                                                                                                                        >  2:651255/1‑140 (MQ=255)       tGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGtt                                                                                                                     >  2:613112/1‑140 (MQ=255)       tGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGtt                                                                                                                     >  1:230960/1‑140 (MQ=255)         gATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTAt                                                                                                                   >  1:907254/1‑140 (MQ=255)              aTACGAACTCGTGACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAACCGTGCGAAAAAATCCGTCAAAACGGGGGTATCGTGACCTGTGAAGCTGGTCCGGTAAAGGGCGACACTACTGTTATCGCGt                                                                                                              <  2:1072908/140‑1 (MQ=255)              aTACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGt                                                                                                              >  1:510925/1‑140 (MQ=255)                     cTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGa                                                                                                       >  2:467564/1‑140 (MQ=255)                     cTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGa                                                                                                       <  1:248772/140‑1 (MQ=255)                           aCTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATcc                                                                                                 >  1:996002/1‑140 (MQ=255)                              gCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGa                                                                                              <  2:907254/140‑1 (MQ=255)                              gCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGa                                                                                              >  1:679852/1‑140 (MQ=255)                                    ggTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGGAAAAATCCGTCAAAACGGAGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTgg                                                                                                        >  2:1091970/1‑124 (MQ=255)                                    ggTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTgg                                                                                                        <  1:1091970/124‑1 (MQ=255)                                    ggTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTa                                                                                        <  1:651255/140‑1 (MQ=255)                                    ggTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTa                                                                                        >  2:184881/1‑140 (MQ=255)                                         caTCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCgg                                                                                               >  1:894984/1‑128 (MQ=255)                                         caTCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCgg                                                                                               <  2:894984/128‑1 (MQ=255)                                         caTCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGa                                                                                              <  1:383027/129‑1 (MQ=255)                                         caTCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGa                                                                                              >  2:383027/1‑129 (MQ=255)                                         caTCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACaaaa                                                                                   >  1:884811/1‑140 (MQ=255)                                         caTCGCGCTTAGCGTAGATAACGCCGCTCAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACaaaa                                                                                   >  1:1088646/1‑140 (MQ=255)                                           tCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAAtt                                                                                 <  2:170058/140‑1 (MQ=255)                                              cgcTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAg                                                                              >  2:1013498/1‑140 (MQ=255)                                              cgcTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAg                                                                              >  2:593087/1‑140 (MQ=255)                                                cTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAAtt                                                                                 <  2:38186/135‑1 (MQ=255)                                                cTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAAtt                                                                                 >  1:38186/1‑135 (MQ=255)                                                   aGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAAt                                                                         <  2:219467/140‑1 (MQ=255)                                                         gATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCg                                                                       <  2:1168890/136‑1 (MQ=255)                                                         gATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCg                                                                       >  1:1168890/1‑136 (MQ=255)                                                              cgccgcTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGaa                                                                     >  1:460781/1‑133 (MQ=255)                                                              cgccgcTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGaa                                                                     <  2:460781/133‑1 (MQ=255)                                                              cgccgcTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAg                                                              >  1:336271/1‑140 (MQ=255)                                                               gccgcTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGa                                                             >  1:866494/1‑140 (MQ=255)                                                               gccgcTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGa                                                             >  1:531163/1‑140 (MQ=255)                                                                ccgcTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGAc                                                            <  2:679852/140‑1 (MQ=255)                                                                 cgcTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACg                                                           >  1:1029664/1‑140 (MQ=255)                                                                 cgcTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACg                                                           <  1:369940/140‑1 (MQ=255)                                                                   cTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGcc                                                         >  1:1070517/1‑140 (MQ=255)                                                                    tGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCg                                                        >  2:3194/1‑140 (MQ=255)                                                                           tGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTAGAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCgg                                                 >  2:1026725/1‑140 (MQ=255)                                                                               aaaaaaTCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTg                                             <  2:884811/140‑1 (MQ=255)                                                                                 aaaaTCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTggg                                           <  2:936211/140‑1 (MQ=255)                                                                                                gggggTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGc                            <  2:336271/140‑1 (MQ=255)                                                                                                gggggTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGc                            <  1:184881/140‑1 (MQ=255)                                                                                                gggggTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGc                            <  1:467564/140‑1 (MQ=255)                                                                                                   ggTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGCCgg                         >  2:591944/1‑140 (MQ=255)                                                                                                   ggTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTCGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGCCgg                         >  1:190107/1‑140 (MQ=255)                                                                                                       aCGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGc                            <  2:150045/133‑1 (MQ=255)                                                                                                       aCGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGc                            >  1:150045/1‑133 (MQ=255)                                                                                                       aCGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGCCGGGCGt                     >  2:620166/1‑140 (MQ=255)                                                                                                        cGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGACAACTAATCTCCTGCCGGGCGTg                    <  1:1026725/140‑1 (MQ=255)                                                                                                                     gCGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGCCGGGCGTGAACTCATCGCGcc       <  1:593087/140‑1 (MQ=255)                                                                                                                      cGGGTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGCCGGGCGTGAACTCATCGCGccc      <  1:139630/140‑1 (MQ=255)                                                                                                                       gggTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGCCGGGCGTGAACTCATCGCGCCCg     <  1:613112/140‑1 (MQ=255)                                                                                                                       gggTCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGCCGGGCGTGAACTCATCGCGCCCg     <  2:1029664/140‑1 (MQ=255)                                                                                                                          tCCGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGCCGGGCGTGAACTCATCGCGCCCGCAt  <  2:531163/140‑1 (MQ=255)                                                                                                                            cGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGCCGGGCGTGAACTCATCGCGCCCg     >  1:1082156/1‑135 (MQ=255)                                                                                                                            cGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGCCGGGCGTGAACTCATCGCGCCCg     <  2:74887/135‑1 (MQ=255)                                                                                                                            cGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGCCGGGCGTGAACTCATCGCGCCCg     >  1:74887/1‑135 (MQ=255)                                                                                                                            cGGTAAAAGGCGGCACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGCCGGGCGTGAACTCATCGCGCCCg     <  2:1082156/135‑1 (MQ=255)                                                                                                                                         |                                                                                                                             GGGGCGTGGATAAATACGAACTCGGCACTGCTTATGGTCACATCGCGCTTAGCGTAGATAACGCCGCTGAAGCGTGCGAAAAAATCCGTCAAAACGGGGGTAACGTGACCCGTGAAGCGGGTCCGGTAAAAGGCGGTACTACGGTTATCGCGTTTGTGGAAGATCCGGACGGTTACAAAATTGAGTTAATCGAAGAGAAAGACGCCGGTCGCGGTCTGGGCAACTAATCTCCTGCCGGGCGTGAACTCATCGCGCCCGCAT  >  NC_000913/1728018‑1728278

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