Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,189,448	C→T	Q1367* (CAG→TAG)	rpoC →	RNA polymerase, beta prime subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,189,448	0	C	T	100.0%	179.2 / NA	50	Q1367* (CAG→TAG)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/28);  total (22/28)

GTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTA  >  NC_000913/4189300‑4189595                                                                                                                                                      |                                                                                                                                                    atCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGCAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGGTATCGTGGGGCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAg                                                                                                                                                   >  2:898545/2‑151 (MQ=255)   ccTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGa                                                                                                                                                 <  2:671131/151‑1 (MQ=255)   ccTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGa                                                                                                                                                 <  1:1150283/151‑1 (MQ=255)   ccTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGa                                                                                                                                                 <  1:651387/151‑1 (MQ=255)    cctCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGAt                                                                                                                                                <  2:924784/149‑1 (MQ=255)    cTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCCGGCAGGTACCGGTTACGCGTACCACTAGGa                                                                                                                                                 <  2:951696/151‑1 (MQ=255)    cTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGAt                                                                                                                                                <  1:1017316/151‑1 (MQ=255)    cTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGAt                                                                                                                                                >  1:933744/1‑151 (MQ=255)      tCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCg                                                                                                                                              <  2:405022/151‑1 (MQ=255)         tctcCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTAt                                                                                                                                           <  2:158210/151‑1 (MQ=255)             cgcgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGt                                                                                                                                       <  1:914325/151‑1 (MQ=255)                 tcATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCTTCGc                                                                                                                                    >  2:689689/2‑150 (MQ=255)                     cGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGCTCGTATGCGTCGCCGtg                                                                                                                                >  2:470/1‑150 (MQ=255)                     cGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGtgc                                                                                                                               <  1:160420/151‑1 (MQ=255)                     cGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGtg                                                                                                                                >  2:914473/1‑150 (MQ=255)                     cGTCCCAGGAGACAACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGCAAGAGAACGTTATCGTGCGTCGTCTGAT‑CCCGGCCGGGACCGGTTACGTGTACCACTAGGATCGTATGCGTCGCCCtg                                                                                                                                >  2:883200/1‑150 (MQ=255)                      gTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGtg                                                                                                                                >  1:604352/1‑149 (MQ=255)                       ttCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAgg                                                                                                                                                  <  1:149539/130‑1 (MQ=255)                       ttCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAgg                                                                                                                                                  >  2:149500/1‑130 (MQ=255)                         ccAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCg                                                                                                                           >  1:427997/1‑151 (MQ=255)                                    cTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTc                                                                                                                 >  2:298564/1‑150 (MQ=255)                                        cgTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCt                                                                                                            <  2:882849/151‑1 (MQ=255)                                              gACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCAcc                                                                                                       >  2:216328/1‑150 (MQ=255)                                                     gcagcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGTGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCCGGTGa                                                                                               >  2:624404/1‑151 (MQ=255)                                                      cggcCGTTGCGGGCAATATCGACGAACTGCGCGGCCTGAAAAAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGGTACGCGTACCCCTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCAcgtcagaagac                                                                                              <  2:1183176/149‑11 (MQ=255)                                                           gTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCCTACCACTAGGATCGTATGCGTCCCCGTGCTGCGGGGGAAGCTCCGGCTGCACCGCAGGTGACTGCag                                                                                         >  1:797692/1‑151 (MQ=255)                                                               cGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCagaag                                                                                      <  2:250123/150‑1 (MQ=255)                                                               cGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCagaag                                                                                      <  1:748260/150‑1 (MQ=255)                                                               cGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGAGACTGCagaag                                                                                      <  2:810995/150‑1 (MQ=255)                                                                 ggCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGAc                                                                                    <  2:933350/150‑1 (MQ=255)                                                                       cgcgACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCAt                                                                                >  1:1076062/1‑148 (MQ=255)                                                                         cgacgaACTGCGCGGCCTGAAAGGGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGGACCGGAGGTGACTGCAGAAGACGCATCTGcc                                                                           <  2:519684/151‑1 (MQ=255)                                                                                                  aaCGTTATCGTGGGTCGTCTGAT‑GCCGGCAGTTACCGGTTACGCTTACTCCTAGGATCTCATGTGTCGCCGTGTGGCTGGTGAAGCTCCGGTTGCACCTCAGGGGGCTGCAGAAGTCGCATCTGCCAGGCTGGCATAACTGCTGAACGc                                                    <  2:309615/149‑1 (MQ=255)                                                                                                  aaCGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGc                                                    <  1:875972/149‑1 (MQ=255)                                                                                                      ttATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTc                                               <  1:378188/150‑1 (MQ=255)                                                                                                      ttATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTc                                               <  1:686134/150‑1 (MQ=255)                                                                                                      ttATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCCTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGt                                                >  2:798051/1‑149 (MQ=255)                                                                                                      ttATCGTGCTTCGTCTGAT‑CCCGGTGTGTACCTGTTACGCGTGCCACTAGAATCGTATGCGTCGCCGTGCAGCCGTTGAAGCTCCGGCTGAACCGCAGGTGACTTCATAAGACGCATCTGCCAGCCTGTCAGAACTGCTGTACGCAGGTc                                               <  2:514010/150‑1 (MQ=255)                                                                                                        aTCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGc                                                                                                                                         >  2:483210/1‑58 (MQ=255)                                                                                                        aTCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGc                                                                                                                                         <  1:483422/58‑1 (MQ=255)                                                                                                        aTCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGGGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTgg                                            >  2:526061/1‑151 (MQ=255)                                                                                                          cGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTggg                                           >  1:796711/1‑150 (MQ=255)                                                                                                          cGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGc                                          <  1:945529/151‑1 (MQ=255)                                                                                                          cGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGc                                          <  1:1007891/151‑1 (MQ=255)                                                                                                                      gAT‑CCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATaa                               >  1:930665/1‑150 (MQ=255)                                                                                                                             ggCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTaa                        >  2:293149/1‑151 (MQ=255)                                                                                                                             ggCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTaa                        >  1:847232/1‑151 (MQ=255)                                                                                                                                 ggTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGt                    >  2:1029060/1‑151 (MQ=255)                                                                                                                                   tACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGtt                   <  1:1025090/150‑1 (MQ=255)                                                                                                                                      cGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGACAGAACTGCTGAACGCAGGTCTGGGGGGTTCTGATAACGAGTAATCGTTAATc               <  2:724684/151‑1 (MQ=255)                                                                                                                                                  ccactaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGt   >  2:459957/1‑151 (MQ=255)                                                                                                                                                  ccactaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGt   >  1:912604/1‑151 (MQ=255)                                                                                                                                                  ccactaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGt   >  1:855233/1‑151 (MQ=255)                                                                                                                                                  ccactaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGt   >  2:973797/1‑151 (MQ=255)                                                                                                                                                  ccactaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGGTCTGATAACGAGTAATCGTTAATCCGCAAATAACGt   >  2:760331/1‑151 (MQ=255)                                                                                                                                                  ccactaGGATCGTATGCGTCGCCGTGCTGCGGGCGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGt   >  2:934100/1‑151 (MQ=255)                                                                                                                                                   cactaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTa  >  1:321679/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAT‑CCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTA  >  NC_000913/4189300‑4189595

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑