Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,894,306	T→G	V119G (GTG→GGG)	yqcE →	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,894,306	0	T	G	98.4%	216.7 / NA	63	V119G (GTG→GGG)	yqcE	MFS transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (35/27);  minor base C (0/1);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTCCGCGATGATCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTG  >  NZ_CP009273/2894185‑2894437                                                                                                                          |                                                                                                                                    aCTTCCGCGATGATCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCa                                                                                                                     <  2:1275156/138‑1 (MQ=255)     ccGCGATGATCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCg                                                                                                                  <  1:623118/137‑1 (MQ=255)      cgcgATGATCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTg                                                                                                               >  1:1188075/1‑139 (MQ=255)             atCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGc                                                                                                         >  2:183779/1‑138 (MQ=255)                  tACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCa                                                                                                    >  2:1161654/1‑138 (MQ=255)                     cGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGc                                                                                                >  1:894249/1‑139 (MQ=255)                       taTTACTGGGTCTGTTAATGGCAACGAGTCCACCGCTATGGGTACGGCTCTGTATTCATATCGCCTTTGCGATAACGACGAGTTTATTGCAGAGGTCTGGGTCGCTTAAAGCTGCATCGTTGCTAGGCGTTCATAGc                                                                                                >  2:1135138/2‑137 (MQ=255)                       gATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGc                                                                                                >  2:622753/1‑137 (MQ=255)                               ggTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGa                                                                                      <  2:898598/139‑1 (MQ=255)                                ggtcgTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGa                                                                                      <  2:591818/134‑1 (MQ=255)                                 tCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTgg                                                                                                          >  1:610950/1‑117 (MQ=255)                                 tCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTgg                                                                                                          <  2:610950/117‑1 (MQ=255)                                       aaTGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTAt                                                                               <  1:1187099/138‑1 (MQ=255)                                         tGGCAACGTATGCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATggg                                                                            <  2:995186/139‑1 (MQ=255)                                         tGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAAtt                                                                                 <  2:1043452/134‑1 (MQ=255)                                         tGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAAtt                                                                                 >  1:1043452/1‑134 (MQ=255)                                         tGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATggg                                                                            <  1:543936/139‑1 (MQ=255)                                          ggCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGc                                                                           >  2:62800/1‑139 (MQ=255)                                           gaaaCGTATCCCCCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGc                                                                           <  1:1135138/136‑1 (MQ=255)                                               cGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCtggat                                                                      >  2:1143331/1‑139 (MQ=255)                                               cGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCtggat                                                                      >  1:101087/1‑139 (MQ=255)                                                          cTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTgc                                                           >  1:148139/1‑139 (MQ=255)                                                          cTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTgc                                                           >  1:1164947/1‑139 (MQ=255)                                                             tgGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTgcg                                                          <  2:335597/137‑1 (MQ=255)                                                               ggTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTgcg                                                          >  1:375234/1‑135 (MQ=255)                                                               ggTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTgcg                                                          <  2:375234/135‑1 (MQ=255)                                                               ggTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTATGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCgcgg                                                       <  1:376355/138‑1 (MQ=255)                                                                 tAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGt                                                    <  1:622753/139‑1 (MQ=255)                                                                 tAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGt                                                    <  2:1164947/139‑1 (MQ=255)                                                                 tAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGt                                                    <  1:62800/139‑1 (MQ=255)                                                                 tAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAAAGCTGTGGTCGGCGTCGATTAAAGCTGCATCGTTACTTGGCGAGCATAGCGAACAAGGAAAAATTATGGGCTTGATGGAAGTGCTGCGCGGCGt                                                    <  2:449894/139‑1 (MQ=255)                                                                   aTGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCg                                                  <  1:995851/139‑1 (MQ=255)                                                                         tGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTaa                                            >  1:1128278/1‑139 (MQ=255)                                                                         tGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTaa                                            >  1:731260/1‑139 (MQ=255)                                                                         tGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTaa                                            >  1:655959/1‑139 (MQ=255)                                                                           tATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATg                                          >  2:398944/1‑139 (MQ=255)                                                                                aGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCt                                     >  1:1075785/1‑139 (MQ=255)                                                                                  aTCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTgg                                   <  1:956463/139‑1 (MQ=255)                                                                                  aTCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTgg                                   <  2:115224/139‑1 (MQ=255)                                                                                  aTCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTgg                                   >  2:502163/1‑139 (MQ=255)                                                                                     gCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCgg                                >  1:1245337/1‑139 (MQ=255)                                                                                     gCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCgg                                >  2:511736/1‑139 (MQ=255)                                                                                     gCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCgg                                >  2:831260/1‑139 (MQ=255)                                                                                     gCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCgg                                >  2:59470/1‑139 (MQ=255)                                                                                     gCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCgg                                >  1:1108350/1‑139 (MQ=255)                                                                                     gCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCgg                                >  2:650564/1‑139 (MQ=255)                                                                                     gCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCgg                                >  1:280422/1‑139 (MQ=255)                                                                                     gCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCgg                                >  1:217027/1‑139 (MQ=255)                                                                                       cTTTGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGgtg                              >  1:108192/1‑139 (MQ=255)                                                                                          tGCGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGtt                            >  2:97698/1‑138 (MQ=255)                                                                                            cGATAACGACGATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTa                          <  1:242009/138‑1 (MQ=255)                                                                                                 acgacgATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCAtgt                    <  1:650564/139‑1 (MQ=255)                                                                                                 acgacgATTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCAtg                     <  2:1188075/138‑1 (MQ=255)                                                                                                     cgaTTTTAATGCTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGt                >  1:983895/1‑139 (MQ=255)                                                                                                                cTGTGGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCt     <  2:1199643/139‑1 (MQ=255)                                                                                                                  gtgGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTa                                                                                >  2:192397/1‑62 (MQ=255)                                                                                                                  gtgGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTa                                                                                <  1:192397/62‑1 (MQ=255)                                                                                                                  gtgGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCttt   >  2:590716/1‑139 (MQ=255)                                                                                                                  gtgGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCttt   >  1:1001838/1‑139 (MQ=255)                                                                                                                   tgGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTgg                                   <  1:215647/106‑1 (MQ=255)                                                                                                                   tgGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTgg                                   >  2:215647/1‑106 (MQ=255)                                                                                                                   tgGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCttta  <  2:1245337/139‑2 (MQ=255)                                                                                                                   tgGTCGGGGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTg  <  2:731260/139‑1 (MQ=255)                                                                                                                          |                                                                                                                                    ACTTCCGCGATGATCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTG  >  NZ_CP009273/2894185‑2894437

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