Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,133,175	G→A	T13I (ACC→ATC)	yghR ←	putative ATP‑binding protein YghR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,133,175	0	G	A	100.0%	101.2 / NA	35	T13I (ACC→ATC)	yghR	putative ATP‑binding protein YghR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/17);  total (18/17)

TGGTTTTCATGCTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTC‑GCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGG  >  NC_000913/3132933‑3133391                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                          tGGTTTTCATGCTTTGGGTTTTATAGGATCTGGTAGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGAAATCGCTTGATTTTGTCGCGGTCTTCGCCGGGGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGCTGTTGCAGCGAATTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTACGGCGCGATTGGACtt                                                                                                                                                                                                                   >  1:117379/1‑250 (MQ=255)  ggTTTTCATGCTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTAt                                                                                                                                                                                                                  >  1:241187/1‑251 (MQ=255)        cATGCTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCCGACGCCAACCAACGGCAATCGTTTGCTTTTGTCGCCGTCTTCGCCGGAGCGCCGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAAc                                                                                                                                                                                                              >  1:374539/1‑249 (MQ=255)           gCTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGGCGGTGTTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACagttt                                                                                                                                                                                                         <  2:130794/251‑1 (MQ=255)            cTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGACCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACagtttg                                                                                                                                                                                                        >  1:244482/1‑251 (MQ=255)              ttGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACagtttg                                                                                                                                                                                                        >  1:356994/1‑249 (MQ=255)              ttGGGTTTTCGAGGAGTAGGCCGCCAGTCTCCGTTCAAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTTGCGCCGTCATCGCCGGAGAGCAGCCCCAGATTGCGCCGCTCGGGTTGCCAGTGTTGTTGCAGCGATTTCACCAGGGCGATGGTCAGAGTGGATTTACCGGTGCCATCACAGCCGACCACCTCAAGCAGCCGCGGAATGTATTTGGGCTGCGGCGCGATTGCACTATGAACattttga                                                                                                                                                                                                       <  2:374539/250‑1 (MQ=255)                              gtGGCCTCCAGTCGCCGTTCCATCCAGACGCCAACCACCGGCAATCGTGTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGGGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCa                                                                                                                                                                                       <  2:12022/249‑1 (MQ=255)                                           gCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAatcatc                                                                                                                                                                         >  2:7761/1‑251 (MQ=255)                                             cGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGCTCGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATccc                                                                                                                                                                       >  1:32960/1‑251 (MQ=255)                                                       agACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCCttt                                                                                                                                                             >  1:271480/1‑251 (MQ=255)                                                                      gcatgctttggatgttttcgcccgcTTCCCCGGTGAGCAGCCCCAGAGAGCGCCGCGCGGTTTGCCATTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGTTGCTATCCCACCCGACCACCGCAATCAGCCCTGGAATGTAGTTGGGCTGCTTCGCGATTGTACTATTAACAGTTTGAGGTTGTATGGCATCCATTCGTTAAAGCAGCCCCTATCCATTTTGGTTAACAGGAAAg                                                                                                                                              <  2:117379/227‑1 (MQ=255)                                                                      gCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGCGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGaa                                                                                                                                                >  1:369691/1‑249 (MQ=255)                                                                               tGATTTTGTCGCAGTCCTCGCCGGCGAGAGGCCAAAGATAGCGCCGCTCGGGTTGCCCGTGTTGTTGCAGCGATTTAAACAGGTCCGTGGTAAGTGTGGATTGACAGGTGCCATCACACCCGACAACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCTCGATAGTACGATGAACAGTCTGCGTTTGACGAGCATCCATTCGTTAAATCATCCACTATCCATTATGGTGATAAGGAACGCGCGGCGt                                                                                                                                      >  2:383813/1‑250 (MQ=255)                                                                                       tCGCCGTCTTCGCCGGAGCGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCCCCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGCGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTgcgc                                                                                                                             >  1:69054/1‑251 (MQ=255)                                                                                         gCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTTGATTTACCGGTGCCATCCCACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGt                                                                                                                           <  2:271480/251‑1 (MQ=255)                                                                                                     ggAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAAAAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGc                                                                                                               >  1:371452/1‑251 (MQ=255)                                                                                                            agcCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCTTAGTGCCGCTTGCAGAAttt                                                                                                        >  2:112658/1‑251 (MQ=255)                                                                                                                               cGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTC‑Gcc                                                                                     <  1:4762/250‑1 (MQ=255)                                                                                                                                 gTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTC‑Gccccc                                                                                  >  1:109806/1‑249 (MQ=255)                                                                                                                                   ttGCCAGTGTTGTTGCAGCGATTTACCCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCGCCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCAGTCGGTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGTGGCGTGAATTGCGCGTAGTGACGCTTGCAGAAGTTCATCGGCGGGGTGTC‑GCCCATc                                                                                 <  2:375250/250‑1 (MQ=255)                                                                                                                                           gttgttGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCCCACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTTGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTC‑GCCTATCCAGTTCaa                                                                         <  2:399833/250‑1 (MQ=255)                                                                                                                                                  cAGCGATGTCACCAGGGCGGTGGTCAGTGTGGATTTACCGGTGCCATCACCCCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGAGTGTC‑GCCCATCCAGTTCAAGGATCttt                                                                 <  2:224207/251‑1 (MQ=255)                                                                                                                                                   aGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGGATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTC‑GCCCATCCAGTTCAAGGATCTTTg                                                                <  1:46774/251‑1 (MQ=255)                                                                                                                                                    gCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTC‑GCCCATCCAGTTCAAGGATCTTTgc                                                               <  1:294861/251‑1 (MQ=255)                                                                                                                                                               aGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTC‑GCCCATCCCGTTCAAGGAGCTTTGCGCCATTGAATg                                                    >  2:140261/1‑251 (MQ=255)                                                                                                                                                                    ggtggtCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTC‑GCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATg                                                    <  2:301527/246‑1 (MQ=255)                                                                                                                                                                    ggtggtCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTG‑GCCCCTCCCGTTCAAGGATCTTTGCGCCATTGAATg                                                    >  1:301527/1‑246 (MQ=255)                                                                                                                                                                          cAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTC‑GCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgcggc                                         <  1:89228/251‑1 (MQ=255)                                                                                                                                                                           aGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTC‑GCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgcggcg                                        <  1:187990/251‑1 (MQ=255)                                                                                                                                                                                  aTTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTC‑GCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGc                                 <  1:196461/251‑1 (MQ=255)                                                                                                                                                                                  aTTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGAAGTTGGGCTGCGGCGCGATTGTCCTATTAACAGTTTGAGTTTGTAGTGAATCAATTCGTTAAATCATCCACTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTC‑GCCCACCCAGCTCAAGGATCTTTTGGCCATTGCATATCAgatggggggcaacggc                                 >  1:257034/1‑236 (MQ=255)                                                                                                                                                                                      aCCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTC‑GCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGAtt                             <  1:81593/251‑1 (MQ=255)                                                                                                                                                                                                             ccGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTC‑GCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtg      <  1:40240/251‑1 (MQ=255)                                                                                                                                                                                                                 aaTCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAAAAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCTGCCC‑TCCAACTCAAGGATCTTTGCGCCCTTTCATCTTAGTTGCGGCGTAACGGCGGTTTTTTCCTGTAACGcatccccgtggtgg  >  1:199852/1‑238 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                          TGGTTTTCATGCTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTC‑GCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGG  >  NC_000913/3132933‑3133391

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