Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	4,165,884	G→A	7.0%	G235G (GGC→GGT)	yybS ←	putative integral inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	4,165,884	0	G	A	7.0%	178.7 / 3.4	57	G235G (GGC→GGT)	yybS	putative integral inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (29/24);  new base A (2/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

CAGCAATCACAGCCGCAGCCTTTGGATATTGTTTCTTATGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCGGGGGAGAATCGCCTTAAGAGCGGCTTAGCGA  >  BS168_reference_genome/4165754‑4166020                                                                                                                                   |                                                                                                                                         cAGCAATCACAGCCGCAGCCTTTGGATATTGTTTCTTATGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTc                                                                                                                                  >  1:251127/1‑139 (MQ=255)      aTCACAGCCGCAGCCTTTGGATATTGTTTCTTATGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGACCTTTTTCTAAAg                                                                                                                             >  3:93866/1‑139 (MQ=255)      aTCACAGCCGCAGCCTTTGGATATTGTTTCTTATGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGACCTTTTTCTAAAg                                                                                                                             >  3:89681/1‑139 (MQ=255)            gCCGCAGCCTTTGGATATTGTTTCTTATGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCaaa                                                                                                                        >  6:24387/1‑138 (MQ=255)                cAGCCTTTGGATATTGTTTCTTATGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATg                                                                                                                    >  4:171164/1‑138 (MQ=255)                             ttGTTTCTTATGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCa                                                                                                                          <  4:200624/119‑1 (MQ=255)                             ttGTTTCTTATGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCa                                                                                                                          >  3:200624/1‑119 (MQ=255)                                tttCTTATGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATc                                                                                                   <  8:222604/139‑1 (MQ=255)                                     tATGGCAATAGAAAAAGATAATTGATAAACCTTGGATGAAAATCTGAAACCCGAGTATAAACTCGCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTTAAGGCAAAATTGAAAGAATGACAATAATCTAATa                                                                                              >  1:121848/1‑139 (MQ=255)                                     tATGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATa                                                                                              >  7:120120/1‑139 (MQ=255)                                      aTGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGAGAATAATCAAATa                                                                                              <  4:3604/138‑1 (MQ=255)                                      aTGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATa                                                                                              <  7:32810/138‑1 (MQ=255)                                      aTGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATa                                                                                              <  5:24387/138‑1 (MQ=255)                                      aTGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATa                                                                                              <  3:171164/138‑1 (MQ=255)                                        ggCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAg                                                                                           >  4:80838/1‑139 (MQ=255)                                         gCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAg                                                                                           <  8:80246/138‑1 (MQ=255)                                         gCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAg                                                                                           <  2:171079/138‑1 (MQ=255)                                          cAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGta                                                                                         <  3:279060/139‑1 (MQ=255)                                          cAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGta                                                                                         <  3:255365/139‑1 (MQ=255)                                             tAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAATCTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAACGGCAAAAATGAAAGCATGAAAAAAATCAAATACAGtacta                                                                                      >  8:232236/1‑139 (MQ=255)                                                         aaTGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCg                                                                          >  6:58506/1‑139 (MQ=255)                                                                aaCCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATa                                                                                              <  8:219143/112‑1 (MQ=255)                                                                aaCCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATa                                                                                              >  7:219143/1‑112 (MQ=255)                                                                aaCCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACtt                                                                   >  8:215525/1‑139 (MQ=255)                                                                aaCCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGCAAACtt                                                                   >  1:99926/1‑139 (MQ=255)                                                                   cTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCaa                                                                <  3:80838/139‑1 (MQ=255)                                                                   cTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCaa                                                                <  4:246258/139‑1 (MQ=255)                                                                   cTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCaa                                                                >  1:222750/1‑139 (MQ=255)                                                                   cTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGACCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCaa                                                                <  4:93866/139‑1 (MQ=255)                                                                   cTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGACCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCaa                                                                <  4:89681/139‑1 (MQ=255)                                                                   cTTCGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTCGCCAATATCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCaa                                                                <  2:121848/139‑1 (MQ=255)                                                                     tGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAAct                                                              >  3:255160/1‑139 (MQ=255)                                                                      ggATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAActc                                                             >  1:161989/1‑139 (MQ=255)                                                                      ggATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAAct                                                              <  6:105493/138‑1 (MQ=255)                                                                      ggATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAAc                                                               <  7:215525/137‑1 (MQ=255)                                                                            aaaTCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTaa                                                        >  4:239799/1‑138 (MQ=255)                                                                                   aaaCCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACTATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCttt                                                <  2:251127/139‑1 (MQ=255)                                                                                          aGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACaa                                                                                                        <  2:13682/76‑1 (MQ=255)                                                                                          aGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACaa                                                                                                        >  1:13682/1‑76 (MQ=255)                                                                                              taAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGaa                                      <  4:199172/138‑1 (MQ=255)                                                                                                  cTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCg                                 <  3:239799/139‑1 (MQ=255)                                                                                                     cctaTTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCgggg                              >  2:186729/4‑139 (MQ=255)                                                                                                     cccATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTTAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTACAACTCTCTAAAAGGCTTTAGATTAGGAATATCgggg                              >  8:235697/1‑139 (MQ=255)                                                                                                     cccATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAg                                                                                           <  5:210932/78‑1 (MQ=255)                                                                                                     cccATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAg                                                                                           >  6:210932/1‑78 (MQ=255)                                                                                                     cccATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCgggg                              >  4:10932/1‑139 (MQ=255)                                                                                                     cccATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCgggg                              >  7:15561/1‑139 (MQ=255)                                                                                                     cccATTTAAAGCGATAGAGTACAGCATTTGGCCTTTATCAAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGTAATCATCAACTCTCTAAAAGGCTTTAGATAAGGAATATCtggc                              >  4:52567/1‑138 (MQ=255)                                                                                                         tttAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCGGGGgaga                          <  2:222750/139‑1 (MQ=255)                                                                                                          ttAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCGGGGGAGaa                         >  1:123958/1‑139 (MQ=255)                                                                                                          ttAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCGGGGGAGaa                         >  5:78246/1‑139 (MQ=255)                                                                                                            aaaGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCGGGGGAGAAt                        <  3:10932/138‑1 (MQ=255)                                                                                                                   tAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCGGGGGAGAATCGCCtt                  >  2:134113/1‑137 (MQ=255)                                                                                                                      atgaCAGCAGTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATCATCAAATACAGTACTACAACAAGTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCGGGGGAGAATCGCCTTAagag             <  3:212130/136‑1 (MQ=255)                                                                                                                             cATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCGGGGGAGAATCGCCTTAAGAGCGGCt        >  4:170127/1‑137 (MQ=255)                                                                                                                             cATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAAT‑ATCAAATACAGTACTACAACACTTTCCGGAAACTTCAACTCTCTTAAAGGCTTTAGATAAGGAATATCGGGGGAGAATCGCCTTAAGAGCGGCTTAg     >  6:73088/1‑139 (MQ=255)                                                                                                                                ttGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCGGGGGAGAATCGCCTTAAGAGCGGCTTAGc    >  3:254446/1‑138 (MQ=255)                                                                                                                                  ggCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCg                                                                          >  2:256067/1‑66 (MQ=255)                                                                                                                                  ggCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCg                                                                          <  1:256067/66‑1 (MQ=255)                                                                                                                                  ggCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTATATTAGGAATATCGGGGGAGAATCGCCTTAagag             <  5:300932/127‑1 (MQ=255)                                                                                                                                  ggCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCGGGGGAGAATCGCCTTAagag             >  6:300932/1‑127 (MQ=255)                                                                                                                                  ggCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCGGGGGAGAATCGCCTTAAGAGCGGCTTAGCGa  >  6:258752/1‑138 (MQ=255)                                                                                                                                  ggCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCGGGGGAGAATCGCCTTAAGAGCGGCTTAGCGa  <  2:161989/138‑1 (MQ=255)                                                                                                                                  ggCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCGGGGGAGAATCGCCTTAAGAGCGGCTTAGCGa  <  1:134113/138‑1 (MQ=255)                                                                                                                                  ggCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGAGGAAACTTCAACTCTCTAATAGGCTTTAGAGTTGGAATATCGGGGGAGAATCGCCTTAAGAGCGGCTTAGCGa  >  1:30529/1‑138 (MQ=255)                                                                                                                                  ggCCTTATTCTAAAGTCAAAAATGAAAGCATGTCAATACTCAAATACAGTTCTACAACACTTTGCGGAAACTTCAACTCTCTAACAGGCTTTAGATTACGAATATCGGGGGAGAATCGCCTTAAGAGCGGCTTAGCGa  <  7:232236/138‑1 (MQ=255)                                                                                                                                   |                                                                                                                                         CAGCAATCACAGCCGCAGCCTTTGGATATTGTTTCTTATGGCAATAGAAAAAGATAAATGATAAACCTTGGATGAAAATCAGAAACCCGAGTATAAACTCCCCATTTAAAGCGATAGAGTACAGCATTTGGCCTTTTTCTAAAGGCAAAAATGAAAGCATGACAATAATCAAATACAGTACTACAACACTTTGCGGAAACTTCAACTCTCTAAAAGGCTTTAGATTAGGAATATCGGGGGAGAATCGCCTTAAGAGCGGCTTAGCGA  >  BS168_reference_genome/4165754‑4166020

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