Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,065,989	C→G	9.1%	V1025V (GTG→GTC)	mdtC ←	multidrug transporter membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,065,989	0	C	G	9.1%	163.9 / 5.0	54	V1025V (GTG→GTC)	mdtC	multidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base C (23/25);  new base G (2/3);  total (25/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 = 1.00e+00

TCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTT  >  NC_002947/4065845‑4066131                                                                                                                                                 |                                                                                                                                                 tCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTGCGTaccgacgc                                                                                                                                              >  7:193875/1‑149 (MQ=255)       tACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTGCGTACCGAcgcc‑gctct                                                                                                                                       <  8:193875/149‑4 (MQ=255)          aGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAAc                                                                                                                                    >  2:381397/1‑149 (MQ=255)          aGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGAGGCGTCGGTACGCACGCC‑TCGCCAAc                                                                                                                                    >  6:257830/1‑149 (MQ=255)            gCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACgg                                                                                                                                  >  3:337246/1‑149 (MQ=255)              aCTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACgg                                                                                                                                  <  6:199379/147‑1 (MQ=255)              aCTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACgg                                                                                                                                  >  5:199379/1‑147 (MQ=255)              aCTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGtt                                                                                                                                <  5:133426/149‑1 (MQ=255)                   gCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAAc                                                                                                                           <  7:170734/149‑1 (MQ=255)                      aGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGt                                                                                                                        >  8:255805/1‑149 (MQ=255)                       ggCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTg                                                                                                                       >  1:70031/1‑149 (MQ=255)                                 aaCCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTGCGTACCGAcgcc‑gct                                                                                                                                         <  3:103436/121‑2 (MQ=255)                                 aaCCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTGCGTACCGAcgcc‑gct                                                                                                                                         >  4:103436/1‑120 (MQ=255)                                 aaCCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGGCATAGCGGGGTTTCCAGAGCGGCGGCGGTACGCACGCCTGCGCCACCGTTTCAACCGTTGgcgcggcgt                                                                                                              <  7:161852/149‑6 (MQ=38)                                   cccGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTc                                                                                                           <  2:70031/149‑1 (MQ=255)                                      ggCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcag                                                                                                        >  1:42516/1‑149 (MQ=255)                                      ggCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcag                                                                                                        <  4:337246/149‑1 (MQ=255)                                       gCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcagg                                                                                                       <  3:156037/149‑1 (MQ=255)                                          ccGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGGCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtac                                                                                                    <  6:417330/149‑1 (MQ=255)                                           cGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtaca                                                                                                   <  1:381397/149‑1 (MQ=255)                                            gAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacag                                                                                                  <  5:200794/149‑1 (MQ=255)                                            gAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacag                                                                                                  <  5:179648/149‑1 (MQ=255)                                             aGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCcaggtacagg                                                                                                 >  2:292239/1‑149 (MQ=255)                                                    gcgAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGtt                                                                                                                                >  5:342199/1‑111 (MQ=255)                                                    gcgAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGtt                                                                                                                                <  6:342199/111‑1 (MQ=255)                                                    gcgAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAacca                                                                                          >  5:315077/1‑149 (MQ=255)                                                       aGCCACGCCCACGGTTCAGTGGTTGCCGGGCACGGTGAAGTGGGGTCTGGCCAAAAGTCATACCGGGGTTTCCAGAGCGGCGTCGGTACGAACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCg                                                                                       <  5:257830/149‑1 (MQ=255)                                                          cacggccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGcg                                                                                    >  3:250576/1‑149 (MQ=255)                                                           acggccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGcgt                                                                                   >  7:197030/1‑149 (MQ=255)                                                               ccacggGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGta                                                                               >  3:23392/1‑149 (MQ=255)                                                                         gTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGa                                                                     <  4:250576/149‑1 (MQ=255)                                                                               gCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGcgtc                                                                                  >  4:267834/1‑130 (MQ=255)                                                                               gCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGcgtc                                                                                  <  3:267834/130‑1 (MQ=255)                                                                                 aGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTgg                                                             <  7:255805/149‑1 (MQ=255)                                                                                   ggCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCt                                                           <  7:151760/149‑1 (MQ=255)                                                                                      aCGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGaa                                                        <  5:114925/149‑1 (MQ=255)                                                                                        ggTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAaca                                                      >  3:205669/1‑149 (MQ=255)                                                                                            aaGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAg                                                 <  4:23392/149‑1 (MQ=255)                                                                                                   gTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCAcc                                           >  5:369450/1‑148 (MQ=255)                                                                                                   gTCTGGGCAAAAGTCATAGAGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCAcc                                           >  4:387345/1‑148 (MQ=255)                                                                                                         gCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCACCgatg                                       <  4:298699/146‑1 (MQ=255)                                                                                                         gCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCACCgatg                                       >  3:298699/1‑146 (MQ=255)                                                                                                         gCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGACGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCACCGatgatg                                    >  1:44794/1‑149 (MQ=255)                                                                                                           aaaaGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCACCGATGATGGt                                  >  6:54190/1‑149 (MQ=255)                                                                                                             aaGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCACCGATGATGGt                                  >  4:385437/1‑147 (MQ=255)                                                                                                               gTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCACCGATGATGGTCAgg                              <  5:54190/149‑1 (MQ=255)                                                                                                               gTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCACCGATGATGGTCAgg                              <  4:205669/149‑1 (MQ=255)                                                                                                                     gCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCACCGATGATGGTCAGGCCCAg                         <  6:369450/148‑1 (MQ=255)                                                                                                                     gCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCACCGATGATGGTCAGGCCCAGc                        >  1:155360/1‑149 (MQ=255)                                                                                                                                ggcgcggcgTGCGTACCGACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGta                                                                               <  2:152032/81‑1 (MQ=38)                                                                                                                                ggcgcggcgTGCGTACCGACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGta                                                                               >  1:152032/4‑84 (MQ=38)                                                                                                                                    cggcgTCGGTACGCACGGC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCt         <  8:197030/149‑1 (MQ=255)                                                                                                                                       cgTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACACCCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTgc       <  1:292239/149‑1 (MQ=255)                                                                                                                                        gTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTgcgc     <  3:387345/149‑1 (MQ=255)                                                                                                                                           ggTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCtt  >  5:369553/1‑149 (MQ=255)                                                                                                                                                gcacgcC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCa                                                  <  7:37455/96‑1 (MQ=255)                                                                                                                                                gcacgcC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCa                                                  >  8:37455/1‑96 (MQ=255)                                                                                                                                                 |                                                                                                                                                 TCAGGGTACAGGCACTGAGCAAGGCCACGCACAACCCGGCACCGAGCAGGCGCGAGCCACGGCCACGGGTCAGTGCTTGCAGGGCACGGTGAAGTGGGGTCTGGGCAAAAGTCATAGCGGGGTTTCCAGAGCGGCGTCGGTACGCACGCC‑GCGCCAACGGTTGAAACGGTGGCGCAGGCGGTCCAGGTACAGGTAAACCACCGGCGTCGTATAAAGGGTGAGGATCTGGCTGAACA‑CCAGGCCACCGATGATGGTCAGGCCCAGCGGCTGACGCATTTCTGCGCCTT  >  NC_002947/4065845‑4066131

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