Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,572	0	A	C	29.3%	75.3 / 12.0	41	V556G (GTG→GGG)	ybaL	putative transporter YbaL
Reads supporting (aligned to +/- strand):  ref base A (26/3);  new base C (0/12);  total (26/15)
Fisher's exact test for biased strand distribution p-value = 5.76e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCT  >  NC_000913/501433‑501710                                                                                                                                            |                                                                                                                                           aTCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTca                                                                                                                                            >  2:55414/1‑140 (MQ=255)  tcgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcac                                                                                                                                           <  1:1104162/140‑1 (MQ=255)  tcgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcac                                                                                                                                           >  1:1527993/1‑140 (MQ=255)                   aGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTg                                                                                                                          >  2:54839/1‑140 (MQ=255)                         ggCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGttt                                                                                                                    >  2:1351394/1‑140 (MQ=255)                          gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCATCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                                   >  1:84303/1‑140 (MQ=255)                            tttCCCGCCGCTTCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCGGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTCCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCca                                                                                                                 <  2:327208/139‑1 (MQ=255)                               cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                              >  1:139359/1‑140 (MQ=255)                                ccgccgcTTCTGGCATTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATTTTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCcagcag                                                                                                             <  2:145242/139‑1 (MQ=37)                                    cgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTcc                                                                                                         >  2:1135195/1‑140 (MQ=255)                                      cATCTGGCTTTTGCCTCTCGATGCTTGATGCGCCGCGGGCGCGTTTTATCATGCCTACTTTTTTTTCATTTTTTACACCCGGCACCCACCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAg                                                                                                       <  2:79089/140‑1 (MQ=255)                                        tCTGGCTTTTGCCTCTCGATGCCTGATGCGACGCGGGCGCTTTTTTTCATCCCTACATTTTTTTCATTTTTTCCACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                                                                     <  2:386853/140‑1 (MQ=255)                                                       tcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAttt                                                                                      >  1:883352/1‑140 (MQ=255)                                                            gCCTGAGGCGCCGCTGGCGCGTTTTATCATGCCTACATTTTTTTCATATTTTCCACCCGGCACCCACCGTTCCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                 <  2:223970/140‑1 (MQ=255)                                                             cttGATGCGACGCGGGCGTGTTTTATCATGCCTACTTTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                                <  2:97542/138‑1 (MQ=255)                                                               tGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCg                                                                              <  1:1351394/140‑1 (MQ=39)                                                                gATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCACCGTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                             <  1:231816/140‑1 (MQ=255)                                                                 aTGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                                            >  1:1475261/1‑140 (MQ=255)                                                                   gCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCa                                                                          >  1:232992/1‑140 (MQ=255)                                                                    cgccgcgGGCGCTTTTTATCATGCCTACATTTTTTTCATTTTTTCCACCCGCCACCCCCCTTTTCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt                                                                         <  2:156676/137‑1 (MQ=255)                                                                     gccGCTGGCGCTTCTTATCATGCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATc                                                                        <  2:1475261/138‑1 (MQ=255)                                                                     gccGCTGGCGCGTTTTATCATCCCTACATTTTTTTCATTTTTTACACCCGCCACCCACCTTTCCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATc                                                                        <  2:790278/138‑1 (MQ=255)                                                                       cGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                      >  1:55417/1‑140 (MQ=255)                                                                           ggCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc                                                                  >  1:1367469/1‑140 (MQ=255)                                                                           ggCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc                                                                  >  1:29050/1‑140 (MQ=255)                                                                            gcgcGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCt                                                                 >  1:949207/1‑140 (MQ=255)                                                                                 tCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                            >  1:756437/1‑140 (MQ=255)                                                                                 tCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                            >  1:84935/1‑140 (MQ=255)                                                                                  tttATCATCCCTACATTTTTTTCATTTTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCg                                                           <  2:183705/139‑1 (MQ=255)                                                                                    tATCATCCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGca                                                         <  2:949207/140‑1 (MQ=255)                                                                                       cATGCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacca                                                      <  2:1285219/140‑1 (MQ=255)                                                                                           ccTACTTTTTTTTAATTTTTTACATCCGGCAACCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                  <  1:534745/140‑1 (MQ=255)                                                                                           ccTACATTTTTTTCTTTTTTTCCGCCCGACACCCGCCTTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                  <  2:1367469/140‑1 (MQ=255)                                                                                           ccTACATTTTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                  <  1:1480866/140‑1 (MQ=255)                                                                                           ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                  >  2:613458/1‑140 (MQ=255)                                                                                              aCATTTTTTTCATTTTTTACACCCGCCACCCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                               <  2:21862/140‑1 (MQ=255)                                                                                                   tttttCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata                                          <  2:29050/140‑1 (MQ=255)                                                                                                    ttttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                         >  1:896365/1‑140 (MQ=255)                                                                                                     tttCATTTTTTACATCCGGCAACCACCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGAtata                                        <  1:1669031/140‑1 (MQ=255)                                                                                                      ttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                       >  1:21862/1‑140 (MQ=255)                                                                                                       tCATTTTTTACACCCGCCACCCCCCTTTCCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                                      <  2:304425/140‑1 (MQ=255)                                                                                                        aatatTTTACATCCGGCACCCCCCGTTTCCCCCTTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                     <  2:139359/139‑1 (MQ=255)                                                                                                         atatTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                    >  1:321798/1‑140 (MQ=255)                                                                                                             tttACACCCGGCAACCCCCGTTTCCCCCGCCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACtt                                <  2:1653293/140‑1 (MQ=255)                                                                                                             tttACACCCGCCACCCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACtt                                <  2:756437/140‑1 (MQ=255)                                                                                                               tACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                              >  1:860593/1‑140 (MQ=255)                                                                                                               tACACCCGCCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                              <  2:1527993/140‑1 (MQ=255)                                                                                                                 cacccGGCACCCCCCTTTTCCCCCCCCACCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATc                            <  1:1135195/137‑1 (MQ=255)                                                                                                                   cccGCCACCCCCCGTTTCCCCCGCCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt                          <  2:1521418/139‑1 (MQ=255)                                                                                                                      ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                       >  2:1198463/1‑140 (MQ=255)                                                                                                                      ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                       >  1:405029/1‑140 (MQ=255)                                                                                                                       gCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                      >  2:342115/1‑140 (MQ=255)                                                                                                                       gCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                      >  2:1324668/1‑140 (MQ=255)                                                                                                                                gTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCGCTTCATCGTCATAATGGGcgcg             <  1:1318793/140‑1 (MQ=255)                                                                                                                                 tttACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg            >  2:376169/1‑140 (MQ=255)                                                                                                                                   ccccccTTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGc          <  2:55417/138‑1 (MQ=255)                                                                                                                                     ccccGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaa        <  2:232992/140‑1 (MQ=255)                                                                                                                                       ccGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata      <  1:54839/140‑1 (MQ=255)                                                                                                                                         gTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat    >  2:139564/1‑140 (MQ=255)                                                                                                                                          tCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc   >  1:276913/1‑140 (MQ=255)                                                                                                                                           caccacCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtct  >  2:1411636/1‑140 (MQ=255)                                                                                                                                            accaccTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg             >  2:141482/1‑128 (MQ=255)                                                                                                                                            accaccTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg             <  1:141482/128‑1 (MQ=255)                                                                                                                                            |                                                                                                                                           ATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCT  >  NC_000913/501433‑501710

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