Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,264,930	A→G	E406G (GAA→GGA)	hemA →	glutamyl‑tRNA reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,264,930	0	A	G	100.0%	218.8 / NA	59	E406G (GAA→GGA)	hemA	glutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/29);  total (30/29)

GCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGA  >  NC_000913/1264794‑1265064                                                                                                                                         |                                                                                                                                       gCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATaacggacg                                                                                                                                     >  2:1142228/1‑140 (MQ=255) gCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATaacggacg                                                                                                                                     <  1:912200/140‑1 (MQ=255) gCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATaacggacg                                                                                                                                     <  2:985864/140‑1 (MQ=255)  ccAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGc                                                                                                                                    <  1:927798/140‑1 (MQ=255)              cGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCGTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTg                                                                                                                        <  2:158776/140‑1 (MQ=255)                   cTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGa                                                                                                                   <  2:652045/140‑1 (MQ=255)                   cTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGa                                                                                                                   <  1:676129/140‑1 (MQ=255)                   cTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGa                                                                                                                   >  2:676129/1‑140 (MQ=255)                     tGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTg                                                                                                                                 <  2:346946/124‑1 (MQ=255)                     tGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTg                                                                                                                                 >  1:346946/1‑124 (MQ=255)                         cAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCt                                                                                                             <  1:1407484/140‑1 (MQ=255)                         cAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCt                                                                                                             <  1:1142228/140‑1 (MQ=255)                         cAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCt                                                                                                             <  2:379827/140‑1 (MQ=255)                            ggcggcGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCgg                                                                                                          >  2:1395028/1‑140 (MQ=255)                            ggcggcGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCgg                                                                                                          >  1:888746/1‑140 (MQ=255)                             gcggcgACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCggg                                                                                                         >  2:55357/1‑140 (MQ=255)                              cggcgACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGc                                                                                                        <  2:229067/140‑1 (MQ=255)                              cggcgACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGc                                                                                                        <  1:132196/140‑1 (MQ=255)                                 cgACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTgg                                                                                                     >  2:1380817/1‑140 (MQ=255)                                  gACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCGCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGa                                                                                                    <  1:373741/140‑1 (MQ=255)                                   aCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAg                                                                                                   <  1:55357/140‑1 (MQ=255)                                      cgcAAGCCATTATGCAGGATCTGGTATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAg                                                                                                <  1:72921/140‑1 (MQ=255)                                      cgcAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAg                                                                                                <  2:77287/140‑1 (MQ=255)                                      cgcAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAg                                                                                                >  2:484263/1‑140 (MQ=255)                                      cgcAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAg                                                                                                <  1:1380817/140‑1 (MQ=255)                                               ttATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGGAGCGGAAcatc                                                                                       >  2:1246972/1‑140 (MQ=255)                                                     aGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTc                                                                                 >  1:957126/1‑140 (MQ=255)                                                     aGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTc                                                                                 <  2:836686/140‑1 (MQ=255)                                                       gATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCtt                                                                               >  2:1075203/1‑140 (MQ=255)                                                         tCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCtttt                                                                             >  2:262529/1‑140 (MQ=255)                                                         tCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCtttt                                                                             >  2:1220925/1‑140 (MQ=255)                                                            ggCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCttttttt                                                                          >  2:233247/1‑140 (MQ=255)                                                             gCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTa                                                                         <  2:888746/140‑1 (MQ=255)                                                               aTGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACa                                                                       <  1:233247/140‑1 (MQ=255)                                                               aTGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACa                                                                       <  1:1075203/140‑1 (MQ=255)                                                               aTGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACa                                                                       <  2:947434/140‑1 (MQ=255)                                                                   aaaCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGt                                                                   >  1:1009156/1‑140 (MQ=255)                                                                   aaaCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTAAATCATTTTCTTTTTTTa                                                                         <  1:329337/134‑1 (MQ=255)                                                                   aaaCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTAAATCATTTTCTTTTTTTa                                                                         >  2:329337/1‑134 (MQ=255)                                                                       tGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCAt                                                               >  1:284490/1‑140 (MQ=255)                                                                         aCTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCAttt                                                             >  2:216915/1‑140 (MQ=255)                                                                           tAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTAc                                                           >  1:489609/1‑140 (MQ=255)                                                                            aaCCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACg                                                          <  1:484263/140‑1 (MQ=255)                                                                              ccGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGcc                                                        >  1:945755/1‑140 (MQ=255)                                                                                  ttGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATg                                                    <  1:477733/140‑1 (MQ=255)                                                                                   tGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGa                                                   <  2:530517/140‑1 (MQ=255)                                                                                    gATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGaa                                                  <  2:957126/140‑1 (MQ=255)                                                                                    gATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGaa                                                  <  1:1220925/140‑1 (MQ=255)                                                                                    gATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGaa                                                  <  1:1094670/140‑1 (MQ=255)                                                                                         aTGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCtt                                             >  2:516720/1‑140 (MQ=255)                                                                                          tGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTc                                            >  1:988062/1‑140 (MQ=255)                                                                                               cAACGAAATCACTTCAACAGGCCGCCCTTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCg                                       >  2:1039750/1‑140 (MQ=255)                                                                                               cAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCg                                       >  1:1054021/1‑140 (MQ=255)                                                                                                aaCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGt                                      >  1:968992/1‑140 (MQ=255)                                                                                                         aCTTCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACt                             >  1:979262/1‑140 (MQ=255)                                                                                                           ttCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTgg                           <  2:933758/140‑1 (MQ=255)                                                                                                            tCAACAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGa                          <  2:988062/140‑1 (MQ=255)                                                                                                              aaCAGGCAACCCGTGACGGGGATACCGGACGCCTCAATATTCTTCGCCACAACCTCGGGCTGGAGTAGCAGTACCTCATTTTCTTTTTTTACAGGGTGCATTGACGCCTATGAAGCCTTCTATCGTTGCCAAGCTGGAAg                        <  2:1305906/140‑1 (MQ=255)                                                                                                                cAGGCCGCCCGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGcc                      >  2:1185126/1‑140 (MQ=255)                                                                                                                        ccGTGACGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGa              >  1:45624/1‑140 (MQ=255)                                                                                                                             aCGGGGATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCC‑TCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCa        >  2:933824/1‑140 (MQ=255)                                                                                                                                 ggATAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATga     >  1:154878/1‑140 (MQ=255)                                                                                                                                    tAACGGACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATgaaga  >  1:321700/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                       GCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGA  >  NC_000913/1264794‑1265064

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