Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	2,459,136	A→C	14.9%	K23N (AAA→AAC)	U712_12535 →	ComE operon protein 4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,459,136	0	A	C	14.9%	157.6 / 2.2	66	K23N (AAA→AAC)	U712_12535	ComE operon protein 4
Reads supporting (aligned to +/- strand):  ref base A (43/12);  new base C (4/7);  total (47/19)
Fisher's exact test for biased strand distribution p-value = 9.61e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCGTCATATTTTATTTTCATAACACATATGGTTTAGAGGATAATAGCTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAGCTCGTTTCTGAAAACGAGATGATTTTTAT  >  CP006881/2459000‑2459270                                                                                                                                         |                                                                                                                                       cTTTCGTCATATTTTATTTTCATAACACAAATGGTTTAGAGGATAATAGCTCAAAGGAGGGGCAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATAATTTATTGAATCaaaacc                                                                                                                                      >  3:264604/1‑137 (MQ=255)                     aaaaaCACATATGGTTTAGAGTATAATACCCCAACGGAGGGTAAACCATTGAAGATAGGCTTTATCGGCACAGGAATTATGGTAACTATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGa                                                                                                                  <  1:84047/136‑1 (MQ=255)                               atGGTTTAGAGGCTAATAGCTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCAGTTATTGAATCATACGCGGCCGATCCCTCAACCATGACAATCACaaa                                                                                                        >  8:47679/1‑139 (MQ=255)                               atGGTTTAGAGGATAATAGCTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATATATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACaaa                                                                                                        >  1:337188/1‑139 (MQ=255)                                     tAGAGGATAATAGCTCAAAGGAGGGGAACCCATTGATGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAACGCGGCCTCTCCCTCAAACATGACAATCACAAACCGca                                                                                                   >  4:310610/1‑138 (MQ=255)                                     tAGAGGATAATAGCTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGca                                                                                                   >  4:29271/1‑138 (MQ=255)                                            taataGCTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAACGCGGCCGCTCCCTCAAACATGACAATCACAAACCGCACAATTGa                                                                                            >  6:315802/1‑138 (MQ=255)                                             aataGCTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGATTCATTTATTGAATCAAAAGCGGCCGCTCCCTCAAACATGACTATCACAAACCGCACAATTGaa                                                                                           >  2:291543/1‑138 (MQ=255)                                             aataGCTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCTATCCCTCAAACATGACAATCACAAACCGCACAATTGaa                                                                                           >  6:367967/1‑138 (MQ=255)                                             aataGCTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAACCATGACAATCACAAACCGCACAATTGaa                                                                                           >  1:338377/1‑138 (MQ=255)                                             aataGCTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATcaca                                                                                                          <  6:185026/123‑1 (MQ=255)                                             aataGCTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATcaca                                                                                                          >  5:185026/1‑123 (MQ=255)                                               taGCTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCTAATCATTTATTGAATCAAAAGCTGCCGATCCCTCAAACATGACAATCACAACCCGCACAATAGaaa                                                                                          >  3:267908/1‑137 (MQ=255)                                                  cTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGTCCGATCCCTCAAACATGACAATCACAAACCGCACAAATGAAAAAGCg                                                                                     >  5:29294/1‑139 (MQ=255)                                                  cTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCg                                                                                     >  5:104865/1‑139 (MQ=255)                                                  cTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATAATTTATTGAATCAAAAGCGTCCGATCCCTCAAACATGACAATCACAAGCAGCACAATTGAAAAAGCg                                                                                     >  7:167094/1‑139 (MQ=255)                                                    cAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGtt                                                                                   >  5:173181/1‑139 (MQ=255)                                                       aggaggGGAAACCAATGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGACGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAATCGTTGc                                                                                 >  1:255392/1‑138 (MQ=255)                                                       aggaggGCAAAACATTGAAGATCGGCTTAATCGGCATTCTAAATAGGGGTACGATCCTTATCCAATCATTTATTGAATCAAACGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGca                                                                                <  3:186831/139‑1 (MQ=255)                                                            gggAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATa                                                                            >  3:111345/1‑138 (MQ=255)                                                            gggAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGACGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATa                                                                            >  5:2543/1‑138 (MQ=255)                                                               aacccATTGACTATTGCCTTTATCGCCACTGTAATAATGGGTACGAGCCTAATCGAAACATTTATTGAATCAAACGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATaaaaa                                                                        <  3:146075/136‑1 (MQ=255)                                                               aaaCCATTGAATATCGGCTTTATCGGCACAGTAAATATGGGTACGATCCTTATCGAATCATTTATTCAAAAAAACGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATaaaaa                                                                        <  3:88176/139‑1 (MQ=255)                                                                 aCCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAACGCGGCCGATCCCTCAGACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAAc                                                                      >  5:271479/1‑139 (MQ=255)                                                                        aaGAAAGGCTAAATCGGCACAGGAATAATTGGCACGTCCCTAATCCACTCTTTTATTTACAAAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTac                                                                <  6:373329/138‑1 (MQ=255)                                                                         agaacggCTTTAACGCCACCGTAAAAATGGGTACCATCCTTAATCAAAAATTTATTCACGAAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTacaa                                                              <  8:248790/134‑1 (MQ=255)                                                                         acatcggCTTTAGCGCCACAGTAAATTTGGGTACGATCCTTATCGAATAATTTATTTAAAAAAACGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTacaa                                                              <  3:238049/134‑1 (MQ=255)                                                                         aGATAGGCTTTATCGGCACAGTAATTATGGGAACGATCCTTATCGAATCATTTATTTAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTacaa                                                              <  5:350719/139‑1 (MQ=255)                                                                             aGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGCTCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAg                                                           >  5:287/1‑138 (MQ=255)                                                                             aGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAg                                                           >  3:147919/1‑138 (MQ=255)                                                                             aGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAg                                                           >  2:331999/1‑138 (MQ=255)                                                                             aGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAg                                                           >  5:329440/1‑138 (MQ=255)                                                                             aGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAg                                                           >  1:369402/1‑138 (MQ=255)                                                                             aGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAg                                                           >  3:202221/1‑138 (MQ=255)                                                                               gCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCa                                                         >  3:398599/1‑138 (MQ=255)                                                                                  ttATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATaaa                                                     >  2:284322/1‑139 (MQ=255)                                                                                    aTCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACg                                                   >  7:197806/1‑139 (MQ=255)                                                                                    aTCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACg                                                   >  1:398541/1‑139 (MQ=255)                                                                                    aTCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACg                                                   >  2:368444/1‑139 (MQ=255)                                                                                         cacaGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTAcc                                              >  4:51512/1‑139 (MQ=255)                                                                                         cacaGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTAcc                                              >  6:98380/1‑139 (MQ=255)                                                                                           acgcaaaGATGGCAACGAACCTAATCTAATCATTTATTCAAACAACAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGa                                            <  3:3835/135‑1 (MQ=255)                                                                                                      ggTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACCTGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAACCGTTACCGAAAGTCTCGaa                                  >  3:160717/1‑138 (MQ=255)                                                                                                      ggTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGaa                                  >  7:74105/1‑138 (MQ=255)                                                                                                      ggTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGaa                                  >  4:22640/1‑138 (MQ=255)                                                                                                        tACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAg                               >  4:111115/1‑139 (MQ=255)                                                                                                        tACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAg                               >  1:183591/1‑139 (MQ=255)                                                                                                        tACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGACTCGAAAAg                               >  2:161543/1‑139 (MQ=255)                                                                                                                  aTCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTacaa                                                              >  2:300756/1‑98 (MQ=255)                                                                                                                  aTCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTacaa                                                              <  1:300756/98‑1 (MQ=255)                                                                                                                  aTCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACAGAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAGCTCGTTTCt                      >  6:134768/1‑138 (MQ=255)                                                                                                                  aTCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTAAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCTCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAGATCGTTTCt                      >  7:133086/1‑138 (MQ=255)                                                                                                                      aaTCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAGCTCGTTTCTGaaa                  >  6:260425/1‑138 (MQ=255)                                                                                                                       aTCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAGCTCGTTTCTGAAAAc                >  4:200780/1‑139 (MQ=255)                                                                                                                        tCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAGCTCGTTTCTGAAAAc                >  7:64690/1‑138 (MQ=255)                                                                                                                           tttATTGAATCAAACGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAGCTCGTTTCTGAAAACgaga            <  5:315802/139‑1 (MQ=255)                                                                                                                               ttGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAGCTCGTTTCTGAAAACGAgatgat        <  3:200780/139‑1 (MQ=255)                                                                                                                                    tCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAACGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCTAAAGTCTCGAAAAGCTCGTTACTGAAAACGAGATGAttatta   >  2:134822/1‑139 (MQ=255)                                                                                                                                    tCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATaaaaaa                                                                       <  7:247702/71‑1 (MQ=255)                                                                                                                                    tCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATaaaaaa                                                                       >  8:247702/1‑71 (MQ=255)                                                                                                                                    tCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCAt                                                        <  8:330069/86‑1 (MQ=255)                                                                                                                                    tCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCAt                                                        >  7:330069/1‑86 (MQ=255)                                                                                                                                    tCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAGCTCGTTTCTGAAAACGAGATGATTTTTa   >  8:80070/1‑139 (MQ=255)                                                                                                                                     cAACAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAGCTCGTTTCTGAAAACGAGATGATTTTTAt  <  1:368444/139‑1 (MQ=255)                                                                                                                                     cAAACGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAGCTCGTTTCTGAAAACGAGATGATTTTTAt  <  1:331999/139‑1 (MQ=255)                                                                                                                                     cAAACGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAGCTCGTTTCTGAAAACGAGATGATTTTTAt  <  6:329440/139‑1 (MQ=255)                                                                                                                                     cAAAAGCGGCCTATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAGCTCGTTTCTGAAAACGAGATGATTTTTAt  <  5:367967/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       CTTTCGTCATATTTTATTTTCATAACACATATGGTTTAGAGGATAATAGCTCAAAGGAGGGGAAACCATTGAAGATAGGCTTTATCGGCACAGGAAATATGGGTACGATCCTTATCGAATCATTTATTGAATCAAAAGCGGCCGATCCCTCAAACATGACAATCACAAACCGCACAATTGAAAAAGCGTTGCACATAAAAAACCGCTACAACAGCATAAACGTTACCGAAAGTCTCGAAAAGCTCGTTTCTGAAAACGAGATGATTTTTAT  >  CP006881/2459000‑2459270

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