Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,575	0	A	C	28.9%	75.7 / 11.6	38	V555G (GTG→GGG)	ybaL	putative transporter YbaL
Reads supporting (aligned to +/- strand):  ref base A (26/1);  new base C (0/11);  total (26/12)
Fisher's exact test for biased strand distribution p-value = 9.97e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTC  >  NC_000913/501436‑501711                                                                                                                                             |                                                                                                                                         gtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTcacca                                                                                                                                          >  2:1395750/1‑140 (MQ=255)    ggTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCt                                                                                                                                       >  1:1019098/1‑140 (MQ=255)                  aTAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTcacca                                                                                                                                          <  2:1429984/123‑1 (MQ=255)                  aTAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTcacca                                                                                                                                          >  1:1429984/1‑123 (MQ=255)                  aTAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGc                                                                                                                         >  1:4652/1‑140 (MQ=255)                  aTAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGc                                                                                                                         >  1:303581/1‑140 (MQ=255)                       gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                                    >  2:921277/1‑140 (MQ=255)                          ttgccgccgcATCGGGCATTTCGCTCTCGGTGCGGGGGGGGCCGCGGGCGCGTCTTACAAGCCCTTCGTTTTTTTCAT‑TTTTAACTCCATCAAACCACTGTTCCCCCGGCCCCCCCCCCACCCGCCGGTGGCGTTTCcag                                                                                                                 <  2:1261011/136‑1 (MQ=255)                          ttCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCACCGGTGGCGTTTCcag                                                                                                                 <  2:1063388/140‑1 (MQ=255)                          ttCACGCCGCACCTGGCATTTGCCTCTCGATGCCTGATGCGCCGCTGGCGCTTCTTATCATGCCTACATTTTTTTCATATTTTACACCCGGC‑ACCCCCCGTTTACCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCcag                                                                                                                 <  2:4608/140‑1 (MQ=255)                          tccccgccgcATCTGGCTTTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCTTTTTATCATGCCTACATATTTTTCATTTTTTACACCCGGC‑AACCACCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCcag                                                                                                                 <  2:303581/136‑1 (MQ=255)                           ccACCCCGCATCTGGCATTTGCCCCTCGTGGCTTGAGGCAGCGCGGGCGCGTTTTACCTTCCCAACTTATTTTTTTTATTTTACACCCGCC‑CACCACCTTTCCCCCCGTCCCCACCCCACCCGCCGGTGGCGTTTCcagc                                                                                                                <  2:1019098/139‑1 (MQ=255)                            cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                               >  2:1287038/1‑140 (MQ=255)                            cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                               >  1:1022764/1‑140 (MQ=255)                            cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                               >  1:1159866/1‑140 (MQ=255)                            cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                               >  1:1246262/1‑140 (MQ=255)                                   cATCTGGCTTTTGCCTCTCGATGCCTGAGGCGACGCTGGCGCGTTTTATCATGCCTACATTTTTTTAATTTTTTCCACCCGGC‑ACCCCCCGTTTCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAg                                                                                                        <  2:705911/140‑1 (MQ=255)                                   cATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAg                                                                                                        >  1:231304/1‑140 (MQ=255)                                   cATCTGGATTTTGACTCTCGATGCCTGAGGCGGCCGTGGCGCGTTTTATCATCCCTACTTTTTTTTCATTTTTAACACCCGAC‑ACCCCCGGTTCCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAg                                                                                                        <  2:419311/140‑1 (MQ=255)                                     tCTGGCATTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGC‑ACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                                                                      <  2:1309378/140‑1 (MQ=255)                                                      gATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCa                                                                                     >  1:1053245/1‑140 (MQ=255)                                                         gCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACACCCGGC‑ACCCCCCGTTCCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                  <  2:443196/140‑1 (MQ=38)                                                         gCCTGAGGCGACGCGGGCGCTTCTTATCATGCCTACATTTTTTTCATTTTTTACACCCGGC‑ACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                  <  2:564331/140‑1 (MQ=255)                                                          ccTGATGCGACGCTGGCGCTTCTTATCATGCCTACATTTTTTTCATTTTTTACATCCGGC‑AACCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                                 <  1:517563/140‑1 (MQ=255)                                                             ggggggACGCTGGGGCGTCTTACCATCCCAACATATTTTTAATATTTGACATCCGCC‑ACCCACCGTTTCCCCCCTCCCCACCGCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                              <  2:1022764/135‑1 (MQ=255)                                                             gATGCGCGGCGGGCGCGTTTTATCATGCCTACATTTTTTTCATATTTTACACCCGGC‑AACCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                              <  2:1159866/140‑1 (MQ=255)                                                             gATGCAACGCTGGCGCTTCTTATCATGCCTACATTTTTTTCATTTTTTACATCCGGC‑ACCCACCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                              <  2:915624/140‑1 (MQ=255)                                                              aTGCGACGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACACCCGGC‑ACCCACCGTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                                             <  1:921277/140‑1 (MQ=255)                                                                   aCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCa                                                                        >  2:412921/1‑140 (MQ=255)                                                                    cGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                       >  2:3841/1‑140 (MQ=255)                                                                    cGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                       >  2:740130/1‑140 (MQ=255)                                                                                 tATCATGCCTACATATTTTTCATTTTTTACACCCGCC‑ACCCACCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGca                                                          <  2:517484/140‑1 (MQ=255)                                                                                  ttCATCCCTACATTTTTTTAATTTTTTCCACCCGGC‑ACCCCCCTTTTCCCCCGCCCCCCCCTCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcac                                                         <  2:1286640/139‑1 (MQ=255)                                                                                   tCATGCCTACATTTTTTTCATATTTTACACCCGGC‑ACCCCCCGTTTCCCCCGCCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc                                                        <  2:231304/140‑1 (MQ=255)                                                                                      tGCCTACATTTTTTTCATATTTTACATCCGGC‑ACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg                                                     <  1:1286598/140‑1 (MQ=255)                                                                                        ccTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                   >  1:352891/1‑140 (MQ=255)                                                                                         cTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTc                                                  >  1:517484/1‑140 (MQ=255)                                                                                          tACATTTTTTTCATTTTTTACACCCGGC‑ACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                                 <  1:1097296/140‑1 (MQ=255)                                                                                          tACATATTTTTCATATTTTACACCCGGC‑ACCCACCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                                 <  1:56010/140‑1 (MQ=255)                                                                                           aCATTTTTTTAATTTTTCACCCCCGCC‑ACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                                <  2:395730/140‑1 (MQ=255)                                                                                            cATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                               >  1:1286640/1‑140 (MQ=255)                                                                                               ttttttCATTTTTTACACCCGCC‑ACCCCCCTTTTACCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                            <  2:1265861/139‑1 (MQ=255)                                                                                                tttttAATATTTTACACCCGGC‑ACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata                                           <  2:4652/140‑1 (MQ=255)                                                                                                  tttCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGAtata                                         >  1:395730/1‑140 (MQ=255)                                                                                                   ttCATATTTTACATCCGGC‑AACCACCGTTTATCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                        >  1:1265861/1‑140 (MQ=255)                                                                                                   ttCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                        >  2:885228/1‑140 (MQ=255)                                                                                                           ttACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc                                >  1:1314610/1‑140 (MQ=255)                                                                                                            tccccccGGC‑ACCCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                               <  2:542188/135‑1 (MQ=255)                                                                                                            tACATCCGGC‑ACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                               <  1:1395750/140‑1 (MQ=255)                                                                                                            tACACCCGCC‑ACCCCCCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                               <  2:1246262/140‑1 (MQ=255)                                                                                                                   ggC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                        >  2:885730/1‑140 (MQ=255)                                                                                                                              gTTTACCACGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg              >  1:744809/1‑140 (MQ=255)                                                                                                                               tttACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg             >  1:131175/1‑140 (MQ=255)                                                                                                                               ttccccccGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg             <  2:576753/136‑1 (MQ=255)                                                                                                                                ttcccccGCCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCggg            <  2:1314610/137‑1 (MQ=255)                                                                                                                                tccccccGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCggg            <  2:915275/137‑1 (MQ=255)                                                                                                                                  cccccGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCATGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCa          <  1:1236123/139‑1 (MQ=255)                                                                                                                                    cccGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaat        <  1:718138/140‑1 (MQ=255)                                                                                                                                       gTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat     >  1:354674/1‑140 (MQ=255)                                                                                                                                       gTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat     >  1:1048959/1‑140 (MQ=255)                                                                                                                                        tCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc    >  1:232209/1‑140 (MQ=255)                                                                                                                                          accaccTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtctc  <  2:1048959/140‑1 (MQ=255)                                                                                                                                             |                                                                                                                                         GTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGC‑AACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTC  >  NC_000913/501436‑501711

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