Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,607,405	A→C	8.7%	R26S (AGA→AGC)	PP_4079 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,607,405	0	A	C	8.7%	155.7 / 2.4	56	R26S (AGA→AGC)	PP_4079	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base C (2/3);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTG  >  NC_002947/4607267‑4607552                                                                                                                                           |                                                                                                                                                      ccAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAgg                                                                                                                                             >  6:368127/1‑149 (MQ=255)        cgtcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCt                                                                                                                                      >  4:271562/1‑149 (MQ=255)          tcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGCGCCACCTTGGTGGCTCTgc                                                                                                                                    <  1:361324/149‑2 (MQ=255)          tcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCTTGGTGGCTCTgc                                                                                                                                    <  6:169104/149‑2 (MQ=255)          tcgATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTg                                                                                                                                     >  8:159364/1‑148 (MQ=255)            gATGAGCTCGATGGCCAAATCAGTCCCGTTGTCCGCTTCGGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCTTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGAc                                                                                                                                  <  1:245135/149‑1 (MQ=255)            gATGACCACGATGGCCAAATCATTCCCGTTGTCCGCTTCCGGTGGCGGTCATGATCCGCCTGGCAGTCTGGATTGCGCCTCTGCTTGTCGTGCTGGTCGGACTGCCAGGTTGCCCGTCCTTGAGCATCTCCACCAAGGTGGCGCTGGAc                                                                                                                                  >  7:142948/1‑149 (MQ=255)            gATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGa                                                                                                                                   <  5:368127/148‑1 (MQ=255)            gATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGCTTGCACGTCCTAGAGCATAGCCCCCAAGGTGGCGCTGGAc                                                                                                                                  >  7:137772/1‑149 (MQ=255)             aTGACCACGATGGCCAAATCAGTCCCGTTGT‑CGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGCGCCACCTTGGTGGCTCTgctcaa                                                                                                                                <  5:324315/149‑6 (MQ=255)                     gATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTa                                                                                                                        <  5:308966/149‑1 (MQ=255)                          ccAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAAt                                                                                                                   >  3:234436/1‑149 (MQ=255)                              aTCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCggg                                                                                                               <  4:197604/149‑1 (MQ=255)                              aTCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCggg                                                                                                               >  7:274120/1‑149 (MQ=255)                               tCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGgc                                                                                                              >  5:72196/1‑149 (MQ=255)                                   tCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGCGCCACCTTGGTGGCTCTgctcaa                                                                                                                                <  7:324378/128‑6 (MQ=255)                                   tCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCACGTTGCACGTCCTTGAGCACCGACACCTTGGTGGCTCTgctcaa                                                                                                                                >  8:324378/1‑123 (MQ=255)                                       gTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAg                                                                                                      >  8:163117/1‑149 (MQ=255)                                          gTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGa                                                                                                    <  4:234436/148‑1 (MQ=255)                                           tCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGAc                                                                                                   >  6:122533/1‑148 (MQ=255)                                             cGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGCGCCGCCCAAGGAcaac                                                                                                >  5:4224/1‑149 (MQ=255)                                              gCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGAcaacc                                                                                               <  3:271562/149‑1 (MQ=255)                                                   cGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGAcaacccaacc                                                                                          <  3:83427/149‑1 (MQ=255)                                                      gggCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACaacccaaccc                                                                                         >  8:310883/1‑147 (MQ=255)                                                       ggCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAg                                                                                      >  6:285783/1‑149 (MQ=255)                                                          ggTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTcg                                                                                   >  6:187823/1‑149 (MQ=255)                                                            tCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTcgcg                                                                                 <  5:122533/149‑1 (MQ=255)                                                             cATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGt                                                                                >  2:39855/1‑149 (MQ=255)                                                               tGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTc                                                                              <  7:310883/149‑1 (MQ=255)                                                               tGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTc                                                                              <  8:274120/149‑1 (MQ=255)                                                                     gccTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACAACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCAt                                                                         >  5:359757/1‑148 (MQ=255)                                                                          gCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGt                                                                   <  7:159364/149‑1 (MQ=255)                                                                                        gcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGt                                                                                <  4:251175/122‑1 (MQ=255)                                                                                        gcCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGt                                                                                >  3:251175/1‑122 (MQ=255)                                                                                         ccTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCaa                                                   <  7:163117/149‑1 (MQ=255)                                                                                             tgctggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAg                                                >  8:39882/1‑148 (MQ=255)                                                                                                 ggtcgtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGTGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAGCAAcc                                           >  2:28058/1‑149 (MQ=255)                                                                                                     gtgctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAGCAACCCGaa                                       <  7:39882/149‑1 (MQ=255)                                                                                                        ctggtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGCCCACTC‑GTTACCAATCGGGCCGCCCAAGGCCAACCCAACCCAAGTCGCGTTCAGCCTCACCGTC‑AGCCCAACGCTCAACCGCAAGCCGAACAg                                    >  8:284727/1‑149 (MQ=255)                                                                                                           gtcgTACTGCCAGGTTGCACTTCCTTGCGCATAGCCACAAAGGTGGCGCTGGACCACTC‑GTTACCACTCGGGCCGCCCAATGACAACCCAACCCAAGTCGCTTTCATCATCAACGTC‑AGCCCAACGCTCAACAGCAACCCGAACAgtg                                  <  8:142948/148‑1 (MQ=255)                                                                                                           gtcgGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAGCAACCCGAACAgtgt                                 >  2:187137/1‑149 (MQ=255)                                                                                                              ggACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTCAAGCCCAACGCTCAACAGCAACCCGAACAGTGTc                                <  1:39855/148‑1 (MQ=255)                                                                                                                aCTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCAt                                                                         <  5:69767/105‑1 (MQ=255)                                                                                                                aCTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCAt                                                                         >  6:69767/1‑105 (MQ=255)                                                                                                                  tGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAGCAACCCGAAc                                      >  7:344800/1‑137 (MQ=255)                                                                                                                  tGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAGCAACCCGAAc                                      <  8:344800/137‑1 (MQ=255)                                                                                                                    ccAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTg                          <  5:187823/147‑1 (MQ=255)                                                                                                                           gCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTCGGTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCAt                  >  1:78957/1‑149 (MQ=255)                                                                                                                           gCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAAc                                                  <  5:177138/116‑1 (MQ=255)                                                                                                                           gCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAAc                                                  >  6:177138/1‑116 (MQ=255)                                                                                                                            cACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATcc                <  2:78957/149‑1 (MQ=255)                                                                                                                             aCGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATccc               <  5:285783/149‑1 (MQ=255)                                                                                                                                tCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTga            <  3:305110/149‑1 (MQ=255)                                                                                                                                  cTTTAGCAGAGCCACCTTGGTGGCGCTGGAAGACTC‑GTTACCAATCGTGGCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCATCGGTCAACAGCATCCCGAACAGTGTCGACGTGGCGGCCATCCCTgaga          <  6:359757/149‑1 (MQ=255)                                                                                                                                     gAGCAGAGACACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCAGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCCACAGCAACCCGAACAGTGTCGACGTGGCGGCCATACCTGAGAGCa       >  5:208774/1‑149 (MQ=255)                                                                                                                                        cAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATcc    >  7:53013/1‑149 (MQ=255)                                                                                                                                          gagCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTg  >  3:163352/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                      CCAGCTTCGTCGATGACCACGATGGCCAAATCAGTCCCGTTGTCCGCTTCCGGGGGCGGTCATGAGCCGCCTGGCCGTCTGGATTGCGCCTCTGCTGGTCGTGCTGGTCGGACTGCCAGGTTGCACGTCCTTGAGCAGAGCCACCAAGGTGGCGCTGGACCACTC‑GTTACCAATCGGGCCGCCCAAGGACAACCCAACCCAAGTCGCGTTCAGCATCAACGTC‑AGCCCAACGCTCAACAGCAACCCGAACAGTGTCGACGTGGCGGCCATCCCTGAGAGCATCCTG  >  NC_002947/4607267‑4607552

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: ‑