Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,613,433	A→G	6.8%	K97E (AAG→GAG)	PP_4084 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,613,433	0	A	G	6.8%	170.0 / 2.3	59	K97E (AAG→GAG)	PP_4084	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base G (2/2);  total (29/30)
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

GGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGC  >  NC_002947/4613288‑4613578                                                                                                                                                  |                                                                                                                                                  ggATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGc                                                                                                                                                >  2:345563/1‑149 (MQ=255)       ggCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCAc                                                                                                                                           <  8:196923/148‑1 (MQ=255)        gCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTa                                                                                                                                         >  1:169211/1‑149 (MQ=255)             ggaggaGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCa                                                                                                                                      <  3:72998/147‑1 (MQ=255)                ggagCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCtt                                                                                                                                  >  3:5086/1‑148 (MQ=255)                   gCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAgg                                                                                                                              <  1:345563/149‑1 (MQ=255)                   gCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCACGGTCGGCGACAAGCCCCACTACCAGCTTCAgg                                                                                                                              >  3:175740/1‑149 (MQ=255)                     aaCTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGTCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAggc                                                                                                                             <  3:196090/148‑1 (MQ=255)                     aaCTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAggcg                                                                                                                            <  1:87662/149‑1 (MQ=255)                             gTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGgcctgcct                                                                                                                    <  2:169211/149‑1 (MQ=255)                                ttACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGa                                                                                                                  >  1:35872/1‑148 (MQ=255)                                   cGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCAGGCCGCCTGaagct                                                                                                              <  7:50588/149‑5 (MQ=255)                                   cGGAAGAAGACAAAGGCGCTGACCATTACCCCCGCGAAAGAGACCCCCCTTGGAAGCAATGAGCATCACGTCGACAGACGGGAAAAGGCCATCCAGATCAAGGTCGGCGACAACCCCCACTACCAGCTTCAGGAAGGCCGCCTGaa                                                                                                                 >  8:163955/1‑145 (MQ=255)                                   aggAAGAAGTGAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCTCGTTGGAACCAATGAGCTTCCCGTCGTCAGTAGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCAGGCCGCCTGaa                                                                                                                 <  7:163955/145‑2 (MQ=255)                                       tgaagTCAAAGGGGCTGACCATTACGCCGGCTTTAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGTGAAATGGCCATCCTTTTCATTGTCGGCGACAAGCCCGACTACCAGCTTCTGGCGTCCTGCCTGATAcctgta                                                                                                          <  4:242529/148‑5 (MQ=255)                                           gTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCACGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGc                                                                                                      <  4:175740/149‑1 (MQ=255)                                              aaaGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCaa                                                                                                    >  8:293778/1‑148 (MQ=255)                                              aaaGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACAACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAAc                                                                                                   >  8:311672/1‑149 (MQ=255)                                                         ccATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACa                                                                                                        >  4:135528/1‑133 (MQ=255)                                                         ccATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACa                                                                                                        <  3:135528/133‑1 (MQ=255)                                                          cATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGc                                                                                       >  5:229055/1‑149 (MQ=255)                                                            ttACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTg                                                                                        >  4:348400/1‑146 (MQ=255)                                                            ttACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTg                                                                                        <  3:348400/146‑1 (MQ=255)                                                            ttACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCt                                                                                      >  1:287491/1‑148 (MQ=255)                                                               cGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAAc                                                                                  >  1:185446/1‑149 (MQ=255)                                                                gCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCAGATACCGGGACAGCAACGGATGCTGGTGCTCAAc                                                                                  >  2:280493/1‑148 (MQ=255)                                                                    gcgAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAg                                                                             >  5:362045/1‑149 (MQ=255)                                                                      gAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGc                                                                            >  1:339953/1‑148 (MQ=255)                                                                              cccGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTggg                                                                   <  6:229055/149‑1 (MQ=255)                                                                               ccGCCTTGGAAGCAATGAGCATCCCTTCGACAGTCGGGAAATGGCCATCCAGTTAATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTCGCTGATaac                                                                                                              >  6:281247/1‑103 (MQ=255)                                                                               ccGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTgg                                                                                          <  8:278231/125‑1 (MQ=255)                                                                               ccGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTgg                                                                                          >  7:278231/1‑125 (MQ=255)                                                                                    ttGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCttt                                                              >  3:304491/1‑148 (MQ=255)                                                                                      ggAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGa                                                            <  2:339953/148‑1 (MQ=255)                                                                                      ggAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGa                                                            <  7:311672/148‑1 (MQ=255)                                                                                      ggAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGa                                                            <  4:304491/148‑1 (MQ=255)                                                                                      ggAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCAGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGa                                                            <  1:280493/148‑1 (MQ=255)                                                                                             aTGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTgacgac                                                    <  7:293778/149‑1 (MQ=255)                                                                                                 gCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCACACTACCAGCTTCAGGCCGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGata                                                >  2:70601/1‑149 (MQ=255)                                                                                                    tCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCTCTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCa                                             <  2:287491/149‑1 (MQ=255)                                                                                                    tCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTgg                                                                                          >  1:29308/1‑104 (MQ=255)                                                                                                    tCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTgg                                                                                          <  2:29308/104‑1 (MQ=255)                                                                                                     cccGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACGAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCACCTACAGCAAGCCTGGGCCTTTGAGTGCCGACGCTATCAg                                            >  4:181187/1‑149 (MQ=255)                                                                                                     cccGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACGAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTCACGCCGATATCAg                                            >  4:175622/1‑149 (MQ=255)                                                                                                         tCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGata                                                <  2:227785/141‑1 (MQ=255)                                                                                                         tCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGata                                                >  1:227785/1‑141 (MQ=255)                                                                                                             cAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGgc                                    <  2:185446/149‑1 (MQ=255)                                                                                                              aGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGgcg                                   >  2:178215/1‑149 (MQ=255)                                                                                                                tCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTgg                                                                    >  1:221969/1‑114 (MQ=255)                                                                                                                tCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTgg                                                                    <  2:221969/114‑1 (MQ=255)                                                                                                                   ggAAATGGCCATCCAGTTCATGGTCGGCGACGAGCCCCCCTTCCAGCTTCAGGCGGCCTGCCTGATACCGGGCCAGCAACGGATGGTGGTGCTGAACTACAGCAAGCCTGGGCCGTTGAGTGACGACGATATCAGCCACTGgcgcgcgc                              <  3:181187/149‑1 (MQ=255)                                                                                                                   ggAAATGGCCATCCAGTTCATGGTCGGCGACGAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAATCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGgcgcgcgc                              <  3:175622/149‑1 (MQ=255)                                                                                                                          gCCATCCAGTTCATTGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAg                        <  3:10537/148‑1 (MQ=255)                                                                                                                          gCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAg                        <  1:70601/148‑1 (MQ=255)                                                                                                                                        ggTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGCCGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTg         >  1:313954/1‑149 (MQ=255)                                                                                                                                        ggTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTg         >  6:184808/1‑149 (MQ=255)                                                                                                                                         gTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGc        <  2:35872/149‑1 (MQ=255)                                                                                                                                              cGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCCATTTGCCTGAg   >  8:160570/1‑149 (MQ=255)                                                                                                                                               gACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGc  <  4:5086/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTCAGGCGGCCTGCCTGATACCGGGACAGCAACGGATGCTGGTGCTCAACTACAGCAAGCCTGGGCCTTTGAGTGACGACGATATCAGCCACTGGCGCGCGCTTAAAGACGAGCTCCGATTTGCCTGAGC  >  NC_002947/4613288‑4613578

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