Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	264,098	G→T	L50L (CTC→CTA)	ykfH ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	264,098	0	G	T	98.2%	169.4 / NA	57	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (32/23);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATTT  >  NC_000913_3_hsa_tpiA/263952‑264243                                                                                                                                                   |                                                                                                                                                   cGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGCGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAg                                                                                                                                                  >  1:1041069/1‑149 (MQ=255)     cTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACt                                                                                                                                              >  2:977716/1‑149 (MQ=255)        ctaGGCGCAGTGGAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCg                                                                                                                                           <  2:165190/147‑1 (MQ=255)                      aTCTCCCGCTGCAGGCCCCGGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAAc                                                                                                                             >  1:761879/1‑149 (MQ=255)                         tcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTc                                                                                                                          <  2:1070935/149‑1 (MQ=255)                           ccGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTccc                                                                                                                        >  1:583922/1‑149 (MQ=255)                              cTGCAGGTCGATGGTAGGGTGGCTCATTACCGGTCTCCACTGTGGTCAGTGCAGGGTGATGCTTCTTACGCAGCTACCGGCCCCATCGCGGTCCTGGCGGCGTTCGGCGTACACCGCCAGCACTCCGGCGATATCCGGAACGTCCCTGc                                                                                                                     <  2:597093/149‑1 (MQ=255)                                  aGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgt                                                                                                                 >  1:875363/1‑149 (MQ=255)                                  aGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgt                                                                                                                 >  2:289223/1‑149 (MQ=255)                                  aGGCCCCTGGTAGGGGTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCAGAACGTCCCTGCCGgt                                                                                                                 <  1:435082/149‑1 (MQ=255)                                     ccccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTa                                                                                                              <  1:949654/149‑1 (MQ=255)                                         gggTAGGGTTGCTCATTACAGTCTTCCACTATTGTCAGTTCCGGGTGATGCTCTTCAGCTTCTTATTGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGa                                                                                                          <  1:1145762/148‑1 (MQ=25)                                          ggTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATGGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGAc                                                                                                         <  2:716993/149‑1 (MQ=255)                                                    cTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCCTCAGGCAGGTCTAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGATAGGCCTCAGCCGATATCCGGAACATCCCTGACGGTGCAATGACGGCGGCTAcc                                                                                               >  1:1115562/1‑149 (MQ=25)                                                     tCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCg                                                                                              >  1:271791/1‑149 (MQ=255)                                                        ttACAGTTCTCCCCTATTGTCAGTCAAGGGTGATGCTCATCAGGCAGGTATTGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATCTCCGGAACGTCCCTCCCGGGGTCATGGCAGACGCTACCGTGc                                                                                           >  2:96484/1‑149 (MQ=255)                                                         tACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCCGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGcc                                                                                          >  2:483729/1‑149 (MQ=255)                                                            aGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCATTCGGCGTATACCGCTAGGACTCCGGCGATATCCGAAAAGTCCCCGCCGGAAAAATTACAGGCGCTACCGTGccccc                                                                                       >  2:883491/1‑146 (MQ=21)                                                             gTTCTCCACTTTTGTCTGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTg                                                                                      <  2:360270/149‑1 (MQ=255)                                                                 tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                            >  2:55613/1‑123 (MQ=255)                                                                 tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                            <  1:55613/123‑1 (MQ=255)                                                                     cTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGc                                                                             <  1:635741/149‑1 (MQ=255)                                                                      tATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTAGCGTGCCCCTGGT‑ATTTccc                                                                            >  1:1002402/1‑149 (MQ=255)                                                                       aTTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCg                                                                           <  2:742515/149‑1 (MQ=255)                                                                        ttGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTGACCgg                                                                          >  2:934512/1‑149 (MQ=255)                                                                          gTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTg                                                                        >  1:839752/1‑149 (MQ=255)                                                                          gTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGG‑ATTTGCCGGta                                                                        >  1:833646/1‑148 (MQ=255)                                                                           gcAGTTCAGGGAGATGCTCATCATGCAGGTATAGTGCCCAGTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTAACTGCAGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGc                                                                       >  2:544570/2‑149 (MQ=21)                                                                            cAGTGCAGGGTGATGCCCAGCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGCTATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCa                                                                      <  1:413932/149‑1 (MQ=21)                                                                               ttCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTa                                                                   <  1:844611/149‑1 (MQ=255)                                                                                       gATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGACGGTGTAATGACAGACGCTACCGTGCCACTGGC‑ATTTGCCGGTCCAGTAGCGAAAGa                                                           >  1:243001/1‑149 (MQ=255)                                                                                        aTGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGGGAAAGAt                                                          >  1:818691/1‑149 (MQ=255)                                                                                        aTGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGAt                                                          >  1:818702/1‑149 (MQ=255)                                                                                            tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCgg                                                      >  1:770013/1‑149 (MQ=255)                                                                                              atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGa                                                    <  2:1119805/149‑1 (MQ=255)                                                                                              atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGa                                                    <  2:761879/149‑1 (MQ=255)                                                                                                 aggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTc                                                 <  2:271791/149‑1 (MQ=255)                                                                                                 aggcaggTATAGGGCCAATTGCGGTCCTGGCGGCGAGCCACGCATACCCCTAGGACTCCGTCGATATCCGGAACGGCCCTGCCCGTGTAATGACAGACGCTACCGGGCCACTGGT‑ATTTCCCCGTCCAGTAGCGAAAGATTCtggcccc                                                 >  2:532653/1‑145 (MQ=18)                                                                                                  ggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCa                                                <  1:289223/149‑1 (MQ=255)                                                                                                   gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAg                                               <  2:270703/149‑1 (MQ=255)                                                                                                   gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAg                                               >  1:551063/1‑149 (MQ=255)                                                                                                   gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACCCAg                                               >  2:84094/1‑149 (MQ=255)                                                                                                             ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCg                                     >  1:998446/1‑149 (MQ=255)                                                                                                              gCCCGTTGCGGTCCGGGGGGCGGTGGGCGTATACCGCTAGGACTCCGGCGACTTCGGGAACGTCCCTGCTGGTGTAATGGCAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGGGAAAGAGTCTGTACTCAGGATGCTGGCgg                                    <  2:1107243/149‑1 (MQ=255)                                                                                                                ccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTTCCa                                                                                         <  1:443080/95‑1 (MQ=255)                                                                                                                ccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGta                                  >  2:629242/1‑149 (MQ=255)                                                                                                                ccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATACGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGta                                  >  1:336942/1‑149 (MQ=255)                                                                                                                ccAGTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTTCCa                                                                                         >  2:443080/1‑95 (MQ=18)                                                                                                                         tCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTG‑CACCGTTGAGGTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCAt                         >  1:1146801/1‑149 (MQ=17)                                                                                                                          ccggccggcgTGCGGCTTAGCCCGCGAGGACTCCGGCCCTTCCCGCAACGGCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCAtt                        <  1:806565/144‑1 (MQ=14)                                                                                                                               cggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCt                   <  2:202232/149‑1 (MQ=255)                                                                                                                               cggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCt                   >  1:202232/1‑149 (MQ=255)                                                                                                                                ggcgTTCGGCAACTAACGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTg                  >  2:438366/1‑149 (MQ=17)                                                                                                                                 ggggctggCGTATACGGATAGGACCGCGGCGAAATCCGGAACGTCCCTGCCGGTGGAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGc                 <  2:1002402/143‑1 (MQ=14)                                                                                                                                              aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAt    <  2:551063/149‑1 (MQ=255)                                                                                                                                               ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCg                                                                                              <  2:1101879/59‑1 (MQ=255)                                                                                                                                               ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCg                                                                                              >  1:1101879/1‑59 (MQ=255)                                                                                                                                                cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACCCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAttt  <  1:84094/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT‑ATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATTT  >  NC_000913_3_hsa_tpiA/263952‑264243

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