Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,717	T→G	A163A (GCT→GCG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,717	0	T	G	88.6%	90.4 / ‑0.6	35	A163A (GCT→GCG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base G (19/12);  total (21/14)
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

CACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGA  >  NC_000913/4618481‑4618950                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            cattCCACTGGTTGCTTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTa                                                                                                                                                                                                                               <  1:9/248‑1 (MQ=255) cattCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTa                                                                                                                                                                                                                               <  1:354455/248‑1 (MQ=255)                gTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTAAGCGACCCGTGatat                                                                                                                                                                                                                <  1:273074/251‑1 (MQ=255)                  cGGCGATGGGACTTCGCTGATGTTTGGTCCGATGGTCGCATCCTGCGGCGGCTTTATTCCGATGATCTGTGGTCGCGGTCTCGGTCATACTGGCTGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGGCAGTTTCCCGGATGACAACCGTTT‑CCTCGAAATTAGTAAAGACGTCGGGGTGGCGATTATCGGTCAGACCATCTCACTGG‑CTCCGGCTGATAAACGTTTCTACGCGCCCCCTGatatc                                                                                                                                                                                                               <  2:323521/250‑2 (MQ=255)                           gACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAAcc                                                                                                                                                                                                      <  1:142713/250‑1 (MQ=255)                           gACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCg                                                                                                                                                                                                     <  1:252259/251‑1 (MQ=255)                                       gTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATc                                                                                                                                                                                          <  1:8671/250‑1 (MQ=255)                                               cGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATc                                                                                                                                                                                 <  1:349293/251‑1 (MQ=255)                                                                     cTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAAAGTTTCTACGCGACCCGTGTTATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGGCTCTATTTTGGCga                                                                                                                                                              >  1:425267/1‑248 (MQ=255)                                                                                    ctctGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGAGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGGAGAAACTTGCGGAAGGTc                                                                                                                                            >  1:3722/1‑251 (MQ=255)                                                                                               gCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGCAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGTTATTCCCGCAACCGTGGACTCCCTCCCGCGGGTCACCGCCTCTATTCTGGCGATGAAACTGGCGGGAGGGCTGGAcgcggtg                                                                                                                                 >  1:97849/1‑248 (MQ=255)                                                                                                 cTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGAAGACAACCGTTT‑CCGCGAAATTATTAAAAACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCCCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTggg                                                                                                                               >  1:367642/1‑250 (MQ=255)                                                                                                    ggTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGGCTGGACGCGCTGGGGAt                                                                                                                            >  1:290559/1‑251 (MQ=255)                                                                                                       cATACTGGCGGTACGCTCGACAAACTAGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCCGTCAGACCAACTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATAACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGaagg                                                                                                                          >  1:213759/1‑247 (MQ=255)                                                                                                          aCTGGCGGTACGCTCGACAAACTGGAACCCATCCCGGGCTTCGACATTTCCCCGGATGACAACCGTTC‑CCGCGAAATTATTAAAGACGTCGGTGTGGCGTTTATCGGCCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGCCCCCACCCCTCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGt                                                                                                                      <  1:85794/251‑1 (MQ=255)                                                                                                                 gTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTg                                                                                                               >  1:63378/1‑251 (MQ=255)                                                                                                                     gCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCAcc                                                                                                                                                                              >  1:66024/1‑184 (MQ=255)                                                                                                                     gCTCGACAAACTGGAATCCAACCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGGCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCAGCCCGCTGATCAcc                                                                                                                                                                              <  2:66024/184‑1 (MQ=255)                                                                                                                                              ggttccGACATTTTCTCGGATGACACCCGTTT‑CCGCAAATTTAGTAAAGCCGGCGGCGTGGCGATTATCGGTCAGACCGCCCCACTGG‑CCCCGGCGGATAAACGTTTCTACGCGCCCCGTGATATTACCGCAACCGTGGCCTCCTTCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTAc                                                                                    <  1:232072/244‑1 (MQ=255)                                                                                                                                                  tcgtcATTTCCCCGGATGACAACCTTTC‑CCGAAAAATTATTAAAGCCGTCGTCGTGGCGATTATCGGTCAGACCAGCTCACTGC‑CTCCGGCGGATAAACTTTTCTACGCGCCCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAAct                                                                               <  1:78211/246‑1 (MQ=255)                                                                                                                                                         tttCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTACGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGGGGCGCGTTTATGCCGACCTCCGAAatctccga                                                                         >  1:68503/1‑246 (MQ=255)                                                                                                                                                            cccGAATGACA‑CCGTTTTCCGCGAAATTATAAAAGGCGCGGGCGGGGCGTTTATCGGCCAGCCCAGCTCACTGT‑CTCCGGCGGAAAACCTTTTCTACGCGCCCCGTGATATTCCCGCAACCGTGGACTCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGccc                                                                   <  1:446112/251‑1 (MQ=255)                                                                                                                                                                            tt‑CCGCGAAATGATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCAGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGGGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGGCGCGGTGGTGATGGGCGTGCAAGTGGGGAGCGGCGCGGTTATGCCGATCTGCGAACTCTC‑TGAAGGCCATGCCGacacggttgt                                                   >  2:468943/1‑242 (MQ=255)                                                                                                                                                                                cgcgAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGCCCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCCCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGACGACCTACGAACTCTC‑TGAAGACCCTGCCGA‑AGCGGTTGTTgg                                               >  1:294957/1‑251 (MQ=255)                                                                                                                                                                                  cgAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTCAGAAGTCTC‑TGAAGCCCTTGCCGA‑GACGCTTGTTGGCg                                             >  1:134121/1‑251 (MQ=255)                                                                                                                                                                                   aaaaTTAGTAAAGACGGCGGCGTGGCGTTTAGGGGACAGACC‑GGCCCCTGGGCCCCGGCGAATAACCGTTTCTACGCGACCCGTGATATTACCGCAACCGGGGACTCCACCCCGCGGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGt                                            <  1:328498/250‑1 (MQ=255)                                                                                                                                                                                       ttattaAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGAACCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTCTTCAGGCGAAGAAACTGGCGGAAGGTCTGTACGCGGTGGTGATGGACGTGGAAGGGGGTAGCCGCGCGTTTATGCCGACCTACGCACTCTC‑TGGAGCCCTTGCCGA‑AGCGCTTGTTGGCGTgggt                                        >  1:62730/1‑249 (MQ=255)                                                                                                                                                                                                        gTGGCGGTTATCGGTCAGACCAGCTCACTGG‑ATCTGGCGGGTAAACGTTTCTATGCGACCCGTGATATTACCGCACCCGTGGACTCCCTCCAGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGCAGGTCTGGACGCGCTGGTGATGGACGAGCAAATGGGTAGCGGCGCGGTTATGCCGACCTATGAAC‑CTCGTGAACCCCTTGCCGAGAGC‑ATTTTTGGCGGGGCTAACGGCTCTggggagcg                       >  2:319973/1‑245 (MQ=255)                                                                                                                                                                                                          ggCGATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGGAGGTCTGGACGCGCTGGTGATGGACGTGAAAGGGGGCAGCGGAGCGTTTATGCCGCCCTACGAACTCTC‑TGAAGCCCTTGCCGG‑AGCGATTGTTGGCGGGGCTAACGGCGCTGGCggcggcc                     >  2:427915/1‑245 (MQ=255)                                                                                                                                                                                                             gATTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGcacc                   >  1:147039/1‑250 (MQ=255)                                                                                                                                                                                                              aTTATCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCTCGTTTATGCCGACCTACGAACTCTC‑TGAAGGCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGGGCTGGCGTGCGcaccac                 >  1:275808/1‑251 (MQ=255)                                                                                                                                                                                                                 aTCGGTCAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAACGTGGGTAGCGGCGCGTTTATGCCGACCTACGGACTCTC‑TGAAGCCCTTGCCGG‑AGCGATTGTTGGCGGGGCGAACGGGGCTGGCGTGCGCACCCCcgc              >  2:187896/1‑251 (MQ=255)                                                                                                                                                                                                                        aGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATCTTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTATTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTc       >  1:406497/1‑251 (MQ=255)                                                                                                                                                                                                                          aCCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGGGCCCCACCGCGCTGCTCAc     >  1:66729/1‑251 (MQ=255)                                                                                                                                                                                                                              gctCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTA‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCACGGCTCACCGa  >  1:406831/3‑250 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            CACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTT‑CCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTC‑TGAAGCCCTTGCCGA‑AGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGA  >  NC_000913/4618481‑4618950

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