Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,892,710	A→G	Y176H (TAT→CAT)	yoaA ←	putative 5' to 3' DNA helicase implicated in DNA repair

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,892,710	0	A	G	100.0%	219.6 / NA	63	Y176H (TAT→CAT)	yoaA	putative 5' to 3' DNA helicase implicated in DNA repair
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/35);  total (28/35)

TCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCGCAG  >  NC_000913/1892577‑1892839                                                                                                                                      |                                                                                                                                  tCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCgg                                                                                                                             >  2:615414/1‑140 (MQ=255)  ccGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCtgg                                                                                                                            >  2:649447/1‑140 (MQ=255)   cGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGc                                                                                                                           >  2:737519/1‑140 (MQ=255)       tCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCTGGCAGTc                                                                                                                       <  1:649447/140‑1 (MQ=255)       tCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTc                                                                                                                       <  1:75259/140‑1 (MQ=255)       tCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTc                                                                                                                       <  2:939644/140‑1 (MQ=255)         cGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGc                                                                                                                     <  1:737519/140‑1 (MQ=255)            gATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGc                                                                                                                  <  1:793317/140‑1 (MQ=255)                 gTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACAACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGa                                                                                                             <  1:1018444/140‑1 (MQ=255)                  ttCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCggg                                                                                                                            <  2:213616/124‑1 (MQ=255)                  ttCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCggg                                                                                                                            >  1:213616/1‑124 (MQ=255)                  ttCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGAc                                                                                                            >  1:939562/1‑140 (MQ=255)                      ccAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAGAGCAATCATTATGCATCGGGCAGTCGCTGCCAAGTCAGtt                                                                                                        >  2:899871/1‑140 (MQ=255)                         ataccACTCACTTTAACCCCCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCGCATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTc                                                                                                     <  2:442022/137‑1 (MQ=255)                           tCCACTCTCTTTAGCCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGt                                                                                                   <  2:402285/140‑1 (MQ=255)                              aCTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTgg                                                                                                <  2:602669/140‑1 (MQ=255)                                            ccaTATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTCCGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGc                                                                                  >  2:908877/1‑140 (MQ=255)                                            ccaTATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGc                                                                                  >  2:438232/1‑140 (MQ=255)                                              atatCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCgg                                                                                >  2:320113/1‑140 (MQ=255)                                              atatCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATAGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCgg                                                                                >  2:837856/1‑140 (MQ=255)                                               tatCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGc                                                                               <  2:812095/140‑1 (MQ=255)                                                 tCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCa                                                                             <  1:110181/140‑1 (MQ=255)                                                 tCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCa                                                                             <  1:615414/140‑1 (MQ=255)                                                       aGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCt                                                                        >  2:370002/1‑139 (MQ=255)                                                       aGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCt                                                                        <  1:370002/139‑1 (MQ=255)                                                                  atgGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCg                                                            >  2:428254/1‑140 (MQ=255)                                                                         accaccaccACATCGGCGTCCATCGCTTTTTTACGTGCTTTGCCCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCt                                                     <  2:939562/140‑1 (MQ=255)                                                                         accaccaccACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCt                                                     >  1:203401/1‑140 (MQ=255)                                                                         accaccaccACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCt                                                     >  1:100748/1‑140 (MQ=255)                                                                             ccaccacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACg                                                 <  1:428254/140‑1 (MQ=255)                                                                              caccacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGc                                                >  1:189465/1‑140 (MQ=255)                                                                              caccacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGc                                                <  1:908877/140‑1 (MQ=255)                                                                               accacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCa                                               >  2:583855/1‑140 (MQ=255)                                                                               accacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCa                                               >  2:498136/1‑140 (MQ=255)                                                                                ccacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAg                                              <  2:100748/140‑1 (MQ=255)                                                                                ccacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAg                                              <  1:320113/140‑1 (MQ=255)                                                                                ccacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAg                                              >  2:520623/1‑140 (MQ=255)                                                                                ccacATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAg                                              >  1:691473/1‑140 (MQ=255)                                                                                  acaTCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGt                                            <  2:806327/140‑1 (MQ=255)                                                                                    aTCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGc                                          >  2:518335/1‑140 (MQ=255)                                                                                      cGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTg                                        >  1:134652/1‑140 (MQ=255)                                                                                        gCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGt                                            <  2:87122/134‑1 (MQ=255)                                                                                        gCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGt                                            >  1:87122/1‑134 (MQ=255)                                                                                         cGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCTAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGata                                     >  2:822794/1‑140 (MQ=255)                                                                                         cGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGata                                     >  1:87325/1‑140 (MQ=255)                                                                                          gTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGatat                                    <  2:134652/140‑1 (MQ=255)                                                                                            ccATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCTAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCa                                  <  1:822794/140‑1 (MQ=255)                                                                                              aTCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGAGCTGATATCAcc                                >  1:412649/1‑140 (MQ=255)                                                                                                  cTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGTCAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCg                            <  2:147330/140‑1 (MQ=255)                                                                                                       ttACGTGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGt                       <  1:497803/140‑1 (MQ=255)                                                                                                            tGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGAt                  <  2:203401/140‑1 (MQ=255)                                                                                                            tGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGAt                  <  1:899871/140‑1 (MQ=255)                                                                                                            tGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGAt                  <  1:498136/140‑1 (MQ=255)                                                                                                            tGCTTTGACCACAAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGAt                  <  2:87325/140‑1 (MQ=255)                                                                                                                       caAAGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAg       <  1:583855/140‑1 (MQ=255)                                                                                                                        aaaGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGAt                  <  1:725674/128‑1 (MQ=255)                                                                                                                        aaaGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGAt                  <  1:522179/128‑1 (MQ=255)                                                                                                                        aaaGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGAt                  >  2:522179/1‑128 (MQ=255)                                                                                                                        aaaGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGAt                  >  2:725674/1‑128 (MQ=255)                                                                                                                         aaGCAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAgcg     >  2:970250/1‑140 (MQ=255)                                                                                                                            cAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCgcag  <  1:518335/140‑1 (MQ=255)                                                                                                                            cAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCgcag  <  1:837856/140‑1 (MQ=255)                                                                                                                            cAATCTTTATGCATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCgcag  <  2:691473/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                  TCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTTTAACCACCATATCCGCCAGAAAGAGATGATGGTTTACCACCACCACATCGGCGTCCATCGCTTTTTTACGTGCTTTGACCACAAAGCAATCTTTATACATCGGGCAGTCGCTGCCAAGACAGTTGTCGTTGGTGCTGGTGACCAGCGGCCACGCCTGTGAATCTTCCGCCACGCTGACGCAGGTGCTGATATCACCATCGACTGTTTGATTAGACCAGGAGCGCAG  >  NC_000913/1892577‑1892839

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