Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,019,460	G→A	V359I (GTC→ATC)	fliK →	flagellar hook‑length control protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,019,460	0	G	A	100.0%	199.1 / NA	62	V359I (GTC→ATC)	fliK	flagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/29);  total (33/29)

AAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGACACAGG  >  NC_000913/2019314‑2019604                                                                                                                                                   |                                                                                                                                                 aacgTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGT‑‑‑CAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATctct                                                                                                                                             <  5:66684/146‑1 (MQ=255)          ttCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTAc                                                                                                                                         >  3:142905/1‑147 (MQ=255)          ttCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACaa                                                                                                                                       <  4:316569/149‑1 (MQ=255)              gTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGc                                                                                                                                   >  4:81632/1‑149 (MQ=255)                   ggCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGACGCTTCCCAGCAACAGCAAAACCAACGCACAGCAAACCATCAACCTCTTGCGGCGGAAGACGACGAAACGCATCCCGGTCCCATCACAAAACAAGGGCatga                                                                                                                               >  6:236906/1‑147 (MQ=255)                     cAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAAc                                                                                                                            <  3:284632/149‑1 (MQ=255)                     cAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAAc                                                                                                                            <  4:142905/149‑1 (MQ=255)                      aaaGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACACCAAACCATGAACCTCTGGCGGGGGAAGACGACGCTACGCTTCCGGTTCCCATCTCAGTACAAGGGCGTGTaa                                                                                                                             >  8:319359/1‑147 (MQ=255)                       aaGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACCGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAg                                                                                                                          >  2:405642/1‑149 (MQ=255)                           aaCATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCaa                                                                                                                      >  3:4642/1‑149 (MQ=255)                             cATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCaa                                                                                                                      >  2:6783/1‑147 (MQ=255)                                  gTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAAc                                                                                                                            >  7:263999/1‑136 (MQ=255)                                  gTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAAc                                                                                                                            <  8:263999/136‑1 (MQ=255)                                   tGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACa                                                                                                                           <  4:305590/136‑1 (MQ=255)                                   tGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACa                                                                                                                           >  3:305590/1‑136 (MQ=255)                                      cGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAGCCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGtt                                                                                                            >  2:38807/1‑148 (MQ=255)                                       gAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCgtgt                                                                                                                               <  4:86870/128‑1 (MQ=255)                                       gAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCgtgt                                                                                                                               >  3:86870/1‑128 (MQ=255)                                       gAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGa                                                                                                          >  6:115853/1‑149 (MQ=255)                                             tttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtttt                                                                                                    >  2:173233/1‑149 (MQ=255)                                             tttAGTGGTCAGCAGCAGGCCGATTCCAAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAccg                                                                                                                 >  7:234761/1‑136 (MQ=255)                                              ttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                                                                                                   >  2:44555/1‑149 (MQ=255)                                               tAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                                                                                                   <  3:81632/148‑1 (MQ=255)                                                  tGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGcc                                                                                                >  1:162226/1‑148 (MQ=255)                                                              gcccGCTTCCCAGCAACAGCAAAGACAACGCACAGCAAACCATGACCCTCTGGCGGGGGACGACGACGAGACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGGCAGAGGt                                                                                   <  6:234337/147‑1 (MQ=255)                                                                       ccAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCgtgt                                                                                                                               <  6:222549/96‑1 (MQ=255)                                                                       ccAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCgtgt                                                                                                                               >  5:222549/1‑96 (MQ=255)                                                                       ccAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGGGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTa                                                                          <  5:131127/149‑1 (MQ=255)                                                                                  aaaGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCtt                                                               <  1:173233/149‑1 (MQ=255)                                                                                            cacaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCtt                                                     <  7:323165/149‑1 (MQ=255)                                                                                               aGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCATCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGAGGCGATGATATTTTCGCCTACCGACAGAGGTAGCACCGAAATCCGCGTCTTTTCCCCGCTttgt                                                  >  2:78205/1‑149 (MQ=255)                                                                                                gCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTttgtt                                                 <  1:44555/149‑1 (MQ=255)                                                                                                gCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTttgtt                                                 <  2:162226/149‑1 (MQ=255)                                                                                                    aCCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTg                                                                                                           <  8:305307/87‑1 (MQ=255)                                                                                                    aCCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTg                                                                                                           >  7:305307/1‑87 (MQ=255)                                                                                                     ccATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCttt                                                              <  8:250422/131‑1 (MQ=255)                                                                                                     ccATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCttt                                                              >  7:250422/1‑131 (MQ=255)                                                                                                      cATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTc                                           <  1:6783/149‑1 (MQ=255)                                                                                                           aCCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGAc                                      <  1:405642/149‑1 (MQ=255)                                                                                                             ctctGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                                                                                                   <  2:219057/86‑1 (MQ=255)                                                                                                             ctctGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                                                                                                   >  1:219057/1‑86 (MQ=255)                                                                                                             ctctGGCGGGGGAAGACGACGATACGCTTCCCGTTACCATCTCTTTACAAGGGCGTGTAACAGGCTAGAGCGGGGTTGATATTTTCGCCTTACGTCAGAGGTAGAACCGTAATCCGCCTCTTTTCCCCGCTTTGTAGCGATCAAGACGc                                    >  6:283411/1‑149 (MQ=255)                                                                                                                 ggCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGa                                >  5:182395/1‑149 (MQ=255)                                                                                                                  gCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCa                                   >  8:168240/1‑145 (MQ=255)                                                                                                                  gCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCa                                   <  7:168240/145‑1 (MQ=255)                                                                                                                  gCGGGGGAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAgg                                 <  8:139415/147‑1 (MQ=255)                                                                                                                      gggAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                                                                                                   >  6:207960/1‑77 (MQ=255)                                                                                                                      gggAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                                                                                                   <  5:207960/77‑1 (MQ=255)                                                                                                                      gggAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAAt                            >  5:105053/1‑148 (MQ=255)                                                                                                                       ggAAGACGACGATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTa                          <  5:115853/149‑1 (MQ=255)                                                                                                                            acgacgATACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGa                     <  4:75759/149‑1 (MQ=255)                                                                                                                               aagaTACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGACGATaa                  <  4:4642/147‑1 (MQ=255)                                                                                                                                cgaTACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCag                                                                                       >  1:114704/1‑79 (MQ=255)                                                                                                                                cgaTACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCag                                                                                       <  2:114704/79‑1 (MQ=255)                                                                                                                                  aTACGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAg                 >  3:173620/1‑147 (MQ=255)                                                                                                                                    aCGCTTCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAg              >  3:232237/1‑148 (MQ=255)                                                                                                                                        ttCCGGTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCGTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGc          >  8:175717/1‑148 (MQ=255)                                                                                                                                             gTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAg           >  6:63352/1‑142 (MQ=255)                                                                                                                                             gTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAg           <  5:63352/142‑1 (MQ=255)                                                                                                                                             gTTCCCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGaca      >  5:273512/1‑147 (MQ=255)                                                                                                                                             gTTCCCATCTCTTTACAAGCGCGTGTAACAGGCAACATCGGCGTTGATATTTTCGCCTATCGTCAGACGTAGCACCTTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGaca      >  6:300796/1‑147 (MQ=255)                                                                                                                                                cccATCTCTTTACAAGGGCGTGTAACAGGCATCAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGACACAgg  <  7:319359/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 AAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGACACAGG  >  NC_000913/2019314‑2019604

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