Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,815,591	A→G	P575P (CCA→CCG)	katE →	catalase HPII, heme d‑containing

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,815,591	0	A	G	100.0%	208.5 / NA	58	P575P (CCA→CCG)	katE	catalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/25);  total (33/25)

GTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTG  >  NC_000913/1815443‑1815735                                                                                                                                                     |                                                                                                                                                 gTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGcc                                                                                                                                                >  1:210661/1‑151 (MQ=255)  tttCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCt                                                                                                                                               >  1:199014/1‑151 (MQ=255)     cAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCg                                                                                                                                            >  2:537429/1‑151 (MQ=255)      aGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCgg                                                                                                                                           >  2:936266/1‑151 (MQ=255)        ttttGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGAc                                                                                                                                         >  1:423951/1‑151 (MQ=255)         tttGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGAc                                                                                                                                         <  1:527342/150‑1 (MQ=255)         tttGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACg                                                                                                                                        >  2:925902/1‑151 (MQ=255)         tttGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACg                                                                                                                                        >  1:722540/1‑151 (MQ=255)          ttGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGAc                                                                                                                                         <  1:459128/149‑1 (MQ=255)          ttGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGt                                                                                                                                       >  1:817832/1‑151 (MQ=255)           tGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTc                                                                                                                                      >  1:811604/1‑151 (MQ=255)                tAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGCTCTCACTCTGGCCCAGGCGGGGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACgg                                                                                                                                 >  2:764699/1‑151 (MQ=255)                   gCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCt                                                                                                                              <  1:319821/151‑1 (MQ=255)                      aaGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGaa                                                                                                                           <  1:190940/151‑1 (MQ=255)                          cgttcgtCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAg                                                                                                                       <  2:423925/151‑1 (MQ=255)                              cgtCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGatc                                                                                                                   >  2:922402/1‑151 (MQ=255)                               gtCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGatcc                                                                                                                  >  1:727975/1‑151 (MQ=255)                                  cGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGatccatc                                                                                                               >  2:21049/1‑151 (MQ=255)                                             gAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGTCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGGCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTg                                                                                                     <  2:385608/150‑1 (MQ=255)                                             gAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTg                                                                                                     <  1:98977/150‑1 (MQ=255)                                                     tgttgACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCAtt                                                                                            >  1:318337/1‑151 (MQ=255)                                                       ttgACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTc                                                                                           <  1:757720/150‑1 (MQ=255)                                                       ttgACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTc                                                                                           <  2:366064/150‑1 (MQ=255)                                                        tgACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTcc                                                                                          >  2:291050/1‑150 (MQ=255)                                                        tgACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCt                                                                                         >  1:474562/1‑151 (MQ=255)                                                         gACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTg                                                                                        >  1:358891/1‑151 (MQ=255)                                                         gACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTg                                                                                        <  2:483053/151‑1 (MQ=255)                                                           ccAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGAc                                                                                      >  1:370231/1‑151 (MQ=255)                                                            cAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACg                                                                                     >  1:742780/1‑151 (MQ=255)                                                             aGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACgg                                                                                    >  1:382850/1‑151 (MQ=255)                                                                          tGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGt                                                                       >  1:822682/1‑151 (MQ=255)                                                                                cACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCCCCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTcg                                                                     <  2:104869/147‑1 (MQ=255)                                                                                cACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTcgc                                                                    <  1:645186/148‑1 (MQ=255)                                                                                  ctctGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGttg                                                                >  2:102019/1‑148 (MQ=255)                                                                                      ggCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGc                                                            <  1:313924/150‑1 (MQ=255)                                                                                           aGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAtttt                                                      >  1:205929/1‑151 (MQ=255)                                                                                               ggTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACtt                                                  >  1:640122/1‑151 (MQ=255)                                                                                                  ggCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAat                                               <  1:463236/151‑1 (MQ=255)                                                                                                     gAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAatgat                                            <  1:606420/151‑1 (MQ=255)                                                                                                         aaTCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGCTTTTACTTAAtgatga                                          >  1:423251/1‑149 (MQ=255)                                                                                                              cGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTga                                     <  2:815132/149‑1 (MQ=255)                                                                                                                gTATCGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTg                                      >  1:397287/1‑146 (MQ=255)                                                                                                                    cGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGc                              >  1:143199/1‑150 (MQ=255)                                                                                                                    cGAACTGACTGACGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCa                             <  1:942536/151‑1 (MQ=255)                                                                                                                        ctgactgaCGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGAcc                         >  2:653741/1‑151 (MQ=255)                                                                                                                        ctgactgaCGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGAcc                         >  1:420260/1‑151 (MQ=255)                                                                                                                          gactgacGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGAcc                         <  2:817786/149‑1 (MQ=255)                                                                                                                          gactgacGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCtt                       >  2:276158/1‑151 (MQ=255)                                                                                                                          gactgacGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCt                        <  1:339255/150‑1 (MQ=255)                                                                                                                             tgacGACCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTg                    >  1:106235/1‑151 (MQ=255)                                                                                                                              gacgacCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCt                     <  1:32015/149‑1 (MQ=255)                                                                                                                              gacgacCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTg                    <  1:931683/150‑1 (MQ=255)                                                                                                                                cgacCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCAGCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGcc                 <  2:835201/151‑1 (MQ=255)                                                                                                                                 gacCAGCTGAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCa                <  2:606696/151‑1 (MQ=255)                                                                                                                                         gAATATCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCaa         <  1:465371/150‑1 (MQ=255)                                                                                                                                            tatCACCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGgcg     >  2:990962/1‑151 (MQ=255)                                                                                                                                               ccccccGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCt   <  1:191681/148‑1 (MQ=255)                                                                                                                                                aCCCCGCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTg  >  2:787953/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGTGAGCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTG  >  NC_000913/1815443‑1815735

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