Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,744,030	T→C	Y136H (TAT→CAT)	yfiB →	lipoprotein YfiB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,744,030	0	T	C	100.0%	121.9 / NA	37	Y136H (TAT→CAT)	yfiB	lipoprotein YfiB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/18);  total (19/18)

ACAAATTGCTCCCGGAAAGCCAGCAACAGATCCAAACCATGGCAGCTAAATTGGCCTCGACAGGGCTAACACATGCCCGTATGGATGGACACACCGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGAC  >  NC_000913/2743806‑2744268                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               aCAAATTGCTCCCGGAAAGCCAGCAACAGATCCAAACCATGGCAGCTAAATTGGCCTCGACAGGGCTAACACATGCCCGTATGGATGGACACACCGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGcc                                                                                                                                                                                                                      <  1:145180/251‑1 (MQ=255)    aaTTGCTCCCGGAAAGCCAGCAACAGATCCAAACCATGGCAGCTAAATTGGCCTCGACAGGGCTAACACATGCCCGTATGGATGGACACACCGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGccc                                                                                                                                                                                                                     >  1:281696/1‑249 (MQ=255)                  aGCCAGCAACAGATCCAAACCATGGCAGCTAAATTGGCCTCGACAGGGCTAACACATGCCCGTATGGATGGACACACCGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCgaga                                                                                                                                                                                                      >  1:290240/1‑250 (MQ=255)                    ccAGCAACAGAGCCAAACCATGGCAGCTAAATTGGCCTCGACAGGGCTAACACATGCCCGTATGGATGGACACACCGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCAAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAAcc                                                                                                                                                                                                   <  2:249508/251‑1 (MQ=255)                       gCAACAGATCCAAACCATGGCAGCTAAATTGGCCTCGACAGGGCTAACACATGCCCGTATGGATGGACACACCGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAAccgcc                                                                                                                                                                                                <  1:361020/251‑1 (MQ=255)                                     ccATGGCAGCTAAATTGGCCTCGACAGGGCTAACACATGCCCGAATGGATGGACACACCGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAgtggtg                                                                                                                                                                                  >  1:227147/1‑251 (MQ=255)                                                    tGGCCTCGACAGGGCTAACACATGCCCGTATGGATGGACACACCGATAACTATGGGGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACACATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGAGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTa                                                                                                                                                                    >  1:50123/1‑250 (MQ=255)                                                     ggCCTCGACAGGGCTAACACATGCCCGTATGGATGGACACACCGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTaa                                                                                                                                                                   >  1:264257/1‑250 (MQ=255)                                                      gCCTCGACAGGGCTAACACATGCCCGTATGGATGGACACACCGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGAGTACTACCCCTTaaaa                                                                                                                                                                 <  1:286272/251‑1 (MQ=255)                                                       ccTCGACAGGGCTAACCCATGCCCGTATGGATGGACACACCGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGTTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAac                                                                                                                                                                <  1:133876/251‑1 (MQ=255)                                                                  cTAACACATGCCCGTATGGATGGACACACCGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAAc                                                                                                                                                      >  1:124995/1‑250 (MQ=255)                                                                    aaCACATGCCCGTATGGATGGACACACCGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAAc                                                                                                                                                      >  1:125548/1‑248 (MQ=255)                                                                        caTGCCCGTATGGATGGACACACCGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATaa                                                                                                                                                >  2:151892/1‑250 (MQ=255)                                                                        aaggCCCGTATGGATGGACACACCGATAACTATGGTGAGGACAGTTACAACGAAGGATTATCATTGAAACGGGCGAATGTCGTGGCCGAGGCATGG‑CTATGGATGGGCAAATTCCACGCGGCAATCTCACCACACAGGGTTGAGGAAAAAAACAGCCCATAGCCAGTAACAAGGCCGCTCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATaaaa                                                                                                                                              <  1:179901/248‑1 (MQ=255)                                                                                      tggaCACACCGTTAAATATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTACACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCTCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGGTAAAAAAGGgccagcca                                                                                                                                  >  2:321030/1‑250 (MQ=255)                                                                                         acacacCGATAACTATGGTGAAGTCAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTg                                                                                                                              <  1:363078/251‑1 (MQ=255)                                                                                              ccGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGAACAAATTCCACGCAGCAATCTCCCCCCACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAg                                                                                                                         <  1:81013/251‑1 (MQ=255)                                                                                                                       cAACGAAGGATTATCATTGAAACGGGCGTATTTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCAAGTAACAAGACCGCCCAGGGCCGCGCCAAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGTCAGCCCTTCTTAACAGGATGTCGCTTaa                                                                                                <  1:188966/251‑1 (MQ=255)                                                                                                                        aaCGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAg                                                                                               <  1:40583/251‑1 (MQ=255)                                                                                                                          cGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCg                                                                                             >  2:245727/1‑251 (MQ=255)                                                                                                                                 ttATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCCTAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCAAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCtt                                                                                      >  2:47691/1‑251 (MQ=255)                                                                                                                                  tATCATTGAAACGGGCGAATTTCGTGGCCGATGCATGGGCTATGGGGGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATc                                                                                        <  1:371777/248‑1 (MQ=255)                                                                                                                                      attaaaaCGGGCGAATTTCGTGGCGGATGAATGGTCTTTGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACACCCCATAGCCAGTAACAAGTCCGCCCAGGGCCGTGGTGAGAACCGCCGCGTCGCAGTTGTTATTAATTCCCCTTAAAACACCAGCAAACGGATAAAAAAGGGCCATCCGATTGGCCAGCCCTTCGTAACAGGATGTCGCTTAAGCGAAATCTTAGTTa                                                                                 <  1:169448/247‑1 (MQ=255)                                                                                                                                           aaCGGGCGAATGTCGTGGCCGATGCCTGGGCTATGGGTGGCCAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGACAAAAACATCCCCTAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCACACAGATCAACAAGGGCCAGCCCATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAgccg                                                                            >  2:371948/1‑248 (MQ=255)                                                                                                                                               ggCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGctct                                                                        >  2:135254/1‑251 (MQ=255)                                                                                                                                                 cGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCttt                                                                      >  2:347655/1‑251 (MQ=255)                                                                                                                                                  gAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGACCCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTg                                                                     >  1:242119/1‑251 (MQ=255)                                                                                                                                                  gAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATCCCACGCAGCAATCTCCCCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTg                                                                     <  2:124995/251‑1 (MQ=255)                                                                                                                                                      gTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATa                                                                  <  1:115279/250‑1 (MQ=255)                                                                                                                                                            gCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGagc                                                             >  2:279565/1‑249 (MQ=255)                                                                                                                                                                       ggCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGt                                                 >  1:82348/1‑250 (MQ=255)                                                                                                                                                                          tATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGc                                             >  1:347603/1‑251 (MQ=255)                                                                                                                                                                                    aaaaTACCACGCAACAACCTCACCACACAGGGTTGAGGAAAAAAACATCCCATAGCCAGTAACAAGCCCGCCCGGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGgt                                   <  2:173131/250‑1 (MQ=255)                                                                                                                                                                                                       tcACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAg                                                                                                                                                    >  1:94782/1‑119 (MQ=255)                                                                                                                                                                                                       tcACCACACAGGGTTTAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAg                                                                                                                                                    <  2:94782/119‑1 (MQ=255)                                                                                                                                                                                                               caGGGTTTAGAAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGcac        <  1:10008/251‑1 (MQ=255)                                                                                                                                                                                                                      tAGGAAAAAAACATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCacgac  <  1:250410/250‑1 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               ACAAATTGCTCCCGGAAAGCCAGCAACAGATCCAAACCATGGCAGCTAAATTGGCCTCGACAGGGCTAACACATGCCCGTATGGATGGACACACCGATAACTATGGTGAAGACAGTTACAACGAAGGCTTATCATTGAAACGGGCGAATGTCGTGGCCGATGCATGGGCTATGGGTGGACAAATTCCACGCAGCAATCTCACCACACAGGGTTTAGGAAAAAAATATCCCATAGCCAGTAACAAGACCGCCCAGGGCCGCGCCGAGAACCGCCGCGTCGCAGTGGTGATTACTACCCCTTAAAACACCAGCAAACAGATAAAAAAGGGCCAGCCAATTGGCCAGCCCTTCTTAACAGGATGTCGCTTAAGCGAAATCTTAGTTAAGACGCTCTTTGATACGAGCAGCCTTACCAGTACGCTCACGCAGGTAGTACAGTTTAGCTTTACGAACAGCACCACGAC  >  NC_000913/2743806‑2744268

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