Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,980,512	G→T	6.5%	Q83K (CAA→AAA)	USA300HOU_1839 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,980,512	0	G	T	6.5%	211.3 / 3.2	61	Q83K (CAA→AAA)	USA300HOU_1839	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/33);  new base T (2/2);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

TTGATTTCCATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGACTCTAAATTTAAATATAT  >  CP000730/1980378‑1980649                                                                                                                                       |                                                                                                                                          ttGATTTCCATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGACGTAAATTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCt                                                                                                                                     <  1:18980/141‑1 (MQ=255)   gATTTCCATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCttt                                                                                                                                   >  1:225140/1‑141 (MQ=255)   gATTTCCATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCttt                                                                                                                                   >  1:250385/1‑141 (MQ=255)      ttCCATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCa                                                                                                                                >  2:127047/1‑141 (MQ=255)        ccATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTTCTCGCTTTTCAAt                                                                                                                              >  1:234334/1‑141 (MQ=255)        ccATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTTCTCGCTTTTCAAt                                                                                                                              >  1:234365/1‑141 (MQ=255)         cATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAAtt                                                                                                                             >  2:77720/1‑141 (MQ=255)           ttAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCATCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAAtttt                                                                                                                           <  2:442653/141‑1 (MQ=255)               tttGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCat                                                                                                                       >  2:90929/1‑141 (MQ=255)               tttGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTAACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCat                                                                                                                       >  1:221507/1‑141 (MQ=255)                            ttCTAAATCAAGATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTGTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTATTCAATTTTTCATTAATTAACGTTTc                                                                                                          >  2:268284/1‑141 (MQ=255)                                   tCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGc                                                                                                   <  2:20010/141‑1 (MQ=255)                                   tCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGc                                                                                                   <  2:286942/141‑1 (MQ=255)                                    cAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCa                                                                                                  >  2:66794/1‑141 (MQ=255)                                    cAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCa                                                                                                  >  1:325007/1‑141 (MQ=255)                                       tatCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAAt                                                                                               >  2:175239/1‑141 (MQ=255)                                       tatCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAAt                                                                                               >  2:36389/1‑141 (MQ=255)                                       tatCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAAt                                                                                               >  2:36388/1‑141 (MQ=255)                                               tGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCa                                                                                       <  1:271141/141‑1 (MQ=255)                                                 agaAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAAt                                                                                     <  2:181048/141‑1 (MQ=255)                                                     gTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCTCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACt                                                                                 >  1:445986/1‑141 (MQ=255)                                                         tttGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGc                                                                             >  2:246745/1‑141 (MQ=255)                                                            gTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCaa                                                                          <  1:44914/141‑1 (MQ=255)                                                                tATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAAtttt                                                                      <  1:328406/141‑1 (MQ=255)                                                                           ggtCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAAttt                                                           <  2:180299/139‑1 (MQ=255)                                                                           gATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTCAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAAttt                                                           <  2:14340/141‑1 (MQ=255)                                                                           gATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAAttt                                                           <  2:139425/141‑1 (MQ=255)                                                                             tCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATttct                                                         <  1:84843/141‑1 (MQ=255)                                                                             tCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATttct                                                         <  2:368886/141‑1 (MQ=255)                                                                             tCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATttct                                                         <  2:49058/141‑1 (MQ=255)                                                                             tCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATttct                                                         <  1:428226/141‑1 (MQ=255)                                                                             tCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATttct                                                         <  1:152077/141‑1 (MQ=255)                                                                              cTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATttctt                                                        >  1:162117/1‑141 (MQ=255)                                                                                   gTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGt                                                   >  1:227149/1‑141 (MQ=255)                                                                                    tCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTc                                                  <  1:246522/141‑1 (MQ=255)                                                                                     cAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCt                                                 >  1:147380/1‑141 (MQ=255)                                                                                      aaCTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCtt                                                <  1:61208/141‑1 (MQ=255)                                                                                               aaCTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTAtt                                       <  1:443696/141‑1 (MQ=255)                                                                                                      gTTGAATCATAGATTCTATTGCTTGCTGTGATTTCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTg                                <  2:234334/141‑1 (MQ=255)                                                                                                      gTTGAATCATAGATTCTATTGCTTGCTGTGATTTCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTg                                <  2:234365/141‑1 (MQ=255)                                                                                                         gAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGgtc                             <  1:436350/141‑2 (MQ=255)                                                                                                            tCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTaa                          <  1:36389/141‑1 (MQ=255)                                                                                                            tCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTaa                          <  1:36388/141‑1 (MQ=255)                                                                                                            tCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTaa                          <  1:175239/141‑1 (MQ=255)                                                                                                              aTAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAAt                        <  2:225140/141‑1 (MQ=255)                                                                                                              aTAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAAt                        <  2:250385/141‑1 (MQ=255)                                                                                                              aTAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAAt                        <  2:399984/141‑1 (MQ=255)                                                                                                              aTAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAAt                        <  2:305760/141‑1 (MQ=255)                                                                                                               tAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATg                       >  2:174466/1‑141 (MQ=255)                                                                                                                  tttCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGa                    <  2:369635/140‑1 (MQ=255)                                                                                                                  aTTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGa                    <  2:172193/141‑1 (MQ=255)                                                                                                                    tCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGGAAATGCTTTTTAGTGACGCGGCGAATCTTGGTTGCAATTTCTTCGGGGCTTTTCGGTATTAAACCTGGTGTAAATGGGAcg                  >  1:344940/1‑140 (MQ=255)                                                                                                                    tCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAACTGCTCTTCAATAACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGAct                  >  1:80319/1‑141 (MQ=255)                                                                                                                    tCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGAct                  >  1:424834/1‑141 (MQ=255)                                                                                                                      tATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGActct                <  1:249119/141‑1 (MQ=255)                                                                                                                         tgcttgctGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGACTCTaaa             <  2:227149/141‑1 (MQ=255)                                                                                                                         tgcttgctGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGACTCTaaa             <  1:198708/141‑1 (MQ=255)                                                                                                                                tgtgATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCCATCCCTTGGCCACTTCTAGTTGCAATTTCTTCGAGCCTAATAGCTCACAAACAAACTGCAAACGGGACTCTAAATTTAAat      >  1:92672/1‑141 (MQ=255)                                                                                                                                tgtgATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGACTCTAAATTTAAat      >  2:177039/1‑141 (MQ=255)                                                                                                                                  tgATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGACTCTAAATTTAAatat    >  1:272533/1‑141 (MQ=255)                                                                                                                                    aTTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGACTCTAAATTTAAatatat  <  2:371727/141‑1 (MQ=255)                                                                                                                                       |                                                                                                                                          TTGATTTCCATTAGTTTGTAATACTTGTTCTAAATCAATATCGATTTGAGAAGTAATTTGTTTTATTGACAATTGATCTTTCGTCAACTTTTGTAACTGCTGTTGAATCATAGATTCTATTGCTTGCTGTGATTGCTCGCTTTTCAATTTTTCATTAATTAAAGTTTCTGTAAGCAAATGCTCTTCAATCACTTGGCCAATTTTAGTTGCAATTTCTTCGCGTCTTTTCGGTATTAAACCTGGTGTAAATGGGACTCTAAATTTAAATATAT  >  CP000730/1980378‑1980649

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