Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,645,315	A→G	9.1%	N223S (AAC→AGC)	PP_4954 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,645,315	0	A	G	9.1%	148.0 / 5.4	56	N223S (AAC→AGC)	PP_4954	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (26/24);  new base G (3/2);  total (29/27)
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.77e-01

CAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGC  >  NC_002947/5645173‑5645462                                                                                                                                               |                                                                                                                                                     cagAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCa                                                                                                                                                  >  6:9241/1‑147 (MQ=255) cagAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCa                                                                                                                                                  <  5:9241/147‑1 (MQ=255) cagAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGc                                                                                                                                                >  5:17126/1‑149 (MQ=255)    agCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCTGCGTTCAGCGCGGCC‑ac                                                                                                                                            <  5:106341/149‑2 (MQ=255)      ccTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGCGGTGGCCGCGCTGAACGCAGCC‑AAgg                                                                                                                                          <  1:8638/149‑1 (MQ=255)          gCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGgaga                                                                                                                                       <  4:283318/148‑1 (MQ=255)          gCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGgaga                                                                                                                                       <  6:172980/148‑1 (MQ=255)              gACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACTGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGatcatc                                                                                                                                  <  8:287798/149‑1 (MQ=255)               aCTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGgaga                                                                                                                                       <  4:263112/143‑1 (MQ=255)               aCTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGgaga                                                                                                                                       >  3:263112/1‑143 (MQ=255)                cTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATcc                                                                                                                                 <  7:276147/148‑1 (MQ=255)                      ttcGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCg                                                                                                                          >  5:90306/1‑149 (MQ=255)                          aGCTGCTCAGCTACGTGAAGCTCAAGTAAATCCTGTTTTTTACGGTCCGCTCGCTCGACGTCTAGGCAAGGCGCCTGGAACTGGCCCAGCAGAAGGGGAAGGTGGTGGAGGGGCTGATCGGTGCC‑AAGGAGATCTTCCTGGAGCGGGTg                                                                                                                      <  8:236219/149‑1 (MQ=255)                                         tGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGccaacc                                                                                                        >  6:60923/1‑148 (MQ=255)                                          gAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGccaacca                                                                                                       >  6:195988/1‑148 (MQ=255)                                             gttCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACt                                                                                                   >  6:90464/3‑149 (MQ=255)                                               tCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTc                                                                                                  >  2:43673/1‑148 (MQ=255)                                               tCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTc                                                                                                  >  8:252379/1‑148 (MQ=255)                                                 aaGGCCAGCCTGTTTGACCCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCCGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCg                                                                                                 >  1:129732/1‑147 (MQ=255)                                                 aaGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGa                                                                                                <  5:195988/148‑1 (MQ=255)                                                 aaGGCCAGCCTGTCTGATCCCGTCCGGCCGCTCGACCACGACGCTCGGCGCCTGGAACTGGCCCAGCAGTACGGCAAGGTGGTGGACGCGCTGAACGCAGCC‑AAGCAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGaa                                                                                               <  6:164/149‑1 (MQ=255)                                                    gCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGt                                                                                            <  5:138587/149‑1 (MQ=255)                                                        gCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGc                                                                                        <  5:60923/149‑1 (MQ=255)                                                                gAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCt                                                                                >  3:165729/1‑149 (MQ=255)                                                                  aCCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGCACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTc                                                                               >  4:172981/1‑148 (MQ=255)                                                                   ccGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTc                                                                               >  3:13846/1‑147 (MQ=255)                                                                    cGGTCCGGACGCTCGACGTCGAGGCACGGCGCCTGGAACTGTCCCAGGAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGTTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAg                                                                            <  6:90306/149‑1 (MQ=255)                                                                      gTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCGAGCCCACCTCGAAGGTTAGCCGCTGCCTCa                                                                              >  8:211064/1‑145 (MQ=255)                                                                      gTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCa                                                                              <  7:211064/145‑1 (MQ=255)                                                                      gTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCt                                                                          >  8:206521/1‑149 (MQ=255)                                                                         cGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGt                                                                        <  1:43673/148‑1 (MQ=255)                                                                          gCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTAc                                                                      <  3:294666/149‑1 (MQ=255)                                                                              gCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGccaa                                                                                                          <  5:201418/109‑1 (MQ=255)                                                                              gCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGccaa                                                                                                          >  6:201418/1‑109 (MQ=255)                                                                                tCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTg                                                                >  1:312403/1‑149 (MQ=255)                                                                                       cGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTTGCCCAgg                                                         <  7:206521/149‑1 (MQ=255)                                                                                          ggCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGt                                                                        >  5:204582/1‑131 (MQ=255)                                                                                          ggCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGt                                                                        <  6:204582/131‑1 (MQ=255)                                                                                             aCGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCTAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATc                                                    >  1:249052/1‑149 (MQ=255)                                                                                                  gcCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAg                                              >  3:88431/1‑149 (MQ=255)                                                                                                      ggAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCgtg                                          >  5:141565/1‑149 (MQ=255)                                                                                                       gAACTGGCCCAGCAGAACGGCAAGGTTGTGGCCGCGCCGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCgt                                           <  2:129732/147‑1 (MQ=255)                                                                                                               ccAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCt                                 <  4:165729/149‑1 (MQ=255)                                                                                                                    agAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCAc                            <  5:90464/149‑1 (MQ=255)                                                                                                                    agAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCAc                            <  4:13846/149‑1 (MQ=255)                                                                                                                              ggtggtGGCCGCGCCGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTAc                   >  2:54891/1‑148 (MQ=255)                                                                                                                               tccttggCTGCGTTCAGCGCGGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAc                         <  7:27213/137‑1 (MQ=255)                                                                                                                               tccttggCTGCGTTCAGCGCGGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAc                         >  8:27213/5‑141 (MQ=255)                                                                                                                                tggtggCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAAc                >  7:9925/1‑149 (MQ=255)                                                                                                                                  ttggCTGCGTTCAGCGCGGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg               >  2:97939/2‑148 (MQ=255)                                                                                                                                  ttggCTGCGTTCAGCGCGGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg               >  4:136751/2‑148 (MQ=255)                                                                                                                                  gtggCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg               >  2:310520/1‑148 (MQ=255)                                                                                                                                    ggCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCAc            >  5:71656/1‑149 (MQ=255)                                                                                                                                       cgcgCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCa                                       >  4:11081/1‑119 (MQ=255)                                                                                                                                       cgcgCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCa                                       <  3:11081/119‑1 (MQ=255)                                                                                                                                       cgcgCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGa         <  4:88431/149‑1 (MQ=255)                                                                                                                                              aaCGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGc  <  2:249052/149‑1 (MQ=255)                                                                                                                                              aaCGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGc  <  6:141565/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                     CAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGC  >  NC_002947/5645173‑5645462

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