Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	23,520	Δ1 bp	6.2%	coding (250/624 nt)	dgk ←	deoxyguanosine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	23,520	0	T	.	6.2%	207.1 / 5.7	64	coding (250/624 nt)	dgk	deoxyguanosine kinase
Reads supporting (aligned to +/- strand):  ref base T (21/38);  new base . (2/2);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

CTTTATGTAAATGATGAAATTCGGCTTTGGCAGATCATCCGTTAGCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTCCTTTATATTGTCATAAAACTTATCAAGGTACGGGTTATC  >  BS168_reference_genome/23383‑23655                                                                                                                                          |                                                                                                                                        cTTTATGTAAATGATGAAATTCGGCTTTGGCAGATCATCCGTTAGCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTAtt                                                                                                                                        <  8:231186/139‑1 (MQ=255)    tATGTAAATGATGAAATTCGGCTTTGGCAGATCATCCGTTAGCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTa                                                                                                                                       >  6:198705/1‑137 (MQ=255)    tATGTAAATGATGAAATTCGGCTTTGGCAGATCATCCGTTAGCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTa                                                                                                                                       <  4:186027/137‑1 (MQ=255)           atgatgAAATTCGGCTTTGGCAGATCATCCGTTAGCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTA‑TACAGGCTGTcc                                                                                                                             >  3:12153/1‑139 (MQ=255)           atgatgAAATTCGGCTTTGGCAGATCATCCGTTAGCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTA‑TACAGGCTGTcc                                                                                                                             >  3:8168/1‑139 (MQ=255)                              gCAGATCCTCCCTTAGCATGTTATATATCTTTTTGTATTTTTCCAGATGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGAt                                                                                                            <  3:45791/138‑1 (MQ=255)                              gCAGATCATCCGTTAGCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGAt                                                                                                            <  3:45795/138‑1 (MQ=255)                                  atcatcCGTTAGCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTg                                                                                                       <  1:288477/139‑1 (MQ=255)                                   tcatcCGTTAGCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTgg                                                                                                      <  2:256416/139‑1 (MQ=255)                                            aGCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCtttttt                                                                                                                      >  6:158432/1‑114 (MQ=255)                                            aGCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCtttttt                                                                                                                      <  5:158432/114‑1 (MQ=255)                                            aGCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCa                                                                                              <  5:113968/138‑1 (MQ=255)                                            aGCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGCTCGGTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCa                                                                                              <  2:309752/138‑1 (MQ=255)                                             gCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCGCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAg                                                                                            <  5:198705/139‑1 (MQ=255)                                             gCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAg                                                                                            <  1:15354/139‑1 (MQ=255)                                             gCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAg                                                                                            <  3:20211/139‑1 (MQ=255)                                             gCAGGTGATATATCTTTTTGTATATTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAg                                                                                            <  2:268889/139‑1 (MQ=255)                                               aGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCt                                                                                          <  2:75314/139‑1 (MQ=255)                                                 gTGATATATCTTTTTGTATTTTTCCAGCTGATTTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTAATACAGTCTGTCCCTTTTTAAAAATGAGATCAAAGGTTTATTCAAGCtct                                                                                        >  5:34628/1‑137 (MQ=255)                                                 gTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTTTCTTCAAGCTGt                                                                                        >  3:256470/1‑139 (MQ=255)                                                 gTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGt                                                                                        >  3:125262/1‑139 (MQ=255)                                                 gTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGt                                                                                        >  8:182835/1‑139 (MQ=255)                                                 gTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGt                                                                                        >  8:133360/1‑139 (MQ=255)                                                  tGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCt                                                                                          >  6:235633/1‑136 (MQ=255)                                                  tGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCt                                                                                          <  5:235633/136‑1 (MQ=255)                                                     tatatCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGt                                                                                    <  6:211202/139‑1 (MQ=255)                                                            ttttGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTAt                                                                                  <  8:306287/134‑1 (MQ=255)                                                            ttttGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTAt                                                                                  >  7:306287/1‑134 (MQ=255)                                                             tttGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGa                                                                            >  5:108277/1‑139 (MQ=255)                                                             tttGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTGTGTATCTGTGa                                                                            >  5:174229/1‑139 (MQ=255)                                                                     tttCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGaa                                                                     >  4:242009/1‑138 (MQ=255)                                                                                     agaCAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTgg                                                    <  3:242009/139‑1 (MQ=255)                                                                                      gaCAAGGGTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCGGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGa                                                   <  1:193207/139‑1 (MQ=255)                                                                                        cAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAg                                                 <  7:182835/139‑1 (MQ=255)                                                                                           ggTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTc                                              >  5:89692/1‑139 (MQ=255)                                                                                            gttcgttcCGCAAAAATCACATTTTTATAAATATGATAATCAGCTA‑TACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCa                                            <  4:8168/139‑1 (MQ=255)                                                                                            gttcgttcCGCAAAAATCACATTTTTATAAATATGATAATCAGCTA‑TACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGGTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCa                                            <  4:12153/139‑1 (MQ=255)                                                                                             ttcgttcCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCa                                            <  4:125262/139‑1 (MQ=255)                                                                                             ttcgttcCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCa                                            <  6:89692/139‑1 (MQ=255)                                                                                             ttcggtcCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCa                                            <  1:59641/139‑1 (MQ=255)                                                                                             ctcgttcCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAACAACATCTCGAGCTGGAAGCTCCa                                            <  8:59778/138‑1 (MQ=255)                                                                                             atcgctcCGCAAAAATCACATTTTTATAAATAAGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCa                                            <  5:72634/138‑1 (MQ=255)                                                                                              tcgttcCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCa                                            <  6:29140/138‑1 (MQ=255)                                                                                              tcgttcCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCa                                            >  5:29140/1‑138 (MQ=255)                                                                                              tcgttcAGCAAAAATGACATTTTTATAAATATGATAATAAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTTGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGTAAAAACATCTCGAGCTGGAAGCTCCa                                            <  6:34628/138‑1 (MQ=255)                                                                                                 ttcCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCt                                                                                                 <  8:183813/82‑1 (MQ=255)                                                                                                 ttcCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCt                                                                                                 >  7:183813/1‑82 (MQ=255)                                                                                                 ttcCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTcc                                        >  8:116934/1‑139 (MQ=255)                                                                                                        aaaaTCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTCCTTTatat                                 <  6:176259/139‑1 (MQ=255)                                                                                                        aaaaTCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTCCTTTatat                                 <  5:219144/139‑1 (MQ=255)                                                                                                              acaTTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTCCTTTATATTGTCAt                           <  7:97098/139‑1 (MQ=255)                                                                                                                          atatGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAg                                                                                            >  8:173482/1‑62 (MQ=255)                                                                                                                          atatGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAg                                                                                            <  7:173482/62‑1 (MQ=255)                                                                                                                             tGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTCCTTTATATTGTCATAAAACTTATCAAGGt            <  7:116934/139‑1 (MQ=255)                                                                                                                                  atcatcTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTCCTTTATATTGTCATAAAACTTATCAAGGTACggg       <  1:41500/139‑1 (MQ=255)                                                                                                                                     aGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTCCTTTATATTGTCATAAAACTTATCa                >  2:183928/1‑127 (MQ=255)                                                                                                                                     aGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTCCTTTATATTGTCATAAAACTTATCa                <  1:183928/127‑1 (MQ=255)                                                                                                                                     aGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTCCTTTATATTGTCATAAAACTTATCAAGGTACGGGtt     <  6:174229/138‑1 (MQ=255)                                                                                                                                     aGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTCCTTTATATTGTCATAAAACTTATCAAGGTACGGGTTa    >  4:232125/1‑139 (MQ=255)                                                                                                                                      gCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAg                                                        <  7:268424/86‑1 (MQ=255)                                                                                                                                      gCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAg                                                        >  8:268424/1‑86 (MQ=255)                                                                                                                                      gCTATTACACGCCGTCACTTTTTAATAAAGTGAGGAGTGGTGTCGTCAAGCTGTGTGTAGCTGTGACAAAGAAACTACATCTCGAGCTGGTAGCTCCAATCCATTATTTTGTCATAAAACTTATCAAGGTACGGGTTAt   <  6:274017/139‑1 (MQ=255)                                                                                                                                       cTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTCCTTTATATTGTCATAAAACTTATc                 <  6:87709/124‑1 (MQ=255)                                                                                                                                       cTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTCCTTTATATTGTCATAAAACTTATc                 >  5:87709/1‑124 (MQ=255)                                                                                                                                       cTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTCCTTTATATTGTCATAAAACTTATCAAGGTACGGGTTATc  >  5:120308/1‑139 (MQ=255)                                                                                                                                          |                                                                                                                                        CTTTATGTAAATGATGAAATTCGGCTTTGGCAGATCATCCGTTAGCAGGTGATATATCTTTTTGTATTTTTCCAGCTGATGTGGAGACAAGGTTCGTTCCGCAAAAATCACATTTTTATAAATATGATAATCAGCTATTACAGGCTGTCCCTTTTTTAAAAAGTGATCACTGGTGTCTTCAAGCTGTTTGTATCTGTGACAAAGAAAAAACATCTCGAGCTGGAAGCTCCATTCCTTTATATTGTCATAAAACTTATCAAGGTACGGGTTATC  >  BS168_reference_genome/23383‑23655

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