Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	4,111,482	G→T	intergenic (+797/+115)	cdh → / ← yiiQ	CDP‑diacylglycerol phosphotidylhydrolase/DUF1454 family putative periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	4,111,482	0	G	T	100.0%	183.0 / NA	57	intergenic (+797/+115)	cdh/yiiQ	CDP‑diacylglycerol phosphotidylhydrolase/DUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (28/29);  total (28/29)

GCGAAGTCGATATAGGCAGTAGGGGGAGCACAAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGGGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATA  >  NC_000913_3_hsa_tpiA/4111335‑4111628                                                                                                                                                    |                                                                                                                                                   gCGAAGTGGTTATAGGTAGTGGGTTGGGCACAAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCTGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTgc                                                                                                                                                  <  2:517829/150‑1 (MQ=255)       tCGATATAGGCAGTAGGGGGAGCACAAACCCTGTCGGCTTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCAtt                                                                                                                                            >  1:763054/1‑150 (MQ=255)           tataGGCAGTAGGGGGAGCACAAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGTACTTCCTGGAGTGCGCCATTTTaa                                                                                                                                        <  1:886787/150‑1 (MQ=255)                gCAGTAGGGGGAGCACAAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATtctc                                                                                                                                   >  2:125371/1‑150 (MQ=255)                     aGGGGGAGCACAAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCAATTTTAATTCTCCACGc                                                                                                                              <  2:763054/150‑1 (MQ=255)                      gggggAGCACAAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCt                                                                                                                             >  1:236439/1‑150 (MQ=255)                      gggggAGCACAAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCt                                                                                                                             >  2:175182/1‑150 (MQ=255)                         ggAGCACAAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTtat                                                                                                                          <  2:639255/150‑1 (MQ=255)                            gCACAAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAg                                                                                                                       <  1:588913/150‑1 (MQ=255)                             cacaAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTTTaatt                                                                                                                      >  2:211154/1‑148 (MQ=255)                             cacaAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGCTCAGAGTGCCGATGAGCTCCCCAAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTGCCCCCAACGAAGAGATTCCCAGAGTGCGCCCTTTTCATTTCCCCCGCTTATaaac                                                                                                                      >  2:309114/1‑148 (MQ=255)                                aaaCCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGaa                                                                                                                   >  2:1005027/1‑150 (MQ=255)                                 aaCCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAAt                                                                                                                  >  1:257122/1‑150 (MQ=255)                                     accTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAg                                                                                                                       <  1:354834/141‑1 (MQ=255)                                     accTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAg                                                                                                                       >  2:354834/1‑141 (MQ=255)                                      cgtCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCACCGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCCCGCTTATAAGCGAATAAAgg                                                                                                             >  2:533219/3‑150 (MQ=255)                                       cTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATGTTAATTCTCCACGCTTATAAGCGAATAAAGGa                                                                                                            <  2:482467/150‑1 (MQ=255)                                          ggTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGa                                                                                                         <  1:509515/150‑1 (MQ=255)                                             gtCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATgg                                                                                                      >  1:42634/1‑150 (MQ=255)                                             gtCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATgg                                                                                                      <  1:211154/150‑1 (MQ=255)                                                            ggCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGgca                                                                                       >  2:80580/1‑150 (MQ=255)                                                               agcgTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCACACAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCCCCACGCTTATAAGCGAATAAGGGAAGTTGGCCGCCCTGCAGTgcagca                                                                                    >  2:56894/2‑150 (MQ=255)                                                                gcgTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGcagcag                                                                                   >  1:559443/1‑150 (MQ=255)                                                                gcgTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGACCCGCAGGGcagcag                                                                                   <  1:80580/150‑1 (MQ=255)                                                                   ttCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTc                                                                                >  1:644739/1‑150 (MQ=255)                                                                      agagTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCtgt                                                                             <  2:599317/150‑1 (MQ=255)                                                                          tGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGaaa                                                                         >  1:895890/1‑150 (MQ=255)                                                                             cGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAg                                                                                           <  2:1016380/129‑1 (MQ=255)                                                                             cGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAg                                                                                           >  1:1016380/1‑129 (MQ=255)                                                                             cGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAg                                                                      <  2:40211/150‑1 (MQ=255)                                                                                   gCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATaga                                                                >  1:835528/1‑150 (MQ=255)                                                                                    cTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATagag                                                               <  2:236439/150‑1 (MQ=255)                                                                                     tCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGa                                                                           <  2:430927/137‑1 (MQ=255)                                                                                     tCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGa                                                                           >  1:430927/1‑137 (MQ=255)                                                                                      cccccAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGAt                                                             <  1:829403/150‑1 (MQ=255)                                                                                        cccAGACTCTGCTTCAGCCCGGTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTACAGAGATTc                                                           >  2:382806/1‑150 (MQ=255)                                                                                            gACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCCGCAGGTCTGTGAAACAGTATAGAGATTCATCg                                                       <  2:224960/150‑1 (MQ=255)                                                                                            gACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCg                                                       <  2:252463/150‑1 (MQ=255)                                                                                                    ttCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCAGCACAAAg                                               >  2:208334/1‑150 (MQ=255)                                                                                                            gTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAg                                                                                           >  2:277831/1‑98 (MQ=255)                                                                                                            gTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAg                                                                                           <  1:277831/98‑1 (MQ=255)                                                                                                               cATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCtg                                                                              >  1:205422/1‑108 (MQ=255)                                                                                                               cATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCtg                                                                              <  2:205422/108‑1 (MQ=255)                                                                                                                aTCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCttttt                                   <  1:42564/150‑1 (MQ=255)                                                                                                                 tCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTg                                  <  2:559443/150‑1 (MQ=255)                                                                                                                             ccccAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCaa                      >  2:976579/1‑150 (MQ=255)                                                                                                                              cccAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTGTGCTTTTTGTCATTTATTCaaa                     >  1:473559/1‑150 (MQ=255)                                                                                                                               ccAACGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAAc                    >  2:29119/1‑150 (MQ=255)                                                                                                                                 aaCGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGTCAGCAGCTCTGTGAAACAGTATAGAGCTTCCTCGGCCCACAGGCTTTGCTTTTTGTCATATATTCAACCCt                  >  2:769845/1‑150 (MQ=255)                                                                                                                                  aCGAAGAACTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCtt                 <  1:29119/150‑1 (MQ=255)                                                                                                                                     aagaaCTTCCTGGAGTGCGCCATTTTAATTCTCCCCGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCaa              <  2:404485/150‑1 (MQ=255)                                                                                                                                     aagaaCTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCAGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCaa              <  1:208334/150‑1 (MQ=255)                                                                                                                                      agaaCTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAg             <  1:533219/150‑1 (MQ=255)                                                                                                                                       gaaCTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCaa              >  1:629265/1‑148 (MQ=255)                                                                                                                                       gaaCTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCaa              <  2:629265/148‑1 (MQ=255)                                                                                                                                         aCTTCCTGGAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGa          <  2:307237/150‑1 (MQ=255)                                                                                                                                                 gAGTGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATa  <  2:835528/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCGAAGTCGATATAGGCAGTAGGGGGAGCACAAACCACCTCGGTGTCGGCCGGCACCTTGGCCGCGTTCAGAGTGCCGATGAGCTCCCCCAGACTCTGCTTCCGCCCGTTCATCTTCCAGTTTCCCCCAACGAAGAACTTCCTGGAGGGCGCCATTTTAATTCTCCACGCTTATAAGCGAATAAAGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATA  >  NC_000913_3_hsa_tpiA/4111335‑4111628

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