Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,133,273	G→A	G234S (GGC→AGC)	flgE →	flagellar hook protein FlgE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,133,273	0	G	A	100.0%	106.9 / NA	37	G234S (GGC→AGC)	flgE	flagellar hook protein FlgE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/18);  total (19/18)

GCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATT  >  NC_000913/1133027‑1133515                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   gCATCGATGCAGGTCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGAGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACTGCATTGCGAAGACAGCGACACCACTGGAATTTAATGCTAATAGcac                                                                                                                                                                                                                                                >  1:616044/1‑251 (MQ=255)    tCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACAtt                                                                                                                                                                                                                                             >  1:97748/1‑251 (MQ=255)     cGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACAtt                                                                                                                                                                                                                                             <  2:42152/250‑1 (MQ=255)     cGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACAtt                                                                                                                                                                                                                                             >  1:266282/1‑250 (MQ=255)       aTGCAGATCAACCTGAATTCCCGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGt                                                                                                                                                                                                                                          >  1:145606/1‑251 (MQ=255)                 aCCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACACAGGATAGCAGTGATCCAAACAGCATTGCGAAGACCGCGACAACACTGGAGTTTAATGCTAATAGCACATTAGTGGATGGTGCg                                                                                                                                                                                                                                >  1:56404/1‑251 (MQ=255)                   cTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGACAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACCACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGAt                                                                                                                                                                                                                              >  2:328235/1‑251 (MQ=255)                     gAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATgg                                                                                                                                                                                                                            >  1:12896/1‑251 (MQ=255)                        ttCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCg                                                                                                                                                                                                                          >  1:332146/1‑250 (MQ=255)                             gTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGaataat                                                                                                                                                                                                                    <  1:600065/251‑1 (MQ=255)                                       tCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAAc                                                                                                                                                                                                            <  2:602850/249‑1 (MQ=255)                                        ccTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTGTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAAGTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCg                                                                                                                                                                                                          <  2:425114/250‑1 (MQ=255)                                                               cAGCAATGCGGATAGCTATAACAAAAAATGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGGCTACGTTGTGAAGACCGGGGATAATAACTGTCAGGGCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCTAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACAGTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCg                                                                                                                                                                                  <  2:591882/251‑1 (MQ=255)                                                                  cAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGcc                                                                                                                                                                                >  2:219374/1‑250 (MQ=255)                                                                                 taACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTAATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCt                                                                                                                                                                >  2:494935/1‑251 (MQ=255)                                                                                    cAAAAAAGGTTGGGTGACTGTGTTCGAAAGTCAGGGTAAGGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGAAATTTAATGCTAATAGCACTTTAGTGGATGGTGCGATGGCGAATAATATCGCACCCGGCGCAATTACCGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTcc                                                                                                                                                             <  2:145606/251‑1 (MQ=255)                                                                                    aaaaaaaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATACCTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATGGCGAAGACAGCGCCAACACTGGAATTAAATGCTAATAGCACATTAGGGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTACCGTTGCAGAACCCGCCACGTTTAGTCTGAGCTTcc                                                                                                                                                             <  1:543247/250‑1 (MQ=255)                                                                                      aaaaaGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTc                                                                                                                                                           >  2:500716/1‑251 (MQ=255)                                                                                                     cTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTa                                                                                                                                                                                                 >  1:212521/1‑198 (MQ=255)                                                                                                     cTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTa                                                                                                                                                                                                 <  2:212521/198‑1 (MQ=255)                                                                                                            cGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGGGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGAGGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTACCGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCg                                                                                                                                     <  1:72682/251‑1 (MQ=255)                                                                                                                           tGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATtg                                                                                                                      >  1:85538/1‑251 (MQ=255)                                                                                                                               catgacatgaGCGTCTATTTTGTGAAGACCGGGGATAATAACGGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACGGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGAGGGCGAATAATATCGAAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTgg                                                                                                                   <  1:241310/250‑1 (MQ=255)                                                                                                                                atgacatgaGCGTCTAATTTGTGAAGACCGGGGATAATAACTGGCAGGACTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGAGACAACACTGGAATTTAATGCTAATAGCACATGAGTGGATGGTGCGATGGCGAATAATATAGCAACCGGCGCAATTAACGGTGCAGAAACCGCCACGTTTAGTCTGAGCTTACTCAACTCCAAGCAGAAAAATACCGGCGCTAACAATATTGTGGc                                                                                                                  >  2:465848/1‑250 (MQ=255)                                                                                                                                     atgaGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAaccac                                                                                                            <  1:500716/251‑1 (MQ=255)                                                                                                                                                    gggAAGACCGGGAATAATAACTGGCAGGTCACCACCAAGGATTGCAGGAATCCAAACAGCATTGGGAAGACAGCCACAACACTGGAATTTAATGCTAATAGCCCATTAGTGGATGGGGCGATGGCGAATAATATCGCATCCGGCGCAATTAACGGTGCAGAACCCGCCTCGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACa                                                                                              <  2:494424/248‑1 (MQ=255)                                                                                                                                                         aacgggggATAATAACTGGCAGGTCTACCACTAGGATAGCAGTGATCCAAACAGCATCGCGAAGACAGCGCCAACACTGGAATTGAATCTTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCACCCGGCGCAATTAACGGGGCAGAGCCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCTAGAACGGCTACAAACCgg                                                                                        <  2:286716/247‑1 (MQ=255)                                                                                                                                                                gATAATAACTGGCAGGTCTACCCCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATAGCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCt                                                                                 <  1:602264/251‑1 (MQ=255)                                                                                                                                                                 ataatctcTGGCAGGTCTACACCCAGGATCGCCGTGATCCAAATAGCATTGCGTTGTCAGCGCCAACACTGGAATTTAATGCTAATAGCTCATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGCAGCAATTTACGGTGCAGTACCCGCCACGTTTAGTCCGAGCTTCCTCAACTCCACTCAACCCAATACCGGCGCTAACAATATTGTGTCCACCACCCAGAAAGGCTACACACCgcacgatatt                                                                                >  2:436636/8‑242 (MQ=255)                                                                                                                                                                  taataaCTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTgg                                                                               >  2:493275/1‑251 (MQ=255)                                                                                                                                                                                                   agATCCAAACAGCATTGCGAAGACAGCGACAAAACTGGAATTTAATGCTAATAGCACCTTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAACTACCGGCGATAAAAATATTGTGTCAAGCACCCAGAACGGCTACAAACCGGGCGATCTGGTCATCTATCAAATCACTGATGACGGTACCGTTg                                              >  2:487372/2‑251 (MQ=255)                                                                                                                                                                                                         aaaaagaattggaaaGACAGCGACAACACTGAAATTTAATGCAAATAGCAAATTAGTGGAGGGTGCGATGGAGAATAAAATCGCAACCGGCGCAATTAACGGTGCAAACCCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCa                                        <  2:9642/238‑1 (MQ=255)                                                                                                                                                                                                              gCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGAGGGCGAATAATATCGCAACCGGCGCAATTAACGGTGTAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTAt                                   <  1:328235/251‑1 (MQ=255)                                                                                                                                                                                                               cATTGCGAAGACAGCGACCCCACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTAtt                                  >  1:155161/1‑251 (MQ=255)                                                                                                                                                                                                                              gACAACACTGGAATTTAATGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAAccc                   >  1:547569/1‑251 (MQ=255)                                                                                                                                                                                                                                caacaCTGGAATTTAATGCAAATAGCACATTAGTGGATGGTGCGATGGCGAATAATACCGCAACCGGCGCAATAACCGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCaa                 <  1:277481/251‑1 (MQ=255)                                                                                                                                                                                                                                               aTGCTAATAGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTACCGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGAtt  <  2:155161/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   GCATCGATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTGAAGACCGGGGATAATAACTGGCAGGTCTACACCCAGGATAGCAGTGATCCAAACAGCATTGCGAAGACAGCGACAACACTGGAATTTAATGCTAATGGCACATTAGTGGATGGTGCGATGGCGAATAATATCGCAACCGGCGCAATTAACGGTGCAGAACCCGCCACGTTTAGTCTGAGCTTCCTCAACTCCATGCAGCAAAATACCGGCGCTAACAATATTGTGGCAACCACCCAGAACGGCTACAAACCGGGCGATCTGGTGAGTTATCAAATCAATGATGACGGTACGGTTGTCGGCAACTATTCCAACGAACAAACCCAACTGCTGGGGCAGATT  >  NC_000913/1133027‑1133515

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