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.3%	162.6 / NA	58	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (25/32);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAT  >  NC_000913_3_hsa_tpiA/263952‑264241                                                                                                                                                   |                                                                                                                                                cGGGCGCCCAGGCGCAGTGTAATCGCCCGTTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATCCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAg                                                                                                                                               <  2:123430/149‑1 (MQ=39)   ggCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGa                                                                                                                                             <  2:85330/149‑1 (MQ=255)   cgCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCTTCGTTCGGCGTATACCGCTAGGa                                                                                                                                             <  1:590982/148‑1 (MQ=255)      tCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTCTTGTCAGTTCAGGGTCATGCTCGTCAGGCACGTATAGGGCCCATTGCGGCCCTGGAGGCGGTCGGCGTCTACCGCGAGGAccc                                                                                                                                          >  1:293744/1‑147 (MQ=255)      tCCCAGGCGCAGTGTAAGCTCCCGCGGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTc                                                                                                                                          <  2:155970/149‑1 (MQ=255)                 gtgtATTCTCTCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATcc                                                                                                                               <  1:360872/149‑1 (MQ=255)                  tgtAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCg                                                                                                                              >  2:586556/1‑149 (MQ=255)                  tgtAATCTCCCGCTGCAGGCCCCTGGTAGCGTTGCTCATTACAGTTCTCCACTATTGTCACTTCACGGTGATGCTCATCAGGCAGGTATAGGCCCCATTGCGGTCCTCGCGGCGTTCGGCGTATACCGCTAGGCCTCCAGCGATATCCg                                                                                                                              >  2:834326/1‑149 (MQ=255)                    tAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGa                                                                                                                            >  1:587029/1‑149 (MQ=255)                     aaTCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGCTATCCGGag                                                                                                                           >  1:166100/1‑148 (MQ=255)                      atctaacGCTGCAGGTCCCGTGTAGGGTTGCTCATTACAGTTCGCCGCTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAAc                                                                                                                          <  2:587029/143‑1 (MQ=255)                                    gCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAAGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtgt                                                                                                            >  1:496271/1‑149 (MQ=255)                                      cccAGGTAGGGGTGTTCATTACGATTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                          <  2:165180/149‑1 (MQ=255)                                        cTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATg                                                                                                        <  2:101665/149‑1 (MQ=255)                                        cTGGTAGGGTTGCTCATTACAGTTCTCCACTA‑TGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGa                                                                                                       >  2:220137/1‑149 (MQ=255)                                         tGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATGGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGa                                                                                                       <  1:723762/149‑1 (MQ=255)                                            tAGGGTTGGGAACTACAGTGCTCCACTATTGTCAGGTCAGGGTGATGCTCATCAGGAAGGTATAGGGGCCATTGCGGTCCTGGCGGCGTTCGGCGTAGACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAg                                                                                                    <  2:373252/149‑1 (MQ=21)                                             aGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGa                                                                                                   <  1:586556/149‑1 (MQ=255)                                                  tGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTACCGGTGTAATGACAGACGCTa                                                                                              >  2:752090/1‑149 (MQ=255)                                                       aTTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGGATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGAAACGTCCCTGCCGGTGTACTGGCAGACGCTACCGTg                                                                                         >  1:306162/1‑149 (MQ=25)                                                            aGTGCTCCACTATTTTCAGGTCAGGGTGATGCTCAGCAGGCAGGTATAGGGCCCGTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACt                                                                                    <  2:254504/149‑1 (MQ=21)                                                                cttcACTAGTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTa                                                                                <  2:562301/146‑1 (MQ=255)                                                                 tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAt                                                                               >  2:413098/1‑149 (MQ=255)                                                                 tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGGAt                                                                               >  2:784352/1‑149 (MQ=255)                                                                   cccgaTTGTCCGTTCCGGGTCATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTCCGGCGTATACCGCTTGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAttt                                                                             <  2:166100/145‑1 (MQ=21)                                                                   cACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTCCCACTGCTAttt                                                                             >  2:381189/1‑149 (MQ=255)                                                                        ttGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCgg                                                                        >  2:436996/1‑149 (MQ=255)                                                                          gTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCg                                                                         >  1:335587/1‑146 (MQ=255)                                                                          gTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCg                                                                         <  2:335587/146‑1 (MQ=255)                                                                           taaGTTCCGGGGGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCGTGGCGGCGTTGGGCGTATACCGCTAGGACTCGGGCGGTATCCGGAACGTGCCTGCCGGTGTAATGACAGTCGCTACCGTGCGAATGGTATTTGTCGGTGc                                                                     >  2:400060/3‑149 (MQ=255)                                                                           tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGc                                                                     >  1:395120/1‑149 (MQ=255)                                                                           tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCATGCCGCTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGc                                                                     >  2:765144/1‑149 (MQ=255)                                                                            aaGTTCAGGGGGGTGCTTATTAGGGTGGTATAGGGCGGATGGCAGTCGTGGCGGGGTTCGGCGTTTACCGGTAGGCCTCCGGCGATATCCGGATCGTCGCTGCCGGTGTTATGACAGACGCTACCGTGCCACTGGTAGTTGCCGGTGCa                                                                    <  2:783013/148‑1 (MQ=16)                                                                                tCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAg                                                                <  1:220137/149‑1 (MQ=255)                                                                                    ggTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTGGGCGTGTACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaa                                                            <  1:823490/149‑1 (MQ=255)                                                                                      tGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAg                                                          >  1:474316/1‑149 (MQ=255)                                                                                          gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGAAGTAGCGAAAGATTc                                                      >  1:220818/1‑149 (MQ=255)                                                                                           cTCATCAGGCAGGTATAGGGGTCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGCGTCCGGCTATATCGGGAACTTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                                     <  2:167459/149‑1 (MQ=18)                                                                                           cTCAGCAGGCAGGTATAGGGCCCCTTGCGGCCCTGGCGGCGTTGGGCGTATACCGCTAGGAATCAAAAATGATCCGGAAACTACCTGGCTGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                                     <  1:413098/149‑1 (MQ=18)                                                                                            tcatcaGGCAGGTATATGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATTACAGACGCTGCCGTGCCCCCGGTGTTTGCCGGTGCCGTACCCAAGGATTCgg                                                    >  2:371139/1‑149 (MQ=255)                                                                                                   gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAg                                             <  2:404442/149‑1 (MQ=255)                                                                                                    caggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGCTTCGGGACTCAgg                                            >  1:74400/1‑149 (MQ=255)                                                                                                     aggTATAGGGCCCATTGCGGTCCTGGCGGCGTTGGGCGTATCCCGCTAGGACTCCGGCGCTATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGa                                           <  2:306162/149‑1 (MQ=255)                                                                                                       gTATAGTGCCCATTGCGTTCCTGGCTGCGGTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCCACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATg                                         <  1:122145/149‑1 (MQ=18)                                                                                                            gggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgt                                                                                                              <  1:110003/75‑1 (MQ=255)                                                                                                            gggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgt                                                                                                              >  2:110003/1‑75 (MQ=255)                                                                                                               cccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGt                                 <  2:786312/149‑1 (MQ=255)                                                                                                                    tGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATc                            <  1:310855/149‑1 (MQ=255)                                                                                                                           ttGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGGACGTCCCTGCCGGTGTAATGGCAGACGCTGCCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTAAGGATGCTGGCGGTATATCGTCATTg                     >  2:646211/2‑149 (MQ=17)                                                                                                                                 gcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGc               <  1:331641/149‑1 (MQ=255)                                                                                                                                 gcgTGCGGTGTATACCGCTAGGACTGCGGCGATATCCGGAACGGCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGc               <  1:381189/149‑1 (MQ=17)                                                                                                                                    ttCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGtt            >  1:610585/1‑149 (MQ=255)                                                                                                                                    ttCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGtt            <  1:436996/149‑1 (MQ=255)                                                                                                                                        gCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCa                                              >  2:198505/1‑111 (MQ=255)                                                                                                                                        gCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCa                                              <  1:198505/111‑1 (MQ=255)                                                                                                                                           tataACGCTAGGACTCCGGGGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGAt     <  2:673919/149‑1 (MQ=255)                                                                                                                                             taCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATaa   >  1:288324/1‑149 (MQ=255)                                                                                                                                              aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAt  <  2:667837/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAT  >  NC_000913_3_hsa_tpiA/263952‑264241

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: ‑