Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,564	0	A	C	50.0%	80.1 / 20.6	58	*559E (TAA→GAA)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (26/3);  new base C (0/29);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 3.30e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATA  >  NC_000913/501421‑501706                                                                                                                                                 |                                                                                                                                               gATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTc                                                                                                                                         >  1:966138/1‑151 (MQ=255)        aaGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCaccacc                                                                                                                                   >  2:333860/1‑150 (MQ=255)        aaGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCt                                                                                                                                  >  1:747652/1‑151 (MQ=255)               cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACccgc                                                                                                                           >  2:547490/1‑151 (MQ=255)                           ggCCAGGTAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCCCCGGTGGCGTTTc                                                                                                               >  2:963890/1‑151 (MQ=255)                               aGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                                                                           >  2:351460/1‑151 (MQ=255)                                  tAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGt                                                                                                        >  2:1002744/1‑151 (MQ=255)                                         tccccgccgcATCTGGCTTTTGCCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATCCCTACATTTTTTTCATATTTTACACCCGGCAC‑CCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCa                                                                                                    <  1:182538/144‑1 (MQ=37)                                            ccgccgcATTGGGTATTTGCCTCCATCTGCCTGATGCAGCGCGGCCGCATCTTATCATGCCACCATTTTTTTTAT‑TTTTACCTCCGGGAACCCACGTTTTCCCCCGTCCCCCCCACACCCGCCGGTGGCGTTTCCAGCAGTTcc                                                                                                     <  2:520907/143‑1 (MQ=255)                                            aCGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                              >  2:241986/1‑151 (MQ=255)                                            aCGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                              >  1:811759/1‑151 (MQ=255)                                               ccgttgcTGGCATTTCTCTCTCGATACCGGTTGCGACGGTGGCGGGTCTTAACAAGCCTCCATATTTTACATATTTTTCATCCGGCAA‑CCTCCTTTTCCCCCGTCCCCCCCTCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                             <  1:1198479/143‑1 (MQ=255)                                               ccgcATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCTTCTTATCATGCCTACTTTTTTTTCATTTTTTACATCCGGCAA‑CCCCCTTTTCCCCCGTCACCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTAc                                                                                           <  1:447647/151‑1 (MQ=38)                                                  cATCTGTCTTTTGTCCCTCGTCGCCGGAGGCGCCGCCGGCGCGTTTTCTCGTCCCTAATTTTTTTTCATTTTTTTCATCCGCCAA‑CCACTGTTGCCCCTGCCGCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                              <  2:1201909/145‑1 (MQ=255)                                                    tCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTCCACCCGGCAA‑CCCCCTTTTACCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                                                                                        <  2:1170954/149‑1 (MQ=255)                                                     cTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTAACCCGTCACCACCTCACCCGCCGGCGGCGTTTCCAGCAGATCCAGCATGGTACGGGCg                                                                                      >  2:792507/1‑150 (MQ=255)                                                      tGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAt                                                                                    >  1:362765/1‑151 (MQ=255)                                                       ggCTTTTGGCTCTCTATCCTGGATGCGACGCTGGCCCGTTTTATCATGCCTACTTTTTTTCCATTTTTTACACCCGCCAC‑CCCCCGTTACCCCCTCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAt                                                                                    <  2:79697/150‑1 (MQ=255)                                                       ggCTTTTGGCTCTCGATGCCTGATGCGACGCGGGCGCTTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCAC‑CCACCTTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAt                                                                                    <  1:845602/150‑1 (MQ=255)                                                       ggCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt                                                                                   >  1:1224315/1‑151 (MQ=255)                                                           tttGTCTCCCGATGCCTTACGCGACGCGGGCGCGTCTTATCTGCCCTACATATTTTTAATATTTTACACCCGCCAA‑CCACCGTTTCCCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAttt                                                                                  <  1:1025726/148‑1 (MQ=255)                                                           tttGGCTCTCGTTGCCTGATGCGCCGCGGGCGCGTCTTATCATGCCTACATTTTTTTCATTTTTTCCACCCGGCAC‑CCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCa                                                                                <  1:125289/150‑1 (MQ=255)                                                           tttGGCTCTCGATGCCTGATGCGCCGCGGGCGCTTCTTTTCATCCCTACATTTTTTTCATTTTTTACACCCGGCAC‑CCCCCTTTTCCCCCGCCCCCCCCGCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc                                                                                 <  2:11358/149‑1 (MQ=255)                                                                   tcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc                                                                       >  2:963978/1‑151 (MQ=255)                                                                         ccTGATGCGACGCGGGCGCGTTTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCCCCTTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                  <  1:345796/150‑1 (MQ=255)                                                                           tGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                               >  1:893082/1‑151 (MQ=255)                                                                               gCGACGCGGGGGGTTTTTATCATGCCTACATATTTTTCATATTTTACACCCGGCAC‑CCCCCGTTTCCCCCGCCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTg                                                            <  2:793154/150‑1 (MQ=255)                                                                                 gACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                         >  1:1128184/1‑151 (MQ=255)                                                                                       ggCGCTTTTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCAA‑CCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGc                                                      <  1:541485/148‑1 (MQ=255)                                                                                       ggCGCGTTTTATCATCCCTACTTTTTTTTAATTTTTTACACCCGCCAC‑CCCCCTTTTACCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcac                                                    <  2:730483/150‑1 (MQ=255)                                                                                         cgctttttATCATGCCTACATATTTTTAATTTTTTCCATCCGCCAA‑CCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacca                                                  <  1:582016/144‑1 (MQ=255)                                                                                                tATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                          >  1:743163/1‑151 (MQ=255)                                                                                                tATCATGCCTACATATTTTTCATATTTTACACCCGGCAC‑CCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTc                                             <  1:862260/148‑1 (MQ=255)                                                                                                tATCATGCCTACATATTTTTCATATTTTACACCCGGCAC‑CCCCCGTTTCCCCCGTCCCCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTc                                             <  2:468504/148‑1 (MQ=255)                                                                                                 aTCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                         >  1:39453/1‑151 (MQ=255)                                                                                                  tCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                         >  1:8469/1‑150 (MQ=255)                                                                                                  tCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                         >  1:534636/1‑150 (MQ=255)                                                                                                   cATGCCTACATTTTTTTCATATTTTCCACCCGGCAC‑CCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                          <  1:918004/148‑1 (MQ=255)                                                                                                       ccTACATTTTTTTCATTTTTTACATCCGGCAC‑CCACCTTTTACCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata                                      <  2:656535/148‑1 (MQ=255)                                                                                                          aCATTTTTTTCATATTTTACACCCGACAA‑CCCCCGTTTCCCCCTTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                   <  1:364875/148‑1 (MQ=255)                                                                                                           cATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                >  2:582377/1‑150 (MQ=255)                                                                                                               tttttCATTTTTTACCTCGGGCAA‑CCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                                  <  1:911214/144‑1 (MQ=255)                                                                                                               tttttCATATTTTCCATCCGGCAA‑CCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                <  1:675793/146‑1 (MQ=255)                                                                                                                  ttctttttttccccccGGCAC‑CCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCAATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCAt                         <  2:188934/136‑1 (MQ=255)                                                                                                                     atatTTTACATCCGCCAC‑CCCCCTTTCCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt                      <  1:1092294/150‑1 (MQ=255)                                                                                                                      tatTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc                     >  1:537697/1‑150 (MQ=255)                                                                                                                        ttttACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                   >  2:881282/1‑150 (MQ=255)                                                                                                                        tttaccacccGCCAC‑CCACCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc                     <  2:11179/140‑1 (MQ=255)                                                                                                                         tttACCCCCGGCAC‑CCCCCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                  <  2:288781/150‑1 (MQ=255)                                                                                                                           tACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                >  2:174722/1‑150 (MQ=255)                                                                                                                           tACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg               >  1:730805/1‑151 (MQ=255)                                                                                                                           tACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg               >  2:278265/1‑151 (MQ=255)                                                                                                                                ccGGCAC‑CCCCCTTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc          <  1:915026/151‑1 (MQ=255)                                                                                                                                ccGCCAC‑CCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg           <  1:327015/150‑1 (MQ=255)                                                                                                                                 cGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg         >  1:539700/1‑151 (MQ=255)                                                                                                                                     caccccccttttCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg         <  2:741004/139‑1 (MQ=255)                                                                                                                                      accccccttttCCCCCATCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGc      <  2:951021/142‑1 (MQ=255)                                                                                                                                         ccccGTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata  <  2:202919/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                               GATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATA  >  NC_000913/501421‑501706

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