Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	69,013	A→G	L269L (CTT→CTC)	SAUSA300_RS00310 ←	carbamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	69,013	0	A	G	100.0%	224.6 / NA	84	L269L (CTT→CTC)	SAUSA300_RS00310	carbamate kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/58);  total (26/58)

ACCGTCTAAAGCGTCCCCTAATCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGAAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCGAGACCTCTTTGGTTTTCTTTTCC  >  NC_007793/68917‑69105                                                                                                 |                                                                                             ccgatCTAAAGCGTCCCCTAATCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGc                                                                                           <  1:1957601/96‑1 (MQ=255) aCCGTCTAAAGCGTCCCCTAATCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGc                                                                                           >  1:2985020/1‑100 (MQ=255)        aaaGCGTCCCCTAATCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTc                                                                                    >  1:867244/1‑100 (MQ=255)        aaaGCGTCCCCTAATCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTc                                                                                    >  1:1253505/1‑100 (MQ=255)             gTCCCCTAATCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCtttg                                                                               <  1:2334784/100‑2 (MQ=255)             gTCCCCTAATCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCtttg                                                                               <  1:1691438/100‑2 (MQ=255)               ccccTAATCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTGTTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGc                                                                             >  1:4255855/1‑100 (MQ=255)               ccccTAATCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGc                                                                             >  1:271176/1‑100 (MQ=255)               ccccTAATCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGc                                                                             >  1:2033924/1‑100 (MQ=255)               ccccTAATCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGc                                                                             >  1:4591169/1‑100 (MQ=255)                 cccaaTCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCaa                                                                           <  1:4225110/97‑1 (MQ=255)                 ccTAGTCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCaa                                                                           >  1:5145183/1‑100 (MQ=255)                 ccTAATCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCaa                                                                           >  1:872414/1‑100 (MQ=255)                    aaTCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAAtt                                                                        >  1:1599646/1‑100 (MQ=255)                    aaTCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAAtt                                                                        >  1:3678693/1‑100 (MQ=255)                    aaTCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAAtt                                                                        <  1:372250/100‑1 (MQ=255)                     aTCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTg                                                                       >  1:594977/1‑100 (MQ=255)                      tCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGa                                                                      <  1:4978963/100‑1 (MQ=255)                      gccTGCTAGAGATGTCATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGa                                                                      >  1:5488385/2‑100 (MQ=255)                       ccTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGAc                                                                     >  1:2019393/1‑100 (MQ=255)                       ccTGCTAGAGATGTAATCAAAACACTGCCTTTAGTA‑TTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATCCTTCCTTTAGCAAATTGAcc                                                                    <  1:4347733/100‑1 (MQ=255)                         tGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCa                                                                   >  1:3859704/1‑100 (MQ=255)                          gCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCac                                                                  >  1:1438499/1‑99 (MQ=255)                          gCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCac                                                                  >  1:852660/1‑99 (MQ=255)                          gCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCAt                                                                  >  1:2949629/1‑100 (MQ=255)                          gCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCAt                                                                  >  1:3773538/1‑100 (MQ=255)                          gCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCAt                                                                  >  1:1718313/1‑100 (MQ=255)                          gCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCAt                                                                  >  1:1456079/1‑100 (MQ=255)                            tAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCagag                                                                <  1:1076413/100‑4 (MQ=255)                            tAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCa                                                                <  1:3676982/100‑1 (MQ=255)                            tAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCa                                                                <  1:2013231/100‑1 (MQ=255)                             agcgaTGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCag                                                               <  1:1179676/97‑1 (MQ=255)                             agagATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCGGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCag                                                               <  1:3165789/100‑1 (MQ=255)                             agagATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGGCCATCag                                                               <  1:235449/100‑1 (MQ=255)                             agagATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCag                                                               >  1:2893691/1‑100 (MQ=255)                             agagATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCag                                                               <  1:910424/100‑1 (MQ=255)                             agagATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCag                                                               <  1:1250542/100‑1 (MQ=255)                             agagATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCag                                                               <  1:691787/100‑1 (MQ=255)                             agagATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCag                                                               <  1:1655963/100‑1 (MQ=255)                             agagATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCag                                                               <  1:5316788/100‑1 (MQ=255)                             agagATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCag                                                               <  1:4923068/100‑1 (MQ=255)                             agagATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCag                                                               <  1:2136283/100‑1 (MQ=255)                             agagATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCag                                                               <  1:3677940/100‑1 (MQ=255)                             agagATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCag                                                               <  1:343802/100‑1 (MQ=255)                             agagATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCACCag                                                               >  1:3699372/1‑100 (MQ=255)                             agagATGTAATCAAAACACTGCCTTTAGTATTTGTTTCAAGAAATTGAAGTGCAGCGTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCag                                                               <  1:4884825/100‑1 (MQ=255)                                     aaTCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGt                                                       <  1:3439290/100‑1 (MQ=255)                                             tacTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTtttcat                                               >  1:3798212/2‑100 (MQ=255)                                                tGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTtttcatttc                                            <  1:4849125/100‑1 (MQ=255)                                                    tttAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATcc                                        <  1:33680/100‑1 (MQ=255)                                                     ttAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCca                                       >  1:4679949/1‑100 (MQ=255)                                                      tAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCGAATTGACCATCAGAGATATGTTTTTTCATTTCATCcac                                      <  1:730312/100‑1 (MQ=255)                                                      tAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCcac                                      <  1:1754963/100‑1 (MQ=255)                                                      tAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCcac                                      <  1:4753398/100‑1 (MQ=255)                                                       aGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCcaca                                     <  1:3712928/100‑1 (MQ=255)                                                        gTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCcacaa                                    <  1:2858117/100‑2 (MQ=255)                                                         tATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAgc                                   <  1:2218780/100‑2 (MQ=255)                                                          aTTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGac                                  <  1:5210805/100‑1 (MQ=255)                                                          aTTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGAAATGTTTTTTCATTTCATCCACAGac                                  <  1:899805/100‑1 (MQ=255)                                                                ttCAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCa                            <  1:3994027/100‑1 (MQ=255)                                                                  cAAGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCAtg                          <  1:308648/100‑2 (MQ=255)                                                                   aaGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCg                         <  1:3383859/100‑1 (MQ=255)                                                                   aaGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCg                         <  1:3348273/100‑1 (MQ=255)                                                                   aaGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCg                         <  1:1295741/100‑1 (MQ=255)                                                                   aaGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCg                         <  1:5511350/100‑1 (MQ=255)                                                                   aaGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCg                         <  1:5506924/100‑1 (MQ=255)                                                                   aaGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCg                         <  1:2423434/100‑1 (MQ=255)                                                                   aaGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCg                         <  1:3999698/100‑1 (MQ=255)                                                                   aaGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCg                         <  1:3481196/100‑1 (MQ=255)                                                                   aaGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCg                         <  1:1912767/100‑1 (MQ=255)                                                                   aaGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCg                         <  1:3833136/100‑1 (MQ=255)                                                                   aaGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCg                         <  1:3871137/100‑1 (MQ=255)                                                                   aaGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCg                         <  1:4608967/100‑1 (MQ=255)                                                                   aaGAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCg                         <  1:1999727/100‑1 (MQ=255)                                                                     gAAATTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCgag                       <  1:1063876/100‑1 (MQ=255)                                                                      aaaTTGAAGTGCAGCTTCAACTTTTGGGAGCATTCTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCgaga                      <  1:4568218/100‑1 (MQ=255)                                                                        aTTGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCGAGAcc                    <  1:3443679/100‑1 (MQ=255)                                                                         ttGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCGAGACct                   <  1:3773729/100‑1 (MQ=255)                                                                         ttGAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCGAGACct                   <  1:3886868/100‑1 (MQ=255)                                                                           gAAGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCGAGACctag                 <  1:5420411/100‑3 (MQ=255)                                                                            aaGTGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCGAGACctcct                <  1:1176451/100‑3 (MQ=255)                                                                               tGCAGCTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCGAGACCTCTTTgg             <  1:291778/100‑1 (MQ=255)                                                                                    cTTCAACTTTTGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCGAGACCTCTTTGGttttc        <  1:2265037/100‑1 (MQ=255)                                                                                          agcctGGGAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCGAGACCTCTTTGGTTTTCTTTTcc  >  1:2864134/5‑100 (MQ=255)                                                                                                 |                                                                                             ACCGTCTAAAGCGTCCCCTAATCCTGCTAGAGATGTAATCAAAACACTGCCTTTAGTATTTTTTTCAAGAAATTGAAGTGCAGCTTCAACTTTTGGAAGCATACTTCCTTTAGCAAATTGACCATCAGAGATATGTTTTTTCATTTCATCCACAGACACTTCATCGAGACCTCTTTGGTTTTCTTTTCC  >  NC_007793/68917‑69105

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