Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	932,340	(A)8→7	intergenic (‑290/‑255)	trxB ← / → lrp	thioredoxin reductase/DNA‑binding transcriptional dual regulator Lrp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	932,333	0	A	.	100.0%	226.4 / NA	53	intergenic (‑283/‑262)	trxB/lrp	thioredoxin reductase/DNA‑binding transcriptional dual regulator Lrp
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (26/27);  total (26/27)

TGCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCTC  >  NC_000913/932194‑932466                                                                                                                                            |                                                                                                                                      tGCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGgaaga                                                                                                                                       >  1:539262/1‑140 (MQ=255) tGCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGgaaga                                                                                                                                       >  1:439845/1‑140 (MQ=255)   cTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaa                                                                                                                                     >  1:114464/1‑140 (MQ=255)   cTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaa                                                                                                                                     >  1:134318/1‑140 (MQ=255)      cTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaa                                                                                                                                  >  2:337352/1‑140 (MQ=255)       ttAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaa                                                                                                                                 <  2:439845/140‑1 (MQ=255)        tAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaac                                                                                                                                >  2:470241/1‑139 (MQ=255)        tAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaac                                                                                                                                >  2:248591/1‑139 (MQ=255)         aaaTTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTGGAGGAAGaaaaaaaca                                                                                                                               <  1:106393/140‑3 (MQ=255)           aTTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGaaaaaaacagt                                                                                                                             <  2:491305/140‑5 (MQ=255)            tttCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTa                                                                                                                           <  1:782365/140‑1 (MQ=255)                ccGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCt                                                                                                                       >  2:267021/1‑140 (MQ=255)                   cTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtat                                                                                                                    <  2:518596/140‑1 (MQ=255)                    tttATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTtata                                                                                                                   <  2:114464/140‑1 (MQ=255)                              gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTatat                                                                                                                  >  2:558322/1‑131 (MQ=255)                              gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTatat                                                                                                                  <  1:558322/131‑1 (MQ=255)                              gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                         >  1:538204/1‑140 (MQ=255)                              gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                         >  1:664762/1‑140 (MQ=255)                              gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                         >  1:359835/1‑140 (MQ=255)                              gATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCAATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAAc                                                                                                         >  2:560109/1‑140 (MQ=255)                                ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACca                                                                                                       >  1:447055/1‑140 (MQ=255)                                ttAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACca                                                                                                       <  1:344508/140‑1 (MQ=255)                                         gaATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTaa                                                                                              >  2:79064/1‑140 (MQ=255)                                               aCGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCa                                                                                        >  1:369598/1‑140 (MQ=255)                                                  ccccTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTCTTCTTATATGCGCATAACCATGCATGTAAATACCATGt                                                                                     >  2:679164/1‑140 (MQ=255)                                                    ccTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGttt                                                                                   >  1:11894/1‑140 (MQ=255)                                                    ccTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGttt                                                                                   >  2:317718/1‑140 (MQ=255)                                                     cTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTCCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa                                                                                  <  2:563605/140‑1 (MQ=255)                                                          aaTGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTg                                                                             <  2:93871/140‑1 (MQ=255)                                                               atatCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGt                                                                        <  1:79064/140‑1 (MQ=255)                                                               atatCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGt                                                                        >  1:782109/1‑140 (MQ=255)                                                                  tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa                                                                                  >  2:655251/1‑127 (MQ=255)                                                                  tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa                                                                                  <  1:655251/127‑1 (MQ=255)                                                                  tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGaa                                                                     >  1:689279/1‑140 (MQ=255)                                                                  tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGaa                                                                     >  1:96690/1‑140 (MQ=255)                                                                  tCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGaa                                                                     >  2:443787/1‑140 (MQ=255)                                                                   cTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGaaa                                                                    <  1:360436/140‑1 (MQ=255)                                                                   cTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGaaa                                                                    <  1:317718/140‑1 (MQ=255)                                                                     ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa                                                                                  <  2:727472/124‑1 (MQ=255)                                                                     ggCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTa                                                                                  >  1:727472/1‑124 (MQ=255)                                                                       cATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTa                                                                >  1:646056/1‑140 (MQ=255)                                                                        aTGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTAc                                                               <  1:184969/140‑1 (MQ=255)                                                                          gttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                             <  1:679164/140‑1 (MQ=255)                                                                          gttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                             <  2:96690/140‑1 (MQ=255)                                                                          gttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                             <  2:447055/140‑1 (MQ=255)                                                                          gttgtACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGt                                                             <  2:134318/140‑1 (MQ=255)                                                                               aCTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGc                                                        <  1:277549/140‑1 (MQ=255)                                                                                cTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCg                                                       >  2:567258/1‑140 (MQ=255)                                                                                cTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTAACGTGCTAGTGAAATCTACGTATGGCg                                                       <  1:102538/140‑1 (MQ=255)                                                                                cTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTAACGTGCTAGTGAAATCTACGTATGGCg                                                       >  2:102538/1‑140 (MQ=255)                                                                                      aTCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACa                                                 <  2:689279/140‑1 (MQ=255)                                                                                        cGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGa                                               >  2:360861/1‑140 (MQ=255)                                                                                          aTGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACg                                             <  2:539262/140‑1 (MQ=255)                                                                                                tGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTc                                       >  2:119705/1‑140 (MQ=255)                                                                                                     tGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGAt                                  >  2:809881/1‑140 (MQ=255)                                                                                                            ccccTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCACGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATa                           >  2:338454/1‑140 (MQ=255)                                                                                                              ccTGGTGTTTTGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGc                         <  1:267021/140‑1 (MQ=255)                                                                                                                       ttGCGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAg                >  2:626709/1‑140 (MQ=255)                                                                                                                          cGAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCa             <  2:664762/140‑1 (MQ=255)                                                                                                                           gAAAACATTCGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCaa            <  1:360861/140‑1 (MQ=255)                                                                                                                                    cGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGtcttct   <  2:552716/140‑1 (MQ=255)                                                                                                                                    cGAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGtcttct   <  2:11894/140‑1 (MQ=255)                                                                                                                                     gAGGAAG‑AAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTtctc  <  1:119705/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                      TGCTACTTAAATTTCCCGCTTTATAAGCCGATTAAATGATGAATAAACGCCCCTGTTAATGAATATCTGGCATGTTGTACTAAAAATCGATGTTTTGCTTTGACAATCCCCTGGTGTTTTGCGAAAACATTCGAGGAAGAAAAAAAACAGTATTCTTATATGCGCATAACCATGCATGTAAATACCATGTTTACCGTGCTAGTGAAATCTACGTATGGCGTGGACAGACGCCATTCGTGATGTCGATAGCTGCCACAAGGCAACGGTCTTCTC  >  NC_000913/932194‑932466

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