Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	557,635	A→T	L36Q (CTG→CAG)	folD ←	bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	557,635	0	A	T	97.2%	111.8 / NA	36	L36Q (CTG→CAG)	folD	bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (13/22);  minor base G (0/1);  total (13/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTT‑TGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCG  >  NC_000913/557407‑557855                                                                                                                                                                                                                                     |                                                                                                                                                                                                                              cccaccttttggacgtaacaaacccccGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCgg                                                                                                                                                                                                          <  1:179642/226‑1 (MQ=255)  cagatatttgcatatactCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGTGCGGGAGACGAACCCGACTTCTTCCCAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCt                                                                                                                                                                                                        <  1:254114/234‑1 (MQ=255)        gttagtctaacacattccccGCCGGTAACGGCAGTTTAACCAGAATGACATCGTGGGTGTTGTCGGCATTCAGGGTATCGATAAGCTCCATCAGCTCCGCTTCGCGGGTGGTTGCCGGTAGGTCATAAGAGCGGGAGACTAACCCGACTGCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGATTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCaaga                                                                                                                                                                                                  <  2:84046/234‑3 (MQ=255)          ttgatgtcaTCATTACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATAc                                                                                                                                                                                                 <  1:99543/242‑1 (MQ=255)          tctggttagaaaaaccccGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATGCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATGTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTTCAATAc                                                                                                                                                                                                 <  2:414958/235‑1 (MQ=255)             aCGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGGATCGATAAGCTCCAGCAGCTCCGCGTCGCTGGTGGTTTCCGGGCGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGAAGGGTTACTACCCACCAGCCCAAAGGCCTGTCCTGGTGCCCGCAGTCCTGCTGCCATACgcgc                                                                                                                                                                                             >  2:340560/1‑251 (MQ=255)                  aTCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCCCAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGa                                                                                                                                                                                         >  1:422262/1‑250 (MQ=255)                   tCAATACCCGCCGGTAACTGCAGTCGAAACAGAATTCCATCGATGGTGCTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGCCATAAGAGCGGGAGACGAACCCGCCTTCTTCCCAAGCCCTGCGGTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCCGCACCACAGCCTGTCCTTGTGACCGCAGTCCGGCTGCCATCCCCGCCTGAAc                                                                                                                                                                                       >  2:165643/1‑251 (MQ=255)                        ccccGCCGGTAACGGCAGTTGAACCAGAATGGCATCGATGCTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCCGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGTAGACGAGCCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACttttt                                                                                                                                                                                  <  1:227173/250‑1 (MQ=255)                        aCCCGCCGGTAACGGCAGTTGACCCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACttttt                                                                                                                                                                                  >  1:406944/1‑251 (MQ=255)                         cccGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTg                                                                                                                                                                                 <  1:235698/251‑1 (MQ=255)                                   aCGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCa                                                                                                                                                                       <  1:364447/251‑1 (MQ=255)                                                  gAATGCCATCGATGGTGGTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGTTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCtgctgc                                                                                                                                                        <  1:226231/251‑1 (MQ=255)                                                           cGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGt                                                                                                                                                <  2:109822/250‑1 (MQ=255)                                                             aTGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGtttt                                                                                                                                             <  1:177542/251‑1 (MQ=255)                                                               ggTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTAc                                                                                                                                           >  1:16989/1‑251 (MQ=255)                                                                 ggttgtCGGCATTAGGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTTCGCAATCGTTTTACCg                                                                                                                                         <  2:406944/250‑1 (MQ=255)                                                                  gttttCGGCATTCAGCTGATCGGTAAGCTCCAGCAGCTCCGCTTCGCTGTTGGTTTCCGTGAGGTCATAAGAGCGGGAGACGACCCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCGCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTTAGCAACTTTAGAGCGCACCTGCTGCGCAATCTGTTTACCGt                                                                                                                                        <  2:156489/247‑1 (MQ=255)                                                                            ttCAGCGAATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGAGGTTTCCGGGAGGTAATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAGCTTTTTGAGCAACTTCAGAGCGCACCTGCGGCGCCATCGTTTTACCGTCAATAATCtt                                                                                                                              >  2:31602/1‑251 (MQ=255)                                                                              cAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATCATCTT‑tg                                                                                                                           >  2:94081/1‑251 (MQ=255)                                                                                      cGATAAGCTCCAGCAGCTCCGATTCTCTGTTGGTTTCCGGTAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTGGCGGTTGCTTTCGACATAAATTTGCGATGCATGGTTACTACCCACCTGCACAACGTCCTTTCCTGGTGCCCGCATTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACCTCAGAGCGCACCTGGTGCGCAATCGTTTTACCGTCAATAATCTT‑TGCTGCCAt                                                                                                                    <  1:449327/250‑1 (MQ=255)                                                                                              cccAGCCGCTGCGCTGCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTT‑TGCTGCCATCagagagag                                                                                                           <  2:305004/250‑1 (MQ=255)                                                                                                                            gggAGGTCATAAGAGCGGGAGACGAACCCGGCTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCGGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGcc                                                                                                                                                                                                                    <  1:323216/117‑1 (MQ=255)                                                                                                                            gggAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGATTGCGACATACCTTTGCGATGCAGGGTTACTACCCACCAGCACCACGGCCTGTCCTGGTGcc                                                                                                                                                                                                                    >  2:323216/1‑117 (MQ=255)                                                                                                                              gAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATTAATTGGCAATGCATGGTTTCTACCCACCAGCACTACGGCCGGTCCTGGTGCCCGCAGTCCGGCTTCAATACGCGCCGTATCTTTTTGAGCAGCGTCATAGCTCACCTGCTGCGCAATCGTGTTACCGTTAAGAATCGG‑TGTTGCCATCAGAGAGAGGATTCTATCTGTTACTGAGATCGAATGGGATg                                                                           <  2:45807/251‑1 (MQ=255)                                                                                                                                      aaGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTT‑TGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAtt                                                                   >  1:119901/1‑251 (MQ=255)                                                                                                                                                gaCGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTT‑TGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGa                                                           >  1:260703/1‑249 (MQ=255)                                                                                                                                                                tcACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTT‑TGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCa                                          <  1:98391/250‑1 (MQ=255)                                                                                                                                                                     aGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCAACAAGCACAACGGCCTGTCCTGGTGCCCACAGTCCGGCTGCAATACGCGCCTGCACTTTTTGGCCACCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATACTCTT‑TGCAGCCATCAGAGAGAGGATTCCATCTGTTACGTCGATCGAAGGGGATGAGCCTATTTTGTCAGAAGCGGGGCGCGTTCTCAGgtttc                                    >  1:248055/1‑251 (MQ=255)                                                                                                                                                                          tGCGGTTGCTGGCGACATAAATTGGCGATGCAGGGGTACTACCCACCAGCACAACGTGCTGGCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCGGCGCTATCGTTTTACCGTCAATAATCTT‑TGCTTCCATCAGAGAGAGGATTCCATCTGTTCCGTAGAGTGAAGGGGAGGCTCCTTTTTTGTCACAAGCGGGGCGCGCTGTCatttgtcgttt                                <  2:432793/250‑11 (MQ=255)                                                                                                                                                                                ggcgtgcGACATGAATTTGCGATGCAGGGTTACTACCCACCAGCACAACTGCCTGTCCTGGGGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTGTTGAGCAACTGCAGAGCGCACCTGCTGCGCAATCGTTTGACCGGCAATAATCTTGTG‑TTCCATCAGAGGGAGGAGTCCAGCTGGTACGTAGATCGAAGGGGATGCGGTTATTGTGTCAGAAGCGTGTCGCGCTGTCAGGTTTGGTTTCAGAttt                         <  2:417831/247‑1 (MQ=255)                                                                                                                                                                                 cattgcGACATAAATTTGCGATTCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATATGCGCCTGAACTTTTTGAGCAACTTCAGAGGGGACCTGCTGCGCAATCGTTTTACCGTCAATAATCTT‑TGATGCCATGAGAGAGAGGGTTACATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTATCCGAGTTCGTTTCAGAtgat                        >  2:268301/3‑248 (MQ=255)                                                                                                                                                                                             aagttGCGATGCAGGGTAACTACCCACAAGCACAACAGCCTGTCCTGGTGCCCAGAGTCCGGCTACAAAACGCGCCTAAACTTTTTGACCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTGCCGTCAATAAACTA‑TGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGCGATGCGACTATTTTGTCAGAAGCGAGTTGCGCTGGCAGGTTTCGTTTCAGATATATCGCGTGAAGCg            <  2:260703/248‑1 (MQ=255)                                                                                                                                                                                                   cGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTT‑TGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACctct      <  1:47567/251‑1 (MQ=255)                                                                                                                                                                                                       gCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTT‑TGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCg  >  1:396862/1‑251 (MQ=255)                                                                                                                                                                                                                ctacCCTCCATCACAACGGCCTGTCCTGGTGCCCTCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCAcc                                                                                                                                                              <  1:239343/87‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                              TCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTT‑TGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCG  >  NC_000913/557407‑557855

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