Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,988,605	A→G	7.3%	T375A (ACC→GCC)	oplA →	5‑oxoprolinase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,988,605	0	A	G	7.3%	148.6 / 2.8	55	T375A (ACC→GCC)	oplA	5‑oxoprolinase A
Reads supporting (aligned to +/- strand):  ref base A (27/24);  new base G (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 = 8.33e-01

CCGATGATCCAGGTGGAGACCCTGGGTGCCGGCGGTGGCTCCATCGCCCACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATCGCCGAACCGCTGGGCATCAGTGTGGAACGCGCGGCATTCGGCAT  >  NC_002947/3988458‑3988753                                                                                                                                                    |                                                                                                                                                       ccGATGATCCAGGTGTAGACTCTGGGTGCCGGCGGTGTCTCCATCGCCCACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGGAGCGCCGGCGCCATGCCTGTGCCGGTGTGCTACGGAAATGGTGGCCGCGAACCCACGGTCAc                                                                                                                                                       <  2:153614/149‑1 (MQ=255)   gatgatCCAGGTGGAGACCCTGGGTGCCGGCGGTGGCTCCATCGCCCACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCAcc                                                                                                                                                      >  1:334197/1‑148 (MQ=255)           aGGTGGAGACCCTGGGTGCCGGCGGTGGCTCCATCGCCCACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAAc                                                                                                                                             <  2:58868/149‑1 (MQ=255)                    cccTGGGTGCCGGCGGTGGCTCCATCGCCCACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCtt                                                                                                                                    <  1:130830/149‑1 (MQ=255)                    cccTGGGTGCCGGCGGTGGCTCCATCGCCCACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCtt                                                                                                                                    <  4:324291/149‑1 (MQ=255)                    cccTGGGTGCCGGCGGTGGCTCCATCGCCCACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCtt                                                                                                                                    <  1:118623/149‑1 (MQ=255)                     ccTGGGTGCCGGCGGTGGCTCCATCGCCCACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTg                                                                                                                                   <  3:117993/149‑1 (MQ=255)                     ccTGGGTGCCGGCGGTGGCTCCATCGCCCACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTg                                                                                                                                   <  1:56262/149‑1 (MQ=255)                          gTGCCGGCGGTGGCTCCATCGCCCACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTAc                                                                                                                              <  8:182689/149‑1 (MQ=255)                                gcggTGGCTCCATCGCCCACGTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGGCGGCGCCATGCCGGGGCCGGTGTGCTACGGCAAGGGGGGCCGCGAACCCACGGTCACCGTCGCCAACCTGGCTCTTGGCTACCTGccc                                                                                                                        <  8:78979/149‑1 (MQ=255)                                   gTGGCTCCATCGCCCACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGAc                                                                                                                     <  6:303161/149‑1 (MQ=255)                                    tGGCTCCATCGCCCACCTGGACGAGGTTGGCAGGCACCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTAACCGACGCCAACCTGGCTGTTCGCTACCTGCCCGACg                                                                                                                    >  6:178698/1‑149 (MQ=255)                                      gCTCCATCGCCCACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACgg                                                                                                                   <  3:39261/148‑1 (MQ=255)                                       cTCCATCGCCCACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGgcg                                                                                                                 >  1:5979/1‑149 (MQ=255)                                                ccACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCGGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGGCACCGACGCCATCCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTgg                                                                                                       <  5:178698/149‑1 (MQ=255)                                                  aCCTGGACGAGTTTGACATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGggc                                                                                                     >  1:185724/1‑149 (MQ=255)                                                    cTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCt                                                                                                                                <  2:299462/121‑1 (MQ=255)                                                    cTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCt                                                                                                                                >  1:299462/1‑121 (MQ=255)                                                     tGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGggcgg                                                                                                   >  1:145062/1‑148 (MQ=255)                                                       gACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAg                                                                                                >  1:339746/1‑149 (MQ=255)                                                            gTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCGCCGACGCCACCCTGGCTCGTGGCTAc                                                                                                                              >  2:257181/1‑115 (MQ=255)                                                            gTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCGCCGACGCCAACCTGGCTCTTGGCTAc                                                                                                                              <  1:262106/115‑1 (MQ=255)                                                            gTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCGCCGACGCCAACCTGGCTCTTGGCTAc                                                                                                                              >  2:262106/1‑115 (MQ=255)                                                            gTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATcc                                                                                           >  5:140857/1‑149 (MQ=255)                                                            gTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCGGGGCCGGTGTGCTACGGCAAGGGGGGCCGCGAACCCACGGTCGCCGACGCCAACCTGGCTCTTGGCTAc                                                                                                                              <  1:257181/115‑1 (MQ=255)                                                             tttGGCATGCTCCAGGTCGGCCCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCg                                                                                          >  4:91287/1‑149 (MQ=255)                                                                   aTGCTCCAGGTCGGAACGCGCAGCGACGGAGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCTAACCCCCGTTCACAGACGCCAACCTGGATCTTTGCTACCTGCCCGCCGGCGC‑AATGCTGGGCGGAAGTATCCGAATGaa                                                                                    >  3:211821/1‑149 (MQ=255)                                                                       tCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGc                                                                                >  3:260751/1‑149 (MQ=255)                                                                                gACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCa                                                                                                 >  4:36549/1‑123 (MQ=255)                                                                                gACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCa                                                                                                 <  3:36549/123‑1 (MQ=255)                                                                                 aCCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCgcgc                                                                      >  5:243982/1‑149 (MQ=255)                                                                                   cgcgCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTg                                                                    >  1:142678/1‑149 (MQ=255)                                                                                   cgcgCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTg                                                                    <  2:185724/149‑1 (MQ=255)                                                                                      gcAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAc                                                                 >  3:12341/1‑149 (MQ=255)                                                                                            ccGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATc                                                          >  3:10386/1‑149 (MQ=255)                                                                                               gcgcCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCgg                                                       <  6:243982/149‑1 (MQ=255)                                                                                                cgcCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGa                                                      <  4:260751/149‑1 (MQ=255)                                                                                                  ccATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTAGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAg                                                     >  8:56686/1‑148 (MQ=255)                                                                                                        cTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGc                                                    >  5:79371/1‑143 (MQ=255)                                                                                                        cTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATc                                              <  2:339746/149‑1 (MQ=255)                                                                                                             ccGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATCGCCg                                          >  2:236611/1‑148 (MQ=255)                                                                                                                  gtgCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATCGCCGAACCGc                                    <  4:12341/149‑1 (MQ=255)                                                                                                                      tACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGa                                                      >  4:244191/1‑127 (MQ=255)                                                                                                                        cGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATCGCCGAACCGCTggg                                <  1:236611/147‑1 (MQ=255)                                                                                                                        cGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATCGCCGAACCGCTGGGCa                              >  7:80214/1‑149 (MQ=255)                                                                                                                           cAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTAGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATCGCCGAACCGCTGGGCATCa                           <  7:56686/149‑1 (MQ=255)                                                                                                                            aaGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATTGCCGAACCGCTGGGCATCAg                          >  2:249375/1‑149 (MQ=255)                                                                                                                               ggTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATCGCCGAACCGCTGGGCATCAgtgt                       <  6:79371/149‑1 (MQ=255)                                                                                                                                 tGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGACGCGCTGGCC‑GCCATCCGGAGCAAGATCGCCGAACCGCTGGGCATCAGtgtg                      >  2:108500/1‑148 (MQ=255)                                                                                                                                   gCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCGACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGACTGCCATCCGGAGCAAGATTGCCGAACCGCTGGGCATCAGTGTgg                     <  1:249375/148‑1 (MQ=255)                                                                                                                                     cgcgAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGCTACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATCGCCGAACCGCTGGGCATCAGTGTGGAAc                  >  1:278844/1‑149 (MQ=255)                                                                                                                                      gcgAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATCGCCGAACCGCTGGGCATCAGTGTGGAAcgc                >  3:326856/1‑149 (MQ=255)                                                                                                                                       cgAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATCGCCGAACCGCTGGGCATCAGTGTGGAAcgcg               <  8:80214/149‑1 (MQ=255)                                                                                                                                            ccACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTTCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATCGCCGAACCGCTGGGCATCAGTGTGGAACGCGCGGCa          >  5:345026/1‑149 (MQ=255)                                                                                                                                                 gTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATCGCCGAACCGCTGGGCATCAGTGTGGAACGCGCGGCATTCgg     >  1:147844/1‑149 (MQ=255)                                                                                                                                                    aCCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGTCGGCAGTATTCGCCTGAACCGCCAGGCCGCTCTGGAC‑GCCATCCGTAGCAAGATCGCCGAACCGCTGGTCATGAGTGTGGAACGCGCGGCATTCGGCAt  <  8:1727/149‑1 (MQ=255)                                                                                                                                                    aCCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATCGCCGAACCGCTGGGCATCAGTGTGGAACGCGCGGCATTCGGCAt  <  6:140857/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       CCGATGATCCAGGTGGAGACCCTGGGTGCCGGCGGTGGCTCCATCGCCCACCTGGACGAGTTTGGCATGCTCCAGGTCGGACCGCGCAGCGCCGGCGCCATGCCTGGGCCGGTGTGCTACGGCAAGGGTGGCCGCGAACCCACGGTCACCGACGCCAACCTGGCTCTTGGCTACCTGCCCGACGGCGC‑ACTGCTGGGCGGCAGTATCCGCCTGAACCGCCAGGCCGCGCTGGAC‑GCCATCCGGAGCAAGATCGCCGAACCGCTGGGCATCAGTGTGGAACGCGCGGCATTCGGCAT  >  NC_002947/3988458‑3988753

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