Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,499,480	2 bp→AC	intergenic (+57/‑106)	trpS → / → zapE	tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,499,480	0	C	A	100.0%	152.8 / NA	53	intergenic (+57/‑107)	trpS/zapE	tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (27/26);  total (27/26)
*	NC_002947	1,499,481	0	A	C	100.0%	153.8 / NA	53	intergenic (+58/‑106)	trpS/zapE	tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (27/26);  total (27/26)

CGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCG‑GTGT‑T‑T‑CGGGCA‑CG‑CCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGA  >  NC_002947/1499347‑1499620                                                                                                                                      ||                                                                                                                                                   cGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑gtg                                                                                                                                     >  1:28336/1‑149 (MQ=255)    aGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGt‑tc                                                                                                                                 <  3:44694/149‑2 (MQ=255)    aGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCGTTGCCG‑GTGt‑tc                                                                                                                                 >  7:354884/1‑148 (MQ=255)     gTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGGGCCCGCGCAAAACGCGCCGCGGTGCCTGGCACCGGCGTTGGCG‑GGGt‑t‑tg                                                                                                                               >  3:184720/1‑148 (MQ=255)      ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGt‑tcgc                                                                                                                               >  8:235043/1‑146 (MQ=255)      ttCGACAACGCTGAAGCCAGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGCCGTTCCCGCGAAACAGGCCCCCCCGTGCATGGCACCCCCTTTGCCG‑GTGt‑tcgc                                                                                                                               >  2:328421/1‑146 (MQ=255)         gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCTTTGCCG‑GTGt‑tcgcggg                                                                                                                            <  7:279511/149‑7 (MQ=37)          aCAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGt‑tcgcgggc                                                                                                                           >  7:210883/1‑142 (MQ=37)            aaCGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGt‑tc                                                                                                                                 >  6:355918/1‑140 (MQ=255)            aaCGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGt‑tc                                                                                                                                 <  5:355918/141‑2 (MQ=255)              cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑cg‑c                                                                                                                    <  6:227143/149‑1 (MQ=37)                cTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCCCGGTGCATGGCACCCGCATTGCCG‑GTGT‑TCG‑CGGGCA‑cg‑ccc                                                                                                                  >  6:232215/1‑148 (MQ=37)                    aGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑cg‑cccgctc                                                                                                              <  6:143446/149‑6 (MQ=25)                                      tCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑At                                                                                            >  5:89525/1‑149 (MQ=21)                                             cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCgg                                                                                     <  7:235043/149‑1 (MQ=21)                                              tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGGCCGCGAAACAGGCGCCGCCGGGCATGGGACCGGGGTTGCCG‑GTGT‑T‑TGCGGGGACCG‑CCGCGCGCACACACGGAGTGGC‑AGTTAGgggg                                                                                    >  8:184573/1‑149 (MQ=255)                                                aaCGAAGCCTTGGCACCCCAGCAGGACTGAGGCTGCCAAACGGTGGGTGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCATGGCACCGGGTTTGCCGCGTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCgg                                                                                     <  1:328421/147‑1 (MQ=21)                                                 aCGAAGCCTTGTCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCt                                                                                 <  1:154173/149‑1 (MQ=21)                                                 aCGAAGCCTGGCCACCACCGCCGCACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCTAAACCGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCt                                                                                 <  3:214173/149‑1 (MQ=18)                                                       ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCCCAAAACGCGCCCCGGGGCATGGCACCCGCGTTTCCC‑GTGT‑TAG‑CGGGGA‑CGCCCGCACCCCCAAACAGAATTGC‑AATAAACCGGCCGATCGc                                                                            >  4:156209/1‑148 (MQ=255)                                                        cTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGTCGCCGCTGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGc                                                                                   <  5:317120/140‑1 (MQ=21)                                                               ccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCCGCTTTGCCG‑GTGT‑TAG‑CGGGCA‑CGCCCGCGCCCCCAGCCTGATTGCC‑ATTAATCGGGCCGATCGCTACAGTGAc                                                                   >  2:294059/1‑149 (MQ=21)                                                                  ctgcagGACGGCGGCTGCCAAACGGTGGGAGCGGGCGGGCCCGCGAAACCGGCGCCGCGGGGCATGGCACCGGCTTTGCCG‑GTGTGG‑G‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGc                                                                <  4:317636/147‑1 (MQ=21)                                                                   agcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAACGCGCCGCGGTGCAAGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGGC‑CGCCCGCGCCCACAGCCTGAATGGC‑ATTAAGCGGGGCGATCGCTAAAGTGAAGGcc                                                               >  2:129224/1‑149 (MQ=12)                                                                    gcagGACTGCGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGGGGTGCATGGCACGGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                              <  2:148559/149‑1 (MQ=21)                                                                             aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCCCAGACTGATTGCC‑ATTAAGCGGGCCTATCGCGACAGTGACGGCCCGTGTTttgg                                                     >  8:296967/1‑148 (MQ=21)                                                                             aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                       <  8:210883/147‑1 (MQ=21)                                                                              ggcggcCAAACGGTGGTCGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCAGGGCACCGGCTTTGCCG‑GTGTGC‑G‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                    <  3:174555/149‑1 (MQ=25)                                                                                    ctatcGGTGGGAGCGGGGGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCg                                              <  1:408250/145‑1 (MQ=21)                                                                                        cGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTaa                                          <  2:28336/149‑1 (MQ=21)                                                                                             gtaGCGGGCGGGCCCGCTAAACTGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑GGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                                     <  8:354884/147‑1 (MQ=21)                                                                                               aGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TGG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                   <  5:232215/149‑1 (MQ=25)                                                                                                gCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                  <  5:341876/149‑1 (MQ=25)                                                                                                 cGGGGGTGCCCGCTAAACAGGGGCCGCGGTGCATGGCACCGGCTTTTCCG‑GTGT‑TCG‑CGGTAA‑CGCCCGCTCCCAAAGACTGATTGCCTCTTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                  <  1:129224/149‑1 (MQ=21)                                                                                                  gggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                                <  3:156209/149‑1 (MQ=37)                                                                                                         cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                          <  6:89525/148‑1 (MQ=37)                                                                                                          ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                          >  2:146626/1‑147 (MQ=37)                                                                                                          ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                          >  3:157894/1‑147 (MQ=37)                                                                                                           cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑CTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTCACGGCCCGTGTTTAGTTGCCTCGATAACGAATAATGACTCCCt                          >  4:317578/1‑146 (MQ=21)                                                                                                            gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                       >  3:250945/1‑148 (MQ=37)                                                                                                            gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                       >  3:364154/1‑148 (MQ=37)                                                                                                               aaaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                        <  5:428203/144‑1 (MQ=37)                                                                                                               aaaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                        >  6:428203/1‑144 (MQ=37)                                                                                                                aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGc                  >  5:197704/1‑149 (MQ=37)                                                                                                                 aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGc                  >  4:394060/1‑148 (MQ=37)                                                                                                                  cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa                >  2:270427/1‑149 (MQ=37)                                                                                                                   aGGCGCCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa                >  4:15621/1‑148 (MQ=37)                                                                                                                     gcgcCGCGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc              >  4:285019/1‑148 (MQ=37)                                                                                                                       gccgcGGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg           <  4:250945/149‑1 (MQ=37)                                                                                                                         cgcgGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg           >  7:48621/1‑147 (MQ=37)                                                                                                                         cgcgGTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCa         <  4:157894/149‑1 (MQ=37)                                                                                                                             gTGCATGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATc     <  1:294059/149‑1 (MQ=38)                                                                                                                                 aTGGCACCGGCTTTGCCG‑GTGT‑TCG‑CGGGCA‑CGCCCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGCGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGa  >  6:18362/1‑148 (MQ=37)                                                                                                                                      ||                                                                                                                                                   CGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCG‑GTGT‑T‑T‑CGGGCA‑CG‑CCGCTCCCACAGACTGATTGCC‑ATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGA  >  NC_002947/1499347‑1499620

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