Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,438,548	T→C	V200A (GTG→GCG)	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,548	0	T	C	100.0%	262.6 / NA	74	V200A (GTG→GCG)	flk	putative flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (41/33);  total (41/33)

CGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGTGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCTTTAGA  >  NC_000913/2438418‑2438680                                                                                                                                   |                                                                                                                                     cgcgAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATg                                                                                                                              <  2:1434933/139‑1 (MQ=255)  gcgAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                             >  2:425547/1‑139 (MQ=255)  gcgAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                             >  1:367871/1‑139 (MQ=255)  gcgAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                             <  1:1443076/139‑1 (MQ=255)          accgaccgTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTaa                                                                                                                     >  2:1627437/1‑139 (MQ=255)           ccgaccgTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAAt                                                                                                                    >  2:517977/1‑139 (MQ=255)           ccgaccgTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAAt                                                                                                                    >  1:1605278/1‑139 (MQ=255)                     tttATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                                          <  1:1785329/139‑1 (MQ=255)                     tttATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                                          <  1:1627437/139‑1 (MQ=255)                             cggcggAGCACAATGCGCTAAAACAGCTGGTGACCAACCTTGCGGCGGCAATGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                                            >  1:359038/1‑113 (MQ=255)                             cggcggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                                            <  2:359038/113‑1 (MQ=255)                             cggcggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATCCCAg                                                                                                               <  1:150998/126‑1 (MQ=255)                             cggcggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATCCCAg                                                                                                               >  2:150998/1‑126 (MQ=255)                              ggcggAGCACAATGCGCTAAAACAGCTGGTGACCAACCTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCAGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAAc                                                                                                 >  1:527366/1‑139 (MQ=255)                               gcggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                                            >  2:1520010/1‑111 (MQ=255)                               gcggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                                            <  1:1520010/111‑1 (MQ=255)                               gcggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                            >  1:403446/1‑127 (MQ=255)                               gcggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                            <  2:403446/127‑1 (MQ=255)                                 ggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCAt                                                                                              >  1:567567/1‑139 (MQ=255)                                   aGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCt                                                                                            >  2:1336371/1‑139 (MQ=255)                                   aGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCt                                                                                            >  1:1425093/1‑139 (MQ=255)                                         aTGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                                      <  2:1605278/139‑1 (MQ=255)                                         aTGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                                      <  2:567567/139‑1 (MQ=255)                                             gcTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                            <  2:1153840/113‑1 (MQ=255)                                             gcTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                            >  1:1153840/1‑113 (MQ=255)                                                aaaaCAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTa                                                                               <  2:1425093/139‑1 (MQ=255)                                                      gCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGgcg                                                                         >  2:213549/1‑139 (MQ=255)                                                         ggTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGa                                                                                                                                    >  2:298948/1‑77 (MQ=255)                                                           tGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCa                                                                    <  2:527366/139‑1 (MQ=255)                                                               cAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAg                                                                                                               >  1:291968/1‑92 (MQ=255)                                                               cAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAg                                                                                                               <  2:291968/92‑1 (MQ=255)                                                                       cggcggCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCATCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCa                                                        >  2:204569/1‑139 (MQ=255)                                                                       cggcggCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTa                                                                                                   >  1:1884268/1‑96 (MQ=255)                                                                       cggcggCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTa                                                                                                   <  2:1884268/96‑1 (MQ=255)                                                                              aaCGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATg                                                                                                                              >  2:495366/1‑62 (MQ=255)                                                                              aaCGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATg                                                                                                                              <  1:495366/62‑1 (MQ=255)                                                                                    ggAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCg                                           >  2:2023243/1‑139 (MQ=255)                                                                                      aaCCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                             <  1:1337608/55‑1 (MQ=255)                                                                                      aaCCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                             >  2:1337608/1‑55 (MQ=255)                                                                                      aaCCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                            >  1:190186/1‑72 (MQ=255)                                                                                      aaCCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                            <  2:190186/72‑1 (MQ=255)                                                                                         ccAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCt                                      <  2:367871/139‑1 (MQ=255)                                                                                          cAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaa                                                                                                           >  2:1780766/1‑69 (MQ=255)                                                                                          cAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaa                                                                                                           <  1:1780766/69‑1 (MQ=255)                                                                                           aGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCAt                                                                                              >  2:1990545/1‑81 (MQ=255)                                                                                           aGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCAt                                                                                              <  1:1990545/81‑1 (MQ=255)                                                                                           aGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCAt                                                                                              <  1:1432284/81‑1 (MQ=255)                                                                                           aGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCAt                                                                                              >  2:1432284/1‑81 (MQ=255)                                                                                            gCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGa                                   >  2:131559/1‑139 (MQ=255)                                                                                            gCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGa                                   >  1:1015940/1‑139 (MQ=255)                                                                                             cAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGaa                                                                                      >  1:520448/1‑86 (MQ=255)                                                                                             cAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                                      <  2:520448/87‑1 (MQ=255)                                                                                               aaCAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                                      <  1:1573908/85‑1 (MQ=255)                                                                                               aaCAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                                      >  2:1573908/1‑85 (MQ=255)                                                                                                aCAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCa                               >  2:1550331/1‑139 (MQ=255)                                                                                                  aGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACt                             <  1:1515837/139‑1 (MQ=255)                                                                                                  aGATCTGACAATCGATGCTGGATCTTGCGCGGGCGAAAGATGGCGAGTTCATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACt                             <  1:706207/139‑1 (MQ=255)                                                                                                  aGATCTGACAATCGATGCTGGAACGTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACt                             >  2:1256577/1‑139 (MQ=255)                                                                                                   gATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTa                            >  1:159512/1‑139 (MQ=255)                                                                                                       tGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGa                        >  1:752265/1‑139 (MQ=255)                                                                                                        ggCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGAt                       <  1:1256577/139‑1 (MQ=255)                                                                                                            aTCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                            >  1:1609008/1‑50 (MQ=255)                                                                                                            aTCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                            <  2:1609008/50‑1 (MQ=255)                                                                                                            aTCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGtt                                                                                                     <  2:1836107/57‑1 (MQ=255)                                                                                                            aTCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGtt                                                                                                     >  1:1836107/1‑57 (MQ=255)                                                                                                              cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                            <  2:192365/48‑1 (MQ=255)                                                                                                              cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                            >  1:192365/1‑48 (MQ=255)                                                                                                              cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCt                 >  2:190389/1‑139 (MQ=255)                                                                                                              cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCGGGCTACAGGCGCGTCAGACGCTCATCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCt                 >  1:912928/1‑139 (MQ=255)                                                                                                                   cTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAAc                                                                                                 >  2:380495/1‑54 (MQ=255)                                                                                                                   cTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAAc                                                                                                 <  1:380495/54‑1 (MQ=255)                                                                                                                      gAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAaacaac         <  1:131559/139‑1 (MQ=255)                                                                                                                           ttCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGCCGCTGGAATCACTACAGATGACGCTAAAACAACCTTTa    >  1:1628398/1‑139 (MQ=255)                                                                                                                             ccGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCTTTAGa  <  1:190389/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                     CGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGTGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCTTTAGA  >  NC_000913/2438418‑2438680

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