Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,438,793	C→T	R282W (CGG→TGG)	flk →	putative flagella assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,438,793	0	C	T	100.0%	201.7 / NA	57	R282W (CGG→TGG)	flk	putative flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/27);  total (30/27)

AGATGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGAT  >  NC_000913/2438656‑2438923                                                                                                                                          |                                                                                                                                   aGATGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAAtgg                                                                                                                                  >  1:660998/1‑140 (MQ=255)   aTGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAAtggcg                                                                                                                                >  2:752774/1‑140 (MQ=255)   aTGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAAtggcg                                                                                                                                >  1:93337/1‑140 (MQ=255)     gACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCggg                                                                                                                              >  2:770783/1‑140 (MQ=255)        gCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTa                                                                                                                           >  1:209754/1‑140 (MQ=255)        gCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTa                                                                                                                           >  2:60592/1‑140 (MQ=255)         cTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTaa                                                                                                                          <  1:770783/140‑1 (MQ=255)           aaaaCAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAAGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAg                                                                                                                        >  1:771748/1‑140 (MQ=255)             aacaacCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAAtggc                                                                                                                                 >  1:349643/1‑129 (MQ=255)             aacaacCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAAtggc                                                                                                                                 <  2:349643/129‑1 (MQ=255)             aacaacCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAggg                                                                                                                      <  2:767384/140‑1 (MQ=255)             aacaacCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAggg                                                                                                                      >  1:767384/1‑140 (MQ=255)               caacCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGAt                                                                                                                    <  2:660998/140‑1 (MQ=255)               caacCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGAt                                                                                                                    <  1:723723/140‑1 (MQ=255)                aacCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGAtt                                                                                                                   <  2:209754/140‑1 (MQ=255)                  ccTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGa                                                                                                                 >  1:433658/1‑140 (MQ=255)                     ttAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATcc                                                                                                              <  2:882111/140‑1 (MQ=255)                     ttAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAAGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATcc                                                                                                              <  2:771748/140‑1 (MQ=255)                      tAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCcg                                                                                                             >  2:634457/1‑140 (MQ=255)                            ccAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGtt                                                                                                       >  1:345042/1‑140 (MQ=255)                            ccAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGtt                                                                                                       >  2:579873/1‑140 (MQ=255)                                tGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGCTATGc                                                                                                   >  1:677465/1‑140 (MQ=255)                                     tGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCg                                                                                              <  1:398861/140‑1 (MQ=255)                                     tGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCg                                                                                              <  1:634457/140‑1 (MQ=255)                                             ttATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATACAACCGCTGCTTaa                                                                                      >  2:295270/1‑140 (MQ=255)                                                    cATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCtttt                                                                               <  2:93337/140‑1 (MQ=255)                                                       atatCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCa                                                                            >  2:53298/1‑140 (MQ=255)                                                          tCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCg                                                                         <  1:713875/140‑1 (MQ=255)                                                          tCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCg                                                                         <  1:53298/140‑1 (MQ=255)                                                          tCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCg                                                                         >  1:61526/1‑140 (MQ=255)                                                          tCCAGCAAAAATATGGTCTTTCTGCGCAATCATCACTTTCTTCTGCCAAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTaa                                                                                      >  1:744701/1‑127 (MQ=255)                                                          tCCAGCAAAAATATGGTCTTTCTGCGCAATCATCACTTTCTTCTGCCAAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTaa                                                                                      <  2:744701/127‑1 (MQ=255)                                                             aGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAAGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCgatg                                                                      >  1:754071/1‑140 (MQ=255)                                                                aaaaaTATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGatgatg                                                                   <  1:947320/140‑1 (MQ=255)                                                                    atatGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGaca                                                               <  1:579873/140‑1 (MQ=255)                                                                       tGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGt                                                            >  2:875388/1‑140 (MQ=255)                                                                        ggTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAAGGATTTATCCGCGGGGTATGCGAGCGCTGCTTTATTCTTTTCCATCGAAGATAGGCACGtt                                                           >  2:212075/1‑140 (MQ=255)                                                                         gTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTg                                                          <  2:754071/140‑1 (MQ=255)                                                                         gTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTg                                                          >  1:952298/1‑140 (MQ=255)                                                                              tCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCaaaa                                                     >  2:1002561/1‑140 (MQ=255)                                                                                  cgcAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATg                                                 >  1:838059/1‑140 (MQ=255)                                                                                    cAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTg                                                          >  2:823889/1‑129 (MQ=255)                                                                                    cAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTg                                                          <  1:823889/129‑1 (MQ=255)                                                                                       tcatcGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGGTTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAAc                                            >  2:722143/1‑140 (MQ=255)                                                                                            gCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTc                                       <  2:61526/140‑1 (MQ=255)                                                                                              tttCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTccc                                     <  2:952298/140‑1 (MQ=255)                                                                                                   tctGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCt                                <  1:875388/140‑1 (MQ=255)                                                                                                   tctGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCt                                <  2:433658/140‑1 (MQ=255)                                                                                                   tctGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCt                                <  1:1002561/140‑1 (MQ=255)                                                                                                     tGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCcgc                                   <  2:655188/135‑1 (MQ=255)                                                                                                     tGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCcgc                                   >  1:655188/1‑135 (MQ=255)                                                                                                          aGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTTATCCGCGTGTTTTGCAACCGCTGCTTAATCCTTTTTCCCCGATGATGGACCCGTTGCAAAATATGGCAACGCGTCCCGCGCTCTGGGTa                         >  2:523530/1‑140 (MQ=255)                                                                                                          aGGCCGAGGATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATa                         >  2:953143/1‑140 (MQ=255)                                                                                                                 ggATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTgg                            >  2:755977/1‑130 (MQ=255)                                                                                                                 ggATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTgg                            <  1:755977/130‑1 (MQ=255)                                                                                                                 ggATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAg                  >  2:36655/1‑140 (MQ=255)                                                                                                                 ggATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAg                  >  2:778109/1‑140 (MQ=255)                                                                                                                 ggATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAg                  >  1:96664/1‑140 (MQ=255)                                                                                                                 ggATATTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGAGTGATCCGCGTGTTATGGAACCGCTGCTTAATCCTTTTCCACCGATGATGGGCACGGTTCAAAATATAGCAACGCGTCCCGCGCTGTTGATACTTTTaa                  >  2:875258/1‑139 (MQ=255)                                                                                                                    tatTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTcg               <  2:345042/140‑1 (MQ=255)                                                                                                                    tatTCTTAATCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTcg               <  2:107013/140‑1 (MQ=255)                                                                                                                             tCAGCTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATcc      <  1:60592/140‑1 (MQ=255)                                                                                                                                 cTTTATCAATGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCACTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGAt  >  1:479048/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                   AGATGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGAT  >  NC_000913/2438656‑2438923

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