Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,123,556	C→T	E167K (GAA→AAA)	grxB ←	reduced glutaredoxin 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,123,556	0	C	T	85.7%	79.4 / 7.6	35	E167K (GAA→AAA)	grxB	reduced glutaredoxin 2
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base T (14/16);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

CACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGC  >  NC_000913/1123328‑1123802                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       cACCAGCGAGAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGttt                                                                                                                                                                                                                                   <  1:272210/250‑1 (MQ=255)  aCCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAACCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGa                                                                                                                                                                                                                                 <  2:168639/251‑1 (MQ=255)                  ccAGCAGAAAAAACTGAAAACGGCGGCAACCCGCTAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAg                                                                                                                                                                                                                 <  2:237395/251‑1 (MQ=255)                           aaaaCTGAAAACAGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGAGTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTGGACGATCAGTTTGTCCAGCGCACGTaa                                                                                                                                                                                                        <  2:429529/251‑1 (MQ=255)                                             aaCCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTTGCCCGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGATGAATATCATCTTCCGAAAGTTTGCCCTTCACGGCGTTCGGGTTTACGGTCAATTTGTCCAACGCCCGTAAAACATCGCTGGTATTCtt                                                                                                                                                                                      >  1:193853/1‑251 (MQ=255)                                                ccGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAAACGCCATTGATGATAACAAAATGCTTTGTGTCTGTTTCGCCATATTATCGCGGGAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCCGTAGCGGGAAATGCTGAATATCATCTTCCCAAAGGTTGCCATTCCCCGCGTTTCGTTTGACGATCAGTTTTTTCAGCACACGTAAATCATCGCTGATATTCTTAAt                                                                                                                                                                                   >  1:343295/1‑251 (MQ=255)                                                  gcgAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCa                                                                                                                                                                                 <  1:156804/251‑1 (MQ=255)                                                  gcgAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCa                                                                                                                                                                                 <  1:178051/251‑1 (MQ=255)                                                  gcgAAGACAGGCTGCCGTGGCGGGTCAGGATTAAAGCGCCATTGATGATAACAAATTGAATGGTGTCTGTTTCGCCATAGTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGGCAGATTACGCAGTAGCGGGAATAGCTGAAGATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATGCTTAATCa                                                                                                                                                                                 <  2:467483/251‑1 (MQ=255)                                                        aCAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCg                                                                                                                                                                            <  1:393643/250‑1 (MQ=255)                                                             gtGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCagag                                                                                                                                                                      <  2:181100/250‑1 (MQ=255)                                                                    ggcggGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATTGCGGTAATAAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGGAAAACTGCATATCATCCTCCGAAAGATTTCCATTCACGGGGGTTGGTTTGAAGATCAGTTTTTCCCGCGCCCGTAAAACATCGCCGGTATTccaatcggaaacgcggcgtggggggg                                                                                                                                                               >  2:97078/1‑226 (MQ=255)                                                                    ggcggGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTc                                                                                                                                                                                        <  1:97078/226‑1 (MQ=255)                                                                            aGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACAGTCAGAGTGGGCCAGCAGGTCg                                                                                                                                                       >  1:117616/1‑251 (MQ=255)                                                                                  aaaTAGCCATTGAAGATAACAAATTGATTTGTGTCTGTTTCGCCATCTTATCGAGGTAATCAGCAACGCGGCTTGGGCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGAGGGAAAAGCAGAATAACATTTTCCGAAAGGTTGCCAATCACGGAGGTCGGTTTGACGGTCAGTTTGTCATACGCACGGAAATCATCGATTATATTCTTTATCAGTACGTCagagaggtgaaacaggtaggcaaa                                                                                                                                                  >  2:374567/1‑230 (MQ=255)                                                                                          aTTGATGATAACAAATTGATTTTTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAAGCATCGCTGATATTCTTAATCAGACCGGCAGAGTGGGCCAGCAGGTCGTCAAAATTACCcgc                                                                                                                                         <  1:320403/251‑1 (MQ=255)                                                                                             gatgatAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGGTTGTCCAGCGCACGTAAATCATCGCTGTTATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCgcgcc                                                                                                                                      >  2:357228/1‑250 (MQ=255)                                                                                                           gATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGt                                                                                                                        >  1:447218/1‑251 (MQ=255)                                                                                                                          cGCCATAATATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTgc                                                                                                         <  2:158122/251‑1 (MQ=255)                                                                                                                              atatTATCGCGGTAATCAGCACCGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCTATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGAGCACGTAAATCATCGCTGATATTCTTAATCAGGCCGTCAGTGTGGGCCCGCAGGTCGGCAAAATTACCAGACCTGGCCTCTTTCTTGTCGACGaaaagtttgagcac                                                                                                     >  2:311081/1‑240 (MQ=255)                                                                                                                                attatCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGcgcgcg                                                                                                    <  1:357228/250‑1 (MQ=255)                                                                                                                                         gTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGGCAGATTTCGCAGTAGCGGGAATAGCGGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTGCGGTTTGACGATCAGTTGGTCCAGCGGACGTAAATCATGGCTGATATGCTTAATCAGACCGTCAAATGGGGCCAGCAGGTCGGCAAAATGACCCGCGCTGGCCTCTGTCTGGTCGACGAAATATTTGCGCGCGGCGGGAGt                                                                                            <  2:59400/249‑1 (MQ=255)                                                                                                                                              cAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGCATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAACTCATCGCTGATATTCTTAATCCGACCGTCAGAGTGGGCCAGCCGGTCGGCAAAATTACCCGCGCTGGCATCTTTCTTGTTGACGAAATTTTTGCGCGCGGGGGGAggggaaaa                                                                                      >  2:376576/1‑243 (MQ=255)                                                                                                                                               aGCAACGCGGCTTGGCCAATTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATCTTTGCGCGCGGCGGGAGTAGAAAACt                                                                                    >  2:382852/1‑251 (MQ=255)                                                                                                                                                    cgcgGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACtcatca                                                                               <  1:298765/251‑1 (MQ=255)                                                                                                                                                           tGGCCAGTGAATTCCGGCCACCAGCGTCAGAGTACGCATGAGCGGGAATAGCTGAATATCATCGTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCGGATATTCTTAAGCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCAACAAATGcc                                                                         <  2:411158/250‑1 (MQ=255)                                                                                                                                                            ggCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCg                                                                        <  1:90689/250‑1 (MQ=255)                                                                                                                                                                         cGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCAGCAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGCTATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCCACAAATGCCGATTTGGCAAAAcgc                                                          >  1:45164/1‑251 (MQ=255)                                                                                                                                                                                  gCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAAcagc                                                 <  1:238442/251‑1 (MQ=255)                                                                                                                                                                                       aGATTACGCAGTAGCGGGAAGTGCGGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGAAATTTGTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGttt                                            <  2:124784/251‑1 (MQ=255)                                                                                                                                                                                                 gTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGCAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTggagt                                    >  1:35905/1‑246 (MQ=255)                                                                                                                                                                                                         aaTAGCTGAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACACGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGGGGAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGGTTGTTGGAGTAGCAATTGAc                          >  2:117656/1‑251 (MQ=255)                                                                                                                                                                                                                gAATATCATCTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTgcg                    <  2:447218/250‑1 (MQ=255)                                                                                                                                                                                                                     tcatcTTCCGAAAGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGccc              >  1:466663/1‑250 (MQ=255)                                                                                                                                                                                                                                 aGTTTGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGCCCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGc  >  1:336899/1‑251 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       CACCAGCGATAAACTCGCCAGCAGAAAAAACTGAAAACGGCGGCAACCCGCGAATACAGGCTGCCGCGGCGGGTCAGGATTAAATCGCCATTGATGATAACAAATTGATTTGTGTCTGTTTCGCCATATTATCGCGGTAATCAGCAACGCGGCTTGGCCAGTTAATTCCGGCTACCAGCGTCAGATTACGCAGTAGCGGGAATAGCTGAATATCATCTTCCGAAAGTTCGCCATTCACGGCGTTCGGTTTGACGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGC  >  NC_000913/1123328‑1123802

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