Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	281,053	T→C	T311A (ACG→GCG)	yagA ←	integrase core domain‑containing protein YagA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	281,053	0	T	C	97.1%	101.3 / ‑6.8	34	T311A (ACG→GCG)	yagA	integrase core domain‑containing protein YagA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (27/6);  minor base G (0/1);  total (27/7)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GAACACTTTTTAACATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGCGGATGATAAGGCCGGGAGTGCCCCACCCGGATACCC  >  NC_000913/280820‑281297                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     gAACACTTTTTAACATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCAGCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACt                                                                                                                                                                                                                                     <  2:91309/251‑1 (MQ=255)      cTTTTTAACATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCTTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCcgc                                                                                                                                                                                                                                >  2:245383/1‑251 (MQ=255)       tttttAACATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACCCCTCGTAGCTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCcgcg                                                                                                                                                                                                                               >  2:206873/1‑251 (MQ=255)             aCATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACg                                                                                                                                                                                                                         >  1:328694/1‑251 (MQ=255)               aTCCTTTACCCATGGAGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGc                                                                                                                                                                                                                       >  2:26871/1‑251 (MQ=255)                        ccctGGTGATCGACTTTTTCTTCAGTTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTGCATCTCCTTCAGCCCGCCCCGTTCTCCCCTGAACGACTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGCCCATCACCCCTTC‑CTCGTATTCCGGTGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACcgc                                                                                                                                                                                                              >  2:340730/4‑251 (MQ=255)                          aTGGTGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACcgcg                                                                                                                                                                                                             >  1:338284/1‑250 (MQ=255)                                              cAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTACATCTCCTTCAGCCCGACCCGTTCTCCCCTGAAAGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATcc                                                                                                                                                                                         >  1:407600/1‑250 (MQ=255)                                                                                       cGTAGCTCCCGTCTCCCTGCATCCCCTCCACCCCGCCCCGTTCCCCCCTGAACGCCTTGCCTGCGCTCAGCCTTCCCCCTTTCCCGCTCAGCTTTCCGCTGATACCCCCTTTCCTGCCCATCCCCCCTCC‑ATCGTATCCCGGGGGCGCTGTGTTCCCGCTGTACTGCCGCGCTGACGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAg                                                                                                                                               <  1:206873/251‑1 (MQ=255)                                                                                        gTAGCTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGt                                                                                                                                              >  1:335431/1‑251 (MQ=255)                                                                                          aGCTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCCCGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTgg                                                                                                                                            >  1:412042/1‑251 (MQ=255)                                                                                              gCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                                                                                                                                                    >  1:41402/1‑223 (MQ=255)                                                                                                cGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCAGCTGTACTGCCGCGATGACGGCTGATACCGCGAGCCCGGTAACGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGa                                                                                                                                        >  1:236624/1‑249 (MQ=255)                                                                                                 gTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTCC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAgg                                                                                                                                     >  2:241850/1‑251 (MQ=255)                                                                                                          cATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGc                                                                                                                             >  2:394672/1‑250 (MQ=255)                                                                                                             ctcCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCg                                                                                                                                                   >  2:375553/1‑225 (MQ=255)                                                                                                             ctcCTTCACCCCGCCCCTTTCCCCCCGGACCGCCTGGCCTGCCCTCAGCCTTCCCCCTTCCCCCCCCAGCTTTCCGCTGATACCCCCTTTCCTAACCACCCCCCCTCC‑ACCGTATCCCGGGGGCGCGGTGTGGCCGCTGTCCGCCCGCGCGGACGTCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCg                                                                                                                                                   <  1:375553/225‑1 (MQ=255)                                                                                                                  ccaccccGCCCCTTTCCCCCCTGACCGCTTTGCCTGCGCTCAGCCTCCCCCCTTCCACCCTCAGCTTCCCGCTGATCTCCCCTTCCCTGCCCCTCCCCCCTTC‑TTCGTATCCGGGGGGCGCTGTGTTGCCGCGGTACCGCCGCGCTGAGGGTTGATCCGGCGGCCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCAcc                                                                                                                     <  2:39766/246‑1 (MQ=255)                                                                                                                                           ccccTTGCCTGCGCTCAGACTCCCCCCTTCCACGCTCAGCTTCCCGCTAATATCCACTTTCCTGACCATCCCCCCTTC‑CTCGTATTCCGGGGGCGCTGTGTTGCCGCTTTACGGCCGCGCGGGCGGCTGATCCCGCGAGCCGGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTgca                                                                                           <  1:332484/249‑1 (MQ=255)                                                                                                                                           cGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTcccagc                                                                                            >  1:150699/1‑247 (MQ=255)                                                                                                                                                   cTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTcc                                                                                   >  1:350117/1‑251 (MQ=255)                                                                                                                                                       gcTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCCCTTCCGCCt                                                                               >  2:222938/1‑251 (MQ=255)                                                                                                                                                                                 gCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGc                                                     >  2:96103/1‑251 (MQ=255)                                                                                                                                                                                        gatggtatCCTCTTTCCTGCCCACCACCCCTCA‑ATCGTATGCGGGGGGCGGTGTTTCGCCGTTGATCTGCCGCGCTGACGGCTGATACCTCGAGCCGGTTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCt                                              <  1:239147/246‑1 (MQ=255)                                                                                                                                                                                        gCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCt                                              >  2:48228/1‑251 (MQ=255)                                                                                                                                                                                          tGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGc                                            >  1:271725/1‑251 (MQ=255)                                                                                                                                                                                           ttatccccatttccctcccctccccccTTCCATCGTTTTCCGGGGGCCCTGTTTTGCCGTTGTCCCGCCGCGCTGACGGCTAATCCCGGGACCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCt                                              <  1:239177/226‑1 (MQ=255)                                                                                                                                                                                             tatCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTc                                         >  1:364344/1‑251 (MQ=255)                                                                                                                                                                                             tatCCACTTTACTGACCATCCCCCCTTC‑ATCGTCTTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCt                                                                                                              >  2:394032/1‑182 (MQ=255)                                                                                                                                                                                                 cACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGCg                                     >  1:387322/1‑251 (MQ=255)                                                                                                                                                                                                    tttCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGCGgat                                  >  2:30512/1‑251 (MQ=255)                                                                                                                                                                                                                 aCCCCTTC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGCGGATGATAAGGCCggg                      >  1:100156/1‑250 (MQ=255)                                                                                                                                                                                                                      ttC‑ATCGTATTCCGGGGGCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGCGGATGATAAGGCCGGGAGTGcc                >  1:412679/1‑251 (MQ=255)                                                                                                                                                                                                                                     ggggCGCTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGCGGATGATAAGGCCGGGAGTGCCCCACCCGGATAccc  >  1:384330/1‑251 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     GAACACTTTTTAACATCCTTTACCCATGGTGATCGACTTTTTCTTCAGGTCGATCACCCCCACTTTCGTGCTGTACCACCACACCTCGTAGCTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTC‑ATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGCGGATGATAAGGCCGGGAGTGCCCCACCCGGATACCC  >  NC_000913/280820‑281297

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