Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,795,485	G→A	A23V (GCC→GTC)	ihfA ←	integration host factor (IHF), DNA‑binding protein, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,795,485	0	G	A	100.0%	224.4 / NA	57	A23V (GCC→GTC)	ihfA	integration host factor (IHF), DNA‑binding protein, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/36);  total (21/36)

GAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCT  >  NC_000913/1795339‑1795632                                                                                                                                                   |                                                                                                                                                    gAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACAt                                                                                                                                                  <  1:76667/150‑1 (MQ=255)  aaTATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCCCGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATc                                                                                                                                                 <  2:110885/150‑1 (MQ=255)     atCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCg                                                                                                                                              >  1:964550/1‑150 (MQ=255)      tCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCg                                                                                                                                              >  1:34727/1‑149 (MQ=255)       ccTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCtt                                                                                                                                           >  1:721098/1‑151 (MQ=255)              cgctttCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAg                                                                                                                                     >  1:67097/4‑150 (MQ=255)                ttttCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGcc                                                                                                                                   <  2:712715/150‑1 (MQ=255)                ttttCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGcc                                                                                                                                   <  1:341510/150‑1 (MQ=255)                ttttCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGcc                                                                                                                                   <  2:460237/150‑1 (MQ=255)                    cGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGc                                                                                                                              <  2:162124/151‑1 (MQ=255)                     gggTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAg                                                                                                                               >  2:753939/1‑149 (MQ=255)                      ggTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCt                                                                                                                             <  1:555231/150‑1 (MQ=255)                      ggTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTCGCTAAGCCCAAGCt                                                                                                                             <  1:715061/150‑1 (MQ=255)                           cGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCaa                                                                                                                       >  2:538751/1‑151 (MQ=255)                             tCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAAc                                                                                                                     <  1:487693/151‑1 (MQ=255)                             tCCCGGGCGTTGATTCTTAGCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAAc                                                                                                                     <  1:576949/151‑1 (MQ=255)                               ccGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAg                                                                                                                   <  1:11021/151‑1 (MQ=255)                                   gCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGata                                                                                                                <  2:386717/150‑1 (MQ=255)                                    cGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATAtt                                                                                                              <  1:757327/151‑1 (MQ=255)                                        gATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGa                                                                                                          >  1:853001/1‑151 (MQ=255)                                        gATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGa                                                                                                          >  2:5159/1‑151 (MQ=255)                                        gATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGa                                                                                                          >  1:550319/1‑151 (MQ=255)                                         aTTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGAc                                                                                                         <  2:27370/151‑1 (MQ=255)                                          ttCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACa                                                                                                        >  1:199205/1‑151 (MQ=255)                                                   cGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGGTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCtt                                                                                               <  2:1053426/151‑1 (MQ=255)                                                                tACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCAt                                                                                  <  2:646501/151‑1 (MQ=255)                                                                 aCCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATa                                                                                 >  1:417614/1‑151 (MQ=255)                                                                   cAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAgg                                                                               <  2:152396/151‑1 (MQ=255)                                                                     aaaCCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGtt                                                                             >  2:896793/1‑151 (MQ=255)                                                                     aaaCCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGtt                                                                             <  2:306180/151‑1 (MQ=255)                                                                     aaaCCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGtt                                                                             >  1:855219/1‑151 (MQ=255)                                                                     aaaCCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGtt                                                                             >  1:695795/1‑151 (MQ=255)                                                                      aaCCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTc                                                                            <  1:554129/151‑1 (MQ=255)                                                                            agagTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATcc                                                                       <  1:778008/150‑1 (MQ=255)                                                                            agagTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTCTTGTAAGCGCCATAGGTTCAATc                                                                        <  1:843640/149‑1 (MQ=255)                                                                             gagTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATcc                                                                       >  1:598449/1‑149 (MQ=255)                                                                                   cACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAatga                                                               <  2:1040949/151‑1 (MQ=255)                                                                                   cACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAatga                                                               >  1:139074/1‑151 (MQ=255)                                                                                        gTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGc                                                            >  1:651237/1‑149 (MQ=255)                                                                                             ccGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAAt                                                     <  1:988990/151‑1 (MQ=255)                                                                                               gTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCg                                                   >  1:509622/1‑151 (MQ=255)                                                                                                  ttCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGctc                                                <  2:392807/151‑1 (MQ=255)                                                                                                        accgccACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGAAATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTaa                                          <  2:462907/149‑1 (MQ=255)                                                                                                                 gATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTa                                  >  2:191337/1‑150 (MQ=255)                                                                                                                   tctcTTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTacac                               <  1:753678/151‑1 (MQ=255)                                                                                                                    ctctTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTAcaca                              <  2:431589/151‑1 (MQ=255)                                                                                                                    ctctTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTAcaca                              <  2:964893/151‑1 (MQ=255)                                                                                                                     tctTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACAt                             <  1:538517/151‑1 (MQ=255)                                                                                                                      ctTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACAtt                            <  2:853292/151‑1 (MQ=255)                                                                                                                      ctTCGAAAAACAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACAtt                            <  2:1035299/151‑1 (MQ=255)                                                                                                                                cAGTTCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGt                  >  2:586073/1‑151 (MQ=255)                                                                                                                                   ttCAACCAGTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGt                  <  1:5158/148‑1 (MQ=255)                                                                                                                                           gTTCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAAt         >  1:223487/1‑149 (MQ=255)                                                                                                                                            ttCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGACAAt         <  2:384797/148‑1 (MQ=255)                                                                                                                                             tCTTTGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATctcc     <  1:798145/151‑1 (MQ=255)                                                                                                                                               tttGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATctcctc   <  1:110640/151‑1 (MQ=255)                                                                                                                                                ttGACATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCctct  <  2:819087/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GAATATCCTCGCCCGTTTTCGGGTTACGTCCCGGGCGTTGATTCTTATCACGCAGATCGAAGTTACCAAAACCAGAGAGTTTCACCTGTTCGCCGTTTTCCAGAGCGCGACGGATCTCTTCGAAAAACAGTTCAACCAGTTCTTTGGCATCCCGCTTGCTAAGCCCAAGCTTATCAAACAGATATTCTGACATTTCAGCTTTTGTAAGCGCCATAGGTTCAATCCCTCAATGATGCCTGGAATCGCTCTTTTAATGCCTCTACACATTTGGCGACGGTAGCGGCAATCTCCTCT  >  NC_000913/1795339‑1795632

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