Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	85,091	Δ1 bp	7.3%	intergenic (+69/+61)	aroE → / ← PP_0075	shikimate 5‑dehydrogenase/choline sulfate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	85,091	0	A	.	7.3%	141.4 / 6.9	55	intergenic (+69/+61)	aroE/PP_0075	shikimate 5‑dehydrogenase/choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base A (26/24);  new base . (2/2);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

TCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAGGTC  >  NC_002947/84952‑85239                                                                                                                                             |                                                                                                                                                       tCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgt                                                                                                                                                <  8:84725/148‑1 (MQ=255) tCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgt                                                                                                                                                <  8:79743/148‑1 (MQ=255)     gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgtg                                                                                                                                               <  1:360221/145‑1 (MQ=255)     gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAg                                                                                                                                           >  6:285450/1‑149 (MQ=255)       gtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCC‑ATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACG‑TACTTTGTGGGAACgg                                                                                                                                        <  1:238685/149‑1 (MQ=18)                 gCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTc                                                                                                                               >  7:272881/1‑149 (MQ=255)                    ggATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTg                                                                                                                                   <  6:142578/142‑1 (MQ=255)                    ggATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTg                                                                                                                                   >  5:142578/1‑142 (MQ=255)                    ggATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTcgc                                                                                                                             >  8:275297/1‑148 (MQ=255)                                ggTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCaaa                                                                                                                <  5:285450/149‑1 (MQ=255)                                ggTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCa                                                                                                                  >  6:331222/1‑147 (MQ=255)                                        cTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTAGGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAc‑‑cc                                                                                                      <  6:176582/149‑1 (MQ=255)                                            aCTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCTTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGaa                                                                                                  <  7:1164/149‑1 (MQ=255)                                              tGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATc                                                                                                >  4:40118/1‑149 (MQ=255)                                               gcgcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGGGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCa                                                                                                             >  6:162708/1‑137 (MQ=255)                                               gcgcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGAGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCa                                                                                                             <  5:162708/137‑1 (MQ=255)                                                cgcCGGCAACTGTCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑GGTCGTACGCCGATCCTGTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATc                                                                                                <  4:222366/147‑1 (MQ=255)                                                cgcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATcc                                                                                               >  5:66270/1‑148 (MQ=255)                                                        aCTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCAt                                                                                      <  5:269369/149‑1 (MQ=255)                                                           cgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCt                                                                                         >  7:256806/1‑143 (MQ=255)                                                           cgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCt                                                                                         <  8:256806/143‑1 (MQ=255)                                                                ggggCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCaaa                                                                              <  6:66270/149‑1 (MQ=255)                                                                            gCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCggg                                                                    <  3:40118/147‑1 (MQ=255)                                                                             ccGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                                  >  7:146094/1‑148 (MQ=255)                                                                              tagggCCGCCAAGCGGCCCTATCGCGACGCAACGCCGTTACCACAAAGTA‑CGTCGTACGCCGCTCCACTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAACCAGAC‑‑CCAGAATCCCCACCTCAATCCTCAACCCCGATCGGGcaca                                                                >  8:139354/3‑149 (MQ=37)                                                                              cGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGc                                                                                                                                  >  8:212209/1‑85 (MQ=255)                                                                                gggCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAAc                                                              <  8:272881/149‑1 (MQ=255)                                                                                   ccGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTcc                                                           <  8:146094/149‑1 (MQ=255)                                                                                      cTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTcc                                                                                   <  8:243762/122‑1 (MQ=255)                                                                                      cTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTcc                                                                                   >  7:243762/1‑122 (MQ=255)                                                                                              gcccaTCGCGACGCAAGGCCGCGCCCACAGAGGATAGGCGTAAGCCG‑TACGTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGCATCCCAACCGAATTCCTCAAACCGGATCGGGCACCACTCCACCCCGTcccg                                                >  4:88614/5‑147 (MQ=9)                                                                                               ccTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAAc                                                                             <  7:214385/119‑1 (MQ=255)                                                                                               ccTATCGCGACGCAAGGCCGTTACCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAAc                                                                             >  8:214385/1‑119 (MQ=255)                                                                                                cTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCt                                              >  3:260588/1‑149 (MQ=255)                                                                                                   tCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCt                                           >  7:204603/1‑149 (MQ=255)                                                                                                     gcgACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTg                                          >  4:124887/1‑148 (MQ=255)                                                                                                      cgACGCAAGGCCGTACCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGTGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCACAATCCACAAACCCGATCGGGCACAACTCCACCCCTTCCAGCGTCTGc                                         >  8:176599/1‑148 (MQ=39)                                                                                                                 cgctCCCACAGAGGATCGGCGTACGACG‑TACTTTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAc                                                         <  6:143270/118‑1 (MQ=18)                                                                                                                 cgctCCCACAGAGGATCGGCGTACGACG‑TACTTTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAc                                                         >  5:143270/4‑121 (MQ=18)                                                                                                                  gTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc                                         <  3:185698/136‑1 (MQ=255)                                                                                                                  gTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc                                         >  4:185698/1‑136 (MQ=255)                                                                                                                   ttCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAc                           >  3:293917/1‑149 (MQ=255)                                                                                                                     cccACAAAGTA‑CGTCGTACGCCGATCCTCTGGGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAcc                          <  7:139354/148‑1 (MQ=255)                                                                                                                      ccACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGa                               <  3:106905/142‑1 (MQ=255)                                                                                                                      ccACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGa                               >  4:106905/1‑142 (MQ=255)                                                                                                                      ccACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTg                        >  5:326430/1‑149 (MQ=255)                                                                                                                       cacaAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGAGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGc                       >  1:44496/1‑149 (MQ=255)                                                                                                                        acaAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCg                      >  6:29826/1‑149 (MQ=255)                                                                                                                        acaAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCg                      <  8:204603/149‑1 (MQ=255)                                                                                                                             gTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACTACCAGAATCCCAACCTCATTCCTCAAACCGGAACGGGCACACCTCCACCCATGCCACGTGGTGCACCACCACGATCACATGCGGca                   >  1:103305/1‑148 (MQ=25)                                                                                                                             gTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTg                 <  4:260588/149‑1 (MQ=255)                                                                                                                                         gCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGc        >  4:176508/1‑147 (MQ=255)                                                                                                                                         gCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCaa      <  3:124887/149‑1 (MQ=255)                                                                                                                                           cGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCaa      >  3:191242/1‑147 (MQ=255)                                                                                                                                           cGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAgg    <  4:293917/149‑1 (MQ=255)                                                                                                                                             aTCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAGGTc  <  5:331222/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                       TCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCGGTGCAAGGTC  >  NC_002947/84952‑85239

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