Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	249,079	C→T	7.0%	intergenic (+27/+22)	PP_0198 → / ← qmcA	putative transporter/membrane protease family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	249,079	0	C	T	7.0%	165.7 / 2.6	57	intergenic (+27/+22)	PP_0198/qmcA	putative transporter/membrane protease family protein
Reads supporting (aligned to +/- strand):  ref base C (30/23);  new base T (2/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GCGCATGCTGCGCACGCCCAAACAGCAGAAGCGTTTGAACCGCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCC  >  NC_002947/248944‑249227                                                                                                                                         |                                                                                                                                                      gcgcATGCTGCGCACGCCCAAACAGCAGAAGCGTTTGAACCGCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTg                                                                                                                                         >  6:100384/1‑149 (MQ=255) gcgcATGCTGCGCACGCCCAAACAGCAGAAGCGTTTGAACCGCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTg                                                                                                                                         >  3:147807/1‑149 (MQ=255)  cgcATGCTGCGCACGCCCAAACAGCAGAAGCGTTTGAACCGCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTgg                                                                                                                                        <  4:147807/149‑1 (MQ=255)                gCCCAAACAGCAGAAGCGTTTGAACCGCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTtc                                                                                                                          >  5:69303/1‑149 (MQ=255)                gCCCAAACAGCAGAAGCGTTTGAACCGCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTtc                                                                                                                          >  6:80640/1‑149 (MQ=255)                   cAAACAGCAGAAGCGTTTGAACCGCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTcc                                                                                                                       <  3:224410/149‑1 (MQ=255)                    aaaCAGCAGAAGCGTTTGAACCGCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCt                                                                                                                      >  8:13517/1‑149 (MQ=255)                       cagcagAAGCGTTTGAACCGCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGt                                                                                                                   >  5:221233/1‑149 (MQ=255)                         gcagAAGCGTTTGAACCGCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTcc                                                                                                                 <  5:80640/149‑1 (MQ=255)                             aaGCGTTTGAACCGCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTccacc                                                                                                              >  3:107013/1‑148 (MQ=255)                                     gAACCGCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCtt                                                                                                      >  2:197210/1‑148 (MQ=255)                                        ccGCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCTATGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGTAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTcag                                                                                                   >  1:242284/1‑149 (MQ=255)                                          gCAGCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTcagcag                                                                                                <  5:100384/149‑1 (MQ=255)                                          gCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTcag                                                                                                   >  2:217130/1‑146 (MQ=255)                                               tttGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCt                                                                                            >  1:312916/1‑148 (MQ=255)                                                  gCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATc                                                                                         <  6:221233/148‑1 (MQ=255)                                                      gggtgttCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCg                                                                                        <  6:225789/142‑1 (MQ=255)                                                      ggTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCg                                                                                        >  5:225789/1‑145 (MQ=255)                                                           ttCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCggg                                                                                >  5:60097/1‑148 (MQ=255)                                                              aTTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATc                                                                                             >  2:7239/1‑132 (MQ=255)                                                              aTTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATc                                                                                             <  1:7239/132‑1 (MQ=255)                                                              aTTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAg                                                                             <  4:107013/148‑1 (MQ=255)                                                              aTTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGa                                                                            >  7:70871/1‑149 (MQ=255)                                                                   tGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATg                                                                       <  2:321287/149‑1 (MQ=255)                                                                    gCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTcag                                                                                                   <  6:158995/120‑1 (MQ=255)                                                                    gCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTcag                                                                                                   >  5:158995/1‑120 (MQ=255)                                                                     cAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATgg                                                                      <  6:69303/148‑1 (MQ=255)                                                                           ccTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGct                                                               <  1:217130/149‑1 (MQ=255)                                                                            cTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGctc                                                              <  7:13517/149‑1 (MQ=255)                                                                               ccTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCtt                                                            <  6:317971/148‑1 (MQ=255)                                                                                    cgACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAACCCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTc                                                                                                     <  1:233214/103‑1 (MQ=255)                                                                                    cgACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACACCTTc                                                                                                     >  2:233214/1‑102 (MQ=255)                                                                                     gACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTAt                                                                                          <  3:242406/112‑1 (MQ=255)                                                                                     gACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTAt                                                                                          >  4:242406/1‑112 (MQ=255)                                                                                           tGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAAt                                                >  5:184734/1‑149 (MQ=255)                                                                                              gccgcGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGt                                             >  2:77968/1‑149 (MQ=255)                                                                                                ctgtatgcGGTCTGTCAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTcc                                           >  3:161091/8‑149 (MQ=255)                                                                                                   gcgcCGGTCTGACAGACCGCATACAGGCCATTCGCGGGTAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCa                                        <  2:242284/149‑1 (MQ=255)                                                                                                     gcCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAgt                                      <  6:60097/149‑1 (MQ=255)                                                                                                       cGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAgtgt                                    <  1:77968/149‑1 (MQ=255)                                                                                                         gTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCt                                  <  2:312916/149‑1 (MQ=255)                                                                                                            tGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCa                               >  2:261406/1‑149 (MQ=255)                                                                                                              aCAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCat                              <  8:70871/148‑1 (MQ=255)                                                                                                                   ccGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGCCGATGGTCGAGCTCTTGTCCCAGCCAATCGTCCCCAGTGTCTGCATAAAgcg                        >  6:230898/1‑149 (MQ=255)                                                                                                                   ccGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAgcg                        <  8:149590/149‑1 (MQ=255)                                                                                                                    cGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAgcg                        >  6:158381/1‑148 (MQ=255)                                                                                                                    cGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAgcg                        >  3:193146/1‑148 (MQ=255)                                                                                                                     agtgatgggCCAATCGCGGGTTTG‑CCCGCGAATGGCCTATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCatat                            >  7:285326/8‑143 (MQ=38)                                                                                                                     agtgatgggCCAATCGCGGGTTTG‑CCCGCGAATGGCCTATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCatat                            <  8:285326/136‑1 (MQ=38)                                                                                                                           aGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATc                                               >  4:141173/1‑118 (MQ=255)                                                                                                                           aGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATc                                               <  3:141173/118‑1 (MQ=255)                                                                                                                           aGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGc                 >  2:94355/1‑148 (MQ=255)                                                                                                                           aGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCa                >  3:286432/1‑149 (MQ=255)                                                                                                                             gCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATc              <  1:197210/149‑1 (MQ=255)                                                                                                                               cATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGc            >  3:120455/1‑149 (MQ=255)                                                                                                                                        gCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTc   >  8:216985/1‑149 (MQ=255)                                                                                                                                        gCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTc   >  1:176693/1‑149 (MQ=255)                                                                                                                                         cAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTcc  <  4:161091/149‑1 (MQ=255)                                                                                                                                         |                                                                                                                                                      GCGCATGCTGCGCACGCCCAAACAGCAGAAGCGTTTGAACCGCACCTTTGCCGGGTTGTTCATTGGTGCAGCCACCTTCCTCGCGACCT‑TGCGCCGCGCGCCGGTCTGACAGACCGCATACAGGCCATTCGCGGGCAAA‑CCCGCGATTGGCCCATCACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCC  >  NC_002947/248944‑249227

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