Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,882,432	C→G	6.9%	L77V (CTG→GTG)	PP_2537 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,882,432	0	C	G	6.9%	175.8 / 4.2	58	L77V (CTG→GTG)	PP_2537	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base G (2/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

CGAACTGGTGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGA  >  NC_002947/2882287‑2882575                                                                                                                                                  |                                                                                                                                                   cGAACTGGTGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCAGGGGCTGCAAAACCGAGATAAACCAAAAGGCCATCGACACGCCTGACGACCGCCAGGTGATGAAGAAGGGCCCCTACTACCAgg                                                                                                                                                 >  1:314776/1‑149 (MQ=255)     cTGGGGCAAAGGAACCCGCAACTGCCAATACGCACCGCACTGCTGGAAGCGTCGGGCATCATGGGCTGCAGCACCGAGATCACCCGCCCGGCCAGCGACAGGCCGGACGACCGCCAGGGGATGATGAAGGCCACGTACTACCTGGCCgg                                                                                                                                             <  3:31607/149‑1 (MQ=255)        gTGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCAGGTAGTACGTGGccttc                                                                                                                                            <  3:15928/147‑4 (MQ=255)        gTGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGccatccc                                                                                                                                          <  8:335130/149‑6 (MQ=255)         tGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCTTCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATg                                                                                                                                         <  1:317721/149‑1 (MQ=255)         tGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGccatcccg                                                                                                                                         >  3:109065/1‑143 (MQ=255)         tGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATg                                                                                                                                         <  7:224170/149‑1 (MQ=255)         agCAAATGAAACCGCAACAGCCAATACGCACCGCACTTATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATg                                                                                                                                         <  7:336969/148‑1 (MQ=255)            aaaTGAAACCGCAACTGCCAATACGCACCGCAGTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCCTCGACAGGCCTGCCGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGc                                                                                                                                       >  2:123400/1‑148 (MQ=255)               tGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAg                                                                                                                                   >  4:189213/1‑149 (MQ=255)               tGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAg                                                                                                                                   >  1:337409/1‑149 (MQ=255)                gAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGt                                                                                                                                  >  5:234384/1‑149 (MQ=255)                   aCCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCAc                                                                                                                               <  4:262294/149‑1 (MQ=255)                     cGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTg                                                                                                                             <  6:234384/149‑1 (MQ=255)                             ccAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCg                                                                                                                                              >  4:156975/1‑124 (MQ=255)                             ccAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCg                                                                                                                                              <  3:156975/124‑1 (MQ=255)                               aaTACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGcc                                                                                                                   >  2:39435/1‑149 (MQ=255)                               aaTACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCGGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGcc                                                                                                                   <  2:49097/149‑1 (MQ=255)                                      accgcacTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAgctgct                                                                                                            >  1:353816/1‑149 (MQ=255)                                      accgcacTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAgctgct                                                                                                            >  3:299875/1‑149 (MQ=255)                                           acTGATGGATGCTTCGGGCATCAGGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGa                                                                                                       <  2:314776/149‑1 (MQ=255)                                           acTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGa                                                                                                       <  4:192786/149‑1 (MQ=255)                                               aTGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑Cg                                                                                                  >  3:333350/1‑149 (MQ=255)                                                 ggAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGt                                                                                                 >  2:29574/1‑148 (MQ=255)                                                  gAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTc                                                                                                <  6:152312/148‑1 (MQ=255)                                                    aGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAgg                                                                                             <  2:353816/149‑1 (MQ=255)                                                     gCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGAGGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGgc                                                                                            <  6:390879/149‑1 (MQ=255)                                                       ttCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGgcgc                                                                                          >  5:142467/1‑149 (MQ=255)                                                                catGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCg                                                                                >  2:254710/1‑149 (MQ=255)                                                                         cATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGaa                                                                        <  5:155351/148‑1 (MQ=255)                                                                              ccGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAg                                                                  >  8:57438/1‑149 (MQ=255)                                                                                     cACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACg                                                                                                                      <  2:341083/92‑1 (MQ=255)                                                                                     cACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACg                                                                                                                      >  1:341083/1‑92 (MQ=255)                                                                                        ccGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGc                                                        >  6:36811/1‑149 (MQ=255)                                                                                        ccGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGc                                                        <  3:189213/149‑1 (MQ=255)                                                                                          gCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATATAAGCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGcc                                                       >  1:324956/1‑149 (MQ=255)                                                                                              ggCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCGCAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTggg                                                   <  1:254710/149‑1 (MQ=255)                                                                                                   tCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGg‑cgc                                            <  2:337409/149‑1 (MQ=255)                                                                                                       cAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAAc                                        >  7:127762/1‑149 (MQ=255)                                                                                                         ggCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACtt                                      <  8:127762/149‑1 (MQ=255)                                                                                                          gCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTg                                     >  7:233448/1‑149 (MQ=255)                                                                                                                        aGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCTTGGGCgg                       >  2:276821/1‑149 (MQ=255)                                                                                                                          gTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACTGCCAGCTGCTGCAGATT‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCTTGGGCGGtt                     <  1:123400/149‑1 (MQ=255)                                                                                                                            gatgatgaAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCTTGGGCGGttt                    >  7:356681/1‑148 (MQ=255)                                                                                                                              tgatgaAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCt                              <  7:364222/136‑1 (MQ=255)                                                                                                                              tgatgaAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCt                              >  8:364222/1‑136 (MQ=255)                                                                                                                                atgaAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGTTA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCCGCGAACTTGGCACCCTTGGGCGGTTTTaa                 <  1:29574/148‑1 (MQ=255)                                                                                                                                atgaAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCTTGGGCGGTTTTAATc               <  4:333350/149‑1 (MQ=255)                                                                                                                                 tgaAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGGGAAGCGTAGCACGGGTTGCAGGGGGGC‑CGCAACATTGGAACCGTTGGGCGGTTTAAAtgt              >  3:220704/1‑147 (MQ=255)                                                                                                                                   ccggCCAGGTAGTACGTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTc            >  7:39963/3‑149 (MQ=255)                                                                                                                                   aaGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGCTA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGACCTTGGCACCCTTGGGCGGTTTTAATCGTc            >  8:295136/1‑149 (MQ=255)                                                                                                                                    cggCCAGGTAGTACGTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCTTGGGCGGt                      <  3:94535/137‑1 (MQ=255)                                                                                                                                    cggCCAGGTAGTACGTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCTTGGGCGGt                      >  4:94535/2‑138 (MQ=255)                                                                                                                                    aGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCt           <  5:301136/149‑1 (MQ=255)                                                                                                                                     ggCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTg          >  4:237090/1‑149 (MQ=255)                                                                                                                                       ccACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGt         >  4:15487/1‑148 (MQ=255)                                                                                                                                          cGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAaca     <  5:36811/149‑1 (MQ=255)                                                                                                                                              ctacCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGC‑CGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGa  <  1:276821/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CGAACTGGTGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATA‑CGTCAGGCGCGCG‑AAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGG‑CGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGA  >  NC_002947/2882287‑2882575

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