Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,589,518	G→A	7.1%	L13L (CTG→CTA)	nicP‑I →	porin‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,589,518	0	G	A	7.1%	164.0 / 2.7	57	L13L (CTG→CTA)	nicP‑I	porin‑like protein
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base A (2/2);  total (26/31)
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.57e-01

CGCTGCTGAGGCTGCGCCAGCCACGTGCCCGGGAGGTGCTGCGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCGGGGCCAACAAACAGTCCAAGGCCGAGGAGTGGGCCCAGGGTT  >  NC_002947/3589370‑3589666                                                                                                                                                     |                                                                                                                                                      cgcTGCTGAGGCTGCGCCAGCCACGTGCCCGGGAGGTGCTGCGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTg                                                                                                                                                       <  7:9605/149‑1 (MQ=255)    tgctgAGGCTGCGCCAGCCACGTGCCCGGGAGGTGCTGCGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCt                                                                                                                                                     >  4:15391/1‑148 (MQ=255)      ctgAGGCTGCGCCAGCCACGTGCCCGGGAGGTGCTGCGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTACTcag                                                                                                                                                  >  5:269044/1‑149 (MQ=255)                  cagccaCGTGCCCGGGAGGTGCTGCGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTACTCAGCCAGCCActgct                                                                                                                                      >  6:275676/1‑149 (MQ=255)                   agccaCGTGCCCGGGAGGTGCTGCGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCActgctg                                                                                                                                     <  1:41014/149‑1 (MQ=255)                    gccaCGTGCCCGGGAGGTGCTGCGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTgctgg                                                                                                                                    <  6:56380/149‑1 (MQ=255)                     ccaCGTGCCCGGGAGGTGCTGCGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCGTAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTgctggc                                                                                                                                   >  1:292249/1‑149 (MQ=255)                          tGCCCGGGAGGTGCTGCGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGcc                                                                                                                               >  8:133890/1‑148 (MQ=255)                          tGCCCGGGAGGTGCTGCGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGcc                                                                                                                               <  7:133890/148‑1 (MQ=255)                             ccGGGAGGTGCTGCGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAg                                                                                                                           <  4:155195/149‑1 (MQ=255)                               gggAGGTGCTGCGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTACTCAGCCAGCCACTGCTGGCTGCCGAAggc                                                                                                                         <  6:269044/149‑1 (MQ=255)                                  aGGTGCTGCGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACGCTGCTCAGCCAGCCACTGCTGGCTGCCGAAggcggc                                                                                                                      <  4:71018/149‑1 (MQ=255)                                        tgcGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCATGGGCCTGACACTGCTCAGCCAGCCACTGCTGGATGCCGAAGGCGGCttcttc                                                                                                                >  3:103441/1‑149 (MQ=255)                                                  cccATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCActg                                                                                                                                        >  2:176909/1‑115 (MQ=255)                                                  cccATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCg                                                                                                      >  2:127592/1‑148 (MQ=255)                                                  cccATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCAGGGCCTGACACTGCTCAGCCAGCCActg                                                                                                                                        <  1:176909/115‑1 (MQ=255)                                                    cATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTACTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAg                                                                                                   <  5:275676/149‑1 (MQ=255)                                                    cATGAACTGCCGGGTTCCCACCGCTTACCCCATAAAACCACGATCAACAGACCTTCTGCATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAg                                                                                                   <  8:101621/149‑1 (MQ=255)                                                     aTGAACTGCCGTGTTCCCAACTCTTTAGCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAgg                                                                                                  <  7:157738/149‑1 (MQ=255)                                                     aTGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACCACCCTTGCGTTGCCCTGTGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAgg                                                                                                  <  7:199861/149‑1 (MQ=255)                                                          cTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCAcc                                                                                             <  1:127592/149‑1 (MQ=255)                                                              cGTTTTCCCAACGCTTAACCCATAAAACCAAGATCAAGAGAGTTTCGCAATGAACCACCCCCCTTGCGTTGCCCTGGGCCTGACTCTGCCCACCCCCCCACTGGTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGa                                                                                         <  2:152096/149‑1 (MQ=255)                                                                tgtTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGCACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCttcttc                                                                                                                >  8:141772/1‑125 (MQ=255)                                                                tgtTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCttcttc                                                                                                                <  7:141772/125‑1 (MQ=255)                                                                         cGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAac                                                                              <  1:286157/149‑1 (MQ=255)                                                                           cTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑Ga                                                                                                          >  7:156202/1‑119 (MQ=255)                                                                           cTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑Ga                                                                                                          <  8:156202/119‑1 (MQ=255)                                                                              tACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAactac                                                                           <  5:210527/147‑1 (MQ=255)                                                                               tcaaaTAAAAACAACATCAACACAGCTTAGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGTCACTGCTCACCCAGCCACTGCTGGCTTCCGAAGGCGGCTTCTTCGAGTGACTCGAAGGCCACCCTGAGTGCCCGTAactacta                                                                         <  1:47659/145‑1 (MQ=255)                                                                                 ccATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑Ga                                                                                                          >  1:97195/1‑113 (MQ=255)                                                                                 ccATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑Ga                                                                                                          <  2:97195/113‑1 (MQ=255)                                                                                   aTAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAg                                                                                                   >  4:168581/1‑118 (MQ=255)                                                                                   aTAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAg                                                                                                   <  3:168581/118‑1 (MQ=255)                                                                                   aTAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGACAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAactacta                                                                         >  2:47659/1‑144 (MQ=255)                                                                                           aaGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGa                                                                                         <  4:84631/120‑1 (MQ=255)                                                                                           aaGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGa                                                                                         >  3:84631/1‑120 (MQ=255)                                                                                                 aaGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTcc                                                      <  2:292249/149‑1 (MQ=255)                                                                                                         ttCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGt                                              >  3:195899/1‑149 (MQ=255)                                                                                                             gAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCtgg                                          >  8:197024/1‑149 (MQ=255)                                                                                                             gAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCggg                                          >  2:182005/1‑149 (MQ=255)                                                                                                              aaTGATCGACACCGCTTGCGTTGCCCTTGGCCTGACAGTTCTCAGCCTGGCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCg                                            <  4:103441/146‑1 (MQ=255)                                                                                                              aaTGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCgggg                                         <  6:165180/149‑1 (MQ=255)                                                                                                              aaTGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCgggg                                         >  7:168586/1‑149 (MQ=255)                                                                                                                      acacCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCGGGGCCaacaaa                                 >  5:20046/1‑149 (MQ=255)                                                                                                                           ccTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGa                                                                                         >  6:138515/1‑88 (MQ=255)                                                                                                                           ccTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGa                                                                                         <  5:138515/88‑1 (MQ=255)                                                                                                                                 gTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCggg                                          <  5:80596/129‑1 (MQ=255)                                                                                                                                 gTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCggg                                          >  6:80596/1‑129 (MQ=255)                                                                                                                                 gTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCGGGGCCAACAAACAGTCCAAGGc                      >  7:273318/1‑149 (MQ=255)                                                                                                                                  ttGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGACTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCGGGGCCAACAAACAGTCCAAGGcc                     >  5:184088/1‑149 (MQ=255)                                                                                                                                    gCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCGGGGCCAACAAACAGTCCAAGGCCg                    <  8:273318/148‑1 (MQ=255)                                                                                                                                          ggCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCGGGGCCAACAAACAGTCCAAGGCCGAGGAGTg             <  4:195899/149‑1 (MQ=255)                                                                                                                                              tGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCggg                                          >  2:289721/1‑116 (MQ=255)                                                                                                                                              tGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCggg                                          <  1:289721/116‑1 (MQ=255)                                                                                                                                               gACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCGGGGCCAACAAACAGTCCAAGGCCGAGGAGTGGGcc         >  3:142010/1‑148 (MQ=255)                                                                                                                                               gACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCGGGGCCAACAAACAGTCCAAGGCCGAGGAGTGGGcc         >  3:137361/1‑148 (MQ=255)                                                                                                                                                acacTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCGGGGCCAACAAACAGTCCAAGGCCGAGGAGTGGGcc         <  3:15391/147‑1 (MQ=255)                                                                                                                                                  acTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAT‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCGGGGCCAACAAACAGTCCAAGGCCGAGGAGTGGGCCCAgg     <  6:20046/149‑1 (MQ=255)                                                                                                                                                     gCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCGGGGCCAACAAACAGTCCAAGGCCGAGGAGTGGGCCCAGGGtt  >  5:195028/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CGCTGCTGAGGCTGCGCCAGCCACGTGCCCGGGAGGTGCTGCGCCGTACCCCATGAACTGCCGTGTTCCCAACGCTTTACCCATAAAAACAAGATCAAGAGAGCTTCGGAATGAACCACACCCCTTGCGTTGCCCTGGGCCTGACACTGCTCAGCCAGCCACTGCTGGCTGCCGAAGGCGGCTTCTTCGAG‑GACTCGAAGGCCACCCTGAGTGCCCGTAACTACTACTTCAGCCGCGATTTTTCCGACATTGTCGGGGCCAACAAACAGTCCAAGGCCGAGGAGTGGGCCCAGGGTT  >  NC_002947/3589370‑3589666

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