Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,066,699	C→T	A372V (GCC→GTC)	agp →	glucose‑1‑phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,066,699	0	C	T	97.2%	100.3 / NA	36	A372V (GCC→GTC)	agp	glucose‑1‑phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/17);  minor base G (0/1);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGACAAAGCTCTGGTCACCGATCGCACCAGCGCACCGAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAGAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGCCCCAAAAAAGCGAATCTTTATCGGTCTTAGCAAAGGCTTTG  >  NC_000913/1066455‑1066936                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              atCGACAAAGCTCTGGTCACCGATCGCACCAGCGCACCGAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATg                                                                                                                                                                                                                                          <  1:48562/250‑1 (MQ=255)  tCGACAAAGCTCTGGTCACCGATCGCACCAGCGCACCGAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCg                                                                                                                                                                                                                                        >  2:100336/1‑251 (MQ=255)      cAAAGCTCTGGTCACCGATCGCACCAGCGCACCGAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGACCAGTTACGTAATGTCGATGCGTTa                                                                                                                                                                                                                                     >  1:307949/1‑250 (MQ=255)          gCTCTGGTCACCGATCGCACCAGCGCACCGAAAATTACAGGGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCt                                                                                                                                                                                                                                <  1:432431/251‑1 (MQ=255)                           cACCAGCGCACCGAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGCTTTGATGAAAATTGAATATGTGTATCAGAGTGCGGCGCAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGc                                                                                                                                                                                                               >  2:230701/1‑251 (MQ=255)                           cACCAGCGCACCGAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTACCCCTGCAGGCCCCTGCGCAGc                                                                                                                                                                                                               >  2:196713/1‑251 (MQ=255)                            aCCAGCGCACCGAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGGTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCg                                                                                                                                                                                                              <  1:137935/251‑1 (MQ=255)                             ccAGCGCACCGAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCgg                                                                                                                                                                                                             >  1:335948/1‑250 (MQ=255)                                  gcACCGAAAATTACAGTGTTGGGTGGGCACGACTCCACCATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCACTCGCGTTTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGAGTGAc                                                                                                                                                                                                        >  2:22416/1‑251 (MQ=255)                                       gAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTgg                                                                                                                                                                                                   >  1:272992/1‑251 (MQ=255)                                                      ggttggGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGCAGCTGGAATTAAGCGGTTGcc                                                                                                                                                                                    >  1:412031/1‑251 (MQ=255)                                                      ggttggGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGc                                                                                                                                                                                     <  1:84158/250‑1 (MQ=255)                                                               cGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACg                                                                                                                                                                           <  1:274085/251‑1 (MQ=255)                                                                     cAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTAGCAGTTGCATGACCAGAACGATCGCACGCCGATTGGCGGCAAAAGCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATg                                                                                                                                                                     <  2:398639/251‑1 (MQ=255)                                                                         aTTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGttt                                                                                                                                                                 <  2:126187/251‑1 (MQ=255)                                                                         aTTGCCCCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTATCCAGCGTGGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATGGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCTTTAACCCGGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATgggtt                                                                                                                                                                 <  2:53555/251‑4 (MQ=255)                                                                            gCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTg                                                                                                                                                                                      <  1:315022/227‑1 (MQ=255)                                                                            gCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTg                                                                                                                                                                                      >  2:315022/1‑227 (MQ=255)                                                                                   tGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATg                                                                                                                                                       >  1:23731/1‑251 (MQ=255)                                                                                       aaCGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATa                                                                                                                                                   <  2:264760/251‑1 (MQ=255)                                                                                        aCGGCGCTGGAGTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGGGTGACGCTGGAATTCAGCGGTTGCCCGATAGACGCTGATGGGTTCTGCCCGATggcaa                                                                                                                                                   >  2:75955/1‑247 (MQ=255)                                                                                                              accAGTTGCATGTCCAGAACGAACGCACGCCGATTGGCGGCATAATCGTGTGCCAGCGTGTGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAAGATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGAGGCGTTAACCCTGCAGGCACCTGATCATCGTGTGACGTTGGAATTAAGCGGTTGCCCGATAGACGGTGATGGTTTCGGCCCGATGGATTAGTTTTATAGCGTGTtgaatgat                                                                                                                            <  2:210338/249‑2 (MQ=255)                                                                                                                cAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAg                                                                                                                           <  1:471132/250‑1 (MQ=255)                                                                                                                         gACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGCTAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAAt                                                                                                                  >  1:20186/1‑250 (MQ=255)                                                                                                                         gACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAAt                                                                                                                  <  1:77139/250‑1 (MQ=255)                                                                                                                                              aTTGGCGGCAAAATCGGTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTAAGTAATGTCGATGCGTTAAACCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGAGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGGATGAAGCGGTGAAATAACAAAAAACTcccccccggg                                                                                             >  2:271377/1‑245 (MQ=255)                                                                                                                                                gggcggcAAAATCGTTTTCCAGCGTGGACATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGGGTATCAGAGTGCGGAACAGTTACGTAATGGCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAGAAAACTCCCCCGCGAGAAg                                                                                          <  2:335948/250‑1 (MQ=255)                                                                                                                                                         aaTCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAGAAAACTCCCCCGCGAGAAGCGGGGGAGt                                                                                 >  2:342872/1‑251 (MQ=255)                                                                                                                                                                      cGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAGAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACg                                                                    <  1:228651/251‑1 (MQ=255)                                                                                                                                                                                                 gATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAGAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGccc                                         <  1:151521/251‑1 (MQ=255)                                                                                                                                                                                                     tgatgaAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAGAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGCCCCaa                                      >  1:482747/1‑250 (MQ=255)                                                                                                                                                                                                                          aTCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAGAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGCCCCAAAAAAGCGAATCTTTATCGGTCt                >  2:55990/1‑251 (MQ=255)                                                                                                                                                                                                                                gTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAGAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGCCCCAAAAAAGCGAATCTTTATCGGTCTTAGCaa          <  1:100336/251‑1 (MQ=255)                                                                                                                                                                                                                                  gCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTtgaat                                                                                                                               >  2:167992/1‑132 (MQ=255)                                                                                                                                                                                                                                  gCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTtgaat                                                                                                                               <  1:167992/132‑1 (MQ=255)                                                                                                                                                                                                                                        cAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAGAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGCCCCAAAAAAGCGAATCTTTATCGGTCTTAGCAAAGGCTTTg  >  1:37922/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              ATCGACAAAGCTCTGGTCACCGATCGCACCAGCGCACCGAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAGAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGCCCCAAAAAAGCGAATCTTTATCGGTCTTAGCAAAGGCTTTG  >  NC_000913/1066455‑1066936

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