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%	163.2 / NA	56	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (28/27);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATTT  >  NC_000913_3_hsa_tpiA/263952‑264243                                                                                                                                                   |                                                                                                                                                   cGGGCTCCCATGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGTTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAg                                                                                                                                                  <  1:369901/149‑1 (MQ=255)       cccAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTAGCTCATTACAGTTCTCCCCTATTGTCGGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCACTTGCGTTCCTGGCGGCGTTTGGCGTATACCGCTAGGACGcc                                                                                                                                            >  1:262849/1‑149 (MQ=255)        ccAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCg                                                                                                                                           <  1:4753/149‑1 (MQ=255)                    tAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGa                                                                                                                               >  2:511905/1‑149 (MQ=255)                         tcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGTACTCCGGCGATATCCGGAACGTc                                                                                                                          >  1:95472/1‑149 (MQ=255)                                gCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTCTTGTCAGATCCGTTATACGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCg                                                                                                                   <  1:718101/149‑1 (MQ=255)                                  aGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGGACGTCCCTGGCGgt                                                                                                                 >  1:133849/1‑149 (MQ=255)                                     ccccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTCCCGGTGTa                                                                                                              >  2:498131/1‑149 (MQ=255)                                       acTGGGAGGGTTGCTCATTTCAGTGCTCCGCTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                            <  1:544495/148‑1 (MQ=255)                                        cTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATg                                                                                                           >  2:552069/1‑149 (MQ=255)                                         tGGTAGGTTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGa                                                                                                          <  1:63028/149‑1 (MQ=255)                                         tGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGa                                                                                                          >  1:323002/1‑149 (MQ=255)                                          ggTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGAc                                                                                                         <  2:520703/149‑1 (MQ=255)                                             aGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGa                                                                                                      <  2:323002/149‑1 (MQ=255)                                                gTTCCTCATTACAGTTCTCCACTATTGTCAGTTCAGGTTGATGCTCATCAGGCCGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGCACGTCCCTGCCGGTGTAATGACAGACGc                                                                                                   <  1:247605/149‑1 (MQ=25)                                                  tGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATTCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATcc                                                                                                                                  >  1:95362/1‑116 (MQ=255)                                                  tGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATcc                                                                                                                                  <  2:95362/116‑1 (MQ=255)                                                    cTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTAcc                                                                                               <  1:552069/149‑1 (MQ=255)                                                      cATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGt                                                                                             <  1:615184/149‑1 (MQ=255)                                                      cATTACAGTCTTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGt                                                                                             <  2:495342/149‑1 (MQ=255)                                                       aTTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGGTAGGACTCCGGCGATATCCTGAACGACCCTGCCGGTGTAATGACAGACGCTACCGTg                                                                                            >  2:430576/1‑149 (MQ=255)                                                                ccccACTAGTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTa                                                                                   <  1:869823/147‑1 (MQ=255)                                                                 tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAAGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAt                                                                                  >  1:763879/1‑149 (MQ=255)                                                                 tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAACAGATGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGGAt                                                                                  >  2:955270/1‑149 (MQ=18)                                                                     cgaTTGTCAGGTCGGGGGGATGCTCGTTAGTCAGGTAGAGGGCCCATGGTGGTCCGGGCGGCGTGCGGCGTATACCGCGAGGACGGCGGCGATATCCGGAACGGCCCTGCAGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑Gc                                                                             <  2:218621/147‑1 (MQ=21)                                                                     cTATTGTCATTTCAGGGTGATGCTCATCAGCCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑Gc                                                                             <  1:498657/149‑1 (MQ=255)                                                                     cTATTGTCAGTTGAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTACTCGGGTCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑Gc                                                                             <  1:761528/149‑1 (MQ=21)                                                                       aTTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGTCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTTGTAGTT‑GCCg                                                                           <  2:903115/149‑1 (MQ=255)                                                                        ttGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCgg                                                                          >  1:224246/1‑149 (MQ=255)                                                                           tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGc                                                                       >  1:206022/1‑149 (MQ=255)                                                                            ccgtgcAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCa                                                                      <  2:393818/144‑1 (MQ=255)                                                                                       gATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTCATCAAAGACTCTACCGTGCCACTGGTATTT‑GCCGCTGCAGTAGCGAAGGa                                                           >  1:970619/1‑149 (MQ=255)                                                                                          gCTAATCAGGCAGGTTTAGGGCCCATTGCGCTCCTGGCGGCGTTCGGCGTATACCGCTAGGAGTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTc                                                        <  2:305218/149‑1 (MQ=21)                                                                                            tcatcaGGCAGGTATAGGGCCCATTGCGGGCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTCgg                                                      <  2:177371/149‑1 (MQ=255)                                                                                              atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAACGATTTGGGa                                                    >  1:801866/1‑149 (MQ=255)                                                                                              atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTCgggc                                                    >  1:777599/1‑148 (MQ=255)                                                                                              atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTCGGGa                                                    >  1:866343/1‑149 (MQ=255)                                                                                              atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAAGCGCTACCGTGCCCCTCGTATTT‑GCCGGTGCAGTAGCGAAAGATTCGGGa                                                    >  1:659663/1‑149 (MQ=21)                                                                                              atcaGGCAGCTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACGCCGGCGATATCCGGAACGTCCCCGCCGGTGTAATGACAAACGCTACCGTGCCACTGCCGTGT‑GCCGGTGCAGTCGCGAAACATTCcggc                                                    >  2:576049/1‑148 (MQ=255)                                                                                                       gTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATg                                           >  1:838549/1‑149 (MQ=255)                                                                                                       gTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATg                                           >  1:294173/1‑149 (MQ=255)                                                                                                            gggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTCGTCTG‑GCAGTAGCGAAAGATTTGGGACTCAGGATGCTGGc                                      <  1:511905/149‑1 (MQ=18)                                                                                                             ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCg                                     <  1:254325/149‑1 (MQ=255)                                                                                                             ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCg                                     <  1:253372/149‑1 (MQ=255)                                                                                                             ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCg                                     <  1:850563/149‑1 (MQ=255)                                                                                                                       ggTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTc                           >  1:720950/1‑149 (MQ=255)                                                                                                                       ggTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTA‑GCTGGTTCAGTAGCGAAAGATTCGGGCCCCAGGATGCTGGCGGTGTATCgcc                           >  1:197448/1‑147 (MQ=17)                                                                                                                             ggcggcGTTCTGCGTATACCGCTAGGACTCCGGCGATATCCGGTACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTACTT‑GCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGcc                     <  2:45181/149‑1 (MQ=255)                                                                                                                               cggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTGCCGTGCCACGCGTATTT‑GCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCt                   >  2:378194/1‑149 (MQ=255)                                                                                                                                   gTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGt               >  1:607522/1‑149 (MQ=255)                                                                                                                                   gTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGAAGTAGAGAAAGATTGGGGACTCAGGATGCTGGCGGTATATCGTCATTGCGCTGCGt               <  1:824903/149‑1 (MQ=17)                                                                                                                                      ccgCGTATACCGCTAGGACTCCTGCGACATCCGGAACGTCCCTGCCGGCGTAATGCCAGACGCTACCGGCCCACTGGTATTT‑GCCGGTGCATTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTa            <  1:120718/147‑1 (MQ=14)                                                                                                                                       ggCGTATCCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGTTTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTAc           <  1:560087/149‑1 (MQ=255)                                                                                                                                        gCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAAACGCTACCGTGCCACTGGGATTC‑GCCGGTGCAGTAGCGAAGGATTCGGGACTGAGGAAGCAGGCGGGAGATCGTCAGTGCCCGGCggtcccg          >  2:662219/1‑143 (MQ=14)                                                                                                                                             taCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATaa     >  2:621392/1‑149 (MQ=255)                                                                                                                                               ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAtt   >  1:533493/1‑149 (MQ=255)                                                                                                                                                cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAACGATTAGGGACTCAGGATGCTGGCGCGATATCGTCATTGCCCTTCGTTTACTGCTAAttt  >  2:785692/1‑149 (MQ=14)                                                                                                                                                   |                                                                                                                                                   CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTT‑GCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATTT  >  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 < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑