Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,207	0	T	C	24.5%	107.1 / 12.4	53	S215G (AGT→GGT) ‡	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (34/6);  new base C (0/13);  total (34/19)
Fisher's exact test for biased strand distribution p-value = 3.22e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556070‑3556345                                                                                                                                          |                                                                                                                                           cTCGGCACTGACGCGCTTTTCACTGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTcc                                                                                                                                          >  1:1250319/1‑140 (MQ=255)               gcTTTTCACCGATGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACCCCCCGCCAGTGTAGCg                                                                                                                            <  2:548352/140‑1 (MQ=255)                         gacgaaAAAAAAGCGTTGGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACCCCCCGCCAGTGTAGCGATTTCACGCt                                                                                                                  <  2:1041154/140‑1 (MQ=255)                             acAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                              >  1:1347977/1‑140 (MQ=255)                                    aGCGTTGGGTAATTTTTTCACTTTCGACTTCAGGCGAGCGGTTTTTCCCCGCCCCCTGGCCCCGGGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                       <  2:1671913/140‑1 (MQ=255)                                    aGCGTTCGGTGTTTTTTTCCTTTCCGACTTAAAGCGAGCCGGTTTTCCCCGCCCCCTGGTCGCGCGGCAGGTTATAAATTTTCTGTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                       <  2:1185352/140‑1 (MQ=255)                                        ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                   >  1:1775730/1‑140 (MQ=255)                                              tattttttCCCTTTCGACTTAAGGTGAGCCGGTATTCCCCGCCCCCGGGCCGGGCGGAAGTTTAAGAATTTTTTGTCAGCGTCGGAAAAGTATCACCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac                                                                                             <  2:1250319/136‑1 (MQ=255)                                                ttttttCCCTTTCGACTTCAAGGGAGACGGTTTTCCCCGGCCCCTGGTCGCGGGGCAGGTTATAAATTTTCTGTTCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                           <  2:1817258/138‑1 (MQ=255)                                                tatTTTCGTTTACGACTCAAGGCGAGACGGTTTTCCCCGCCCCCGGGTCGCGCGGCAGGTTAAAAATATTCTTTCAATGAGGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                           <  2:743168/140‑1 (MQ=255)                                                   tttCACTTTCGACTTCAAGCGAGACGGTATTA‑CCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTa                                                                                       >  2:1555485/1‑140 (MQ=255)                                                     tCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                      >  1:1240444/1‑140 (MQ=255)                                                     tCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                      >  1:1952742/1‑140 (MQ=255)                                                           tCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                >  1:795092/1‑140 (MQ=255)                                                                  caggCGAGACGGTTTTCCCCGGCCCGTGGTCGCGCGGCGGGTTAGAAATTTTCTGTTTAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                         <  2:1251296/137‑1 (MQ=255)                                                                  aaaGCGAGACGGTTTTCCCCGGCCCGGGGTCGCGCGGCGGGTTATAAATTTTCTTTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                         <  2:1095712/139‑1 (MQ=255)                                                                    aagagggcGGTATCCCCCGCCCCCGGGTCGCGCGGCGGTTTATAAATTTTCTTTTCAGGCAGGAAAACCCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                       <  2:487500/133‑1 (MQ=255)                                                                    aGCGAGCCGTTTTTCCCCGCCCCCTGGCCGCGCGGCAGTTTATGAATATTCTGTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                       <  2:29814/140‑1 (MQ=255)                                                                     gCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                                      >  1:307835/1‑140 (MQ=255)                                                                      ggagaCGTTATTCCCCGCCCCTTGGCCGCGGGGCAGGTTATGAATTTTTTGTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                     <  2:307835/139‑1 (MQ=255)                                                                         gaCGGTATTACCCGGCCCCTCGTGGCGCGGCAGGTTGTGAATATTCTGTTCATGTAGTGCAAAACTCCCCGCCAGTGTAGCGCGTTCACGCTCAGCACCGTGGCGCGGCCCACCAGCTAATAGAACTTCTCCACTCATCt                                                                  >  2:1170215/1‑140 (MQ=255)                                                                          gcGGTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTAGGAATATTCTGTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                 <  1:11597/139‑1 (MQ=255)                                                                            gttATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTTTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca                                                               <  2:1421469/138‑1 (MQ=255)                                                                               tttccccGCCCCCTGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACaa                                                            <  2:43392/136‑1 (MQ=255)                                                                                      ggCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACaactgtctc                                                     >  1:487500/1‑133 (MQ=255)                                                                                      ggCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                     >  1:2028211/1‑140 (MQ=255)                                                                                        ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  1:882650/1‑140 (MQ=255)                                                                                        ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  1:12843/1‑140 (MQ=255)                                                                                        ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  2:974081/1‑140 (MQ=255)                                                                                            tGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                               >  1:1229775/1‑140 (MQ=255)                                                                                             ggccgcgggGCAGTTTATAAATTTTCTGTCAATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCg                                              <  2:1347977/137‑1 (MQ=255)                                                                                             ggTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCg                                              <  2:1243099/140‑1 (MQ=255)                                                                                                  cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg                                         >  1:861454/1‑140 (MQ=255)                                                                                                  cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg                                         >  2:590654/1‑140 (MQ=255)                                                                                                  cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg                                         >  2:83785/1‑140 (MQ=255)                                                                                                    cgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                       >  2:1489872/1‑140 (MQ=255)                                                                                                       gggtttATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                    <  2:1229775/136‑1 (MQ=255)                                                                                                       gggtttATAAATATTCTGTTCATGAAGAAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                    <  2:703839/136‑1 (MQ=255)                                                                                                       cAGGTTAGAAATCTTGTGTCCAGGCCGGGAAACCCGCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAGGGCGCGGCACACCAGCTAATAGAACTTCGCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                    <  1:93498/140‑1 (MQ=255)                                                                                                        ggtttATGAATTTTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                   <  2:12843/137‑1 (MQ=255)                                                                                                                   atTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg                        >  2:1111540/1‑140 (MQ=255)                                                                                                                    ttCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAcct                       >  1:703839/1‑138 (MQ=255)                                                                                                                    ttCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       >  1:1938793/1‑140 (MQ=255)                                                                                                                    ttCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       >  2:1030716/1‑140 (MQ=255)                                                                                                                    ttCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       >  1:1192936/1‑140 (MQ=255)                                                                                                                     tCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                      >  1:1301290/1‑140 (MQ=255)                                                                                                                     tCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                      >  2:531462/1‑140 (MQ=255)                                                                                                                       tGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                    >  2:1890902/1‑140 (MQ=255)                                                                                                                         ttCATGCAGGGAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  <  1:1338172/140‑1 (MQ=255)                                                                                                                         ttCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  1:1405789/1‑140 (MQ=255)                                                                                                                         ttCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  2:1338172/1‑140 (MQ=255)                                                                                                                           ttgcgaGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTCCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg                <  2:882650/135‑1 (MQ=255)                                                                                                                            aTGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg               >  1:959937/1‑140 (MQ=255)                                                                                                                            aTGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg               >  2:779790/1‑140 (MQ=255)                                                                                                                              gCAGGGAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc             <  1:779790/140‑1 (MQ=255)                                                                                                                                 gggAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag          >  1:276264/1‑140 (MQ=255)                                                                                                                                   gAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  2:981873/1‑130 (MQ=255)                                                                                                                                   gAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  <  1:981873/130‑1 (MQ=255)                                                                                                                                       aaCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc    <  2:1405789/140‑1 (MQ=255)                                                                                                                                       aaCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc    >  1:1173897/1‑140 (MQ=255)                                                                                                                                         cTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt  >  1:987090/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                           CTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556070‑3556345

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