Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,555,068	G→A	11.8%	P380P (CCG→CCA)	sbcB →	exodeoxyribonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,555,068	0	G	A	11.8%	151.0 / 9.5	60	P380P (CCG→CCA)	sbcB	exodeoxyribonuclease I
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base A (3/4);  total (27/33)
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.95e-02

AACAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGA  >  NC_002947/1554923‑1555210                                                                                                                                                  |                                                                                                                                               aaaaGGAGCAATGGCAAGAAAAGCTGGAGAACATCTATGGCAAGGACGAGTTAGCCCATAGCGAGGACCCGGAACAGCAGTTTTATGAGGGAGTTATAGGTGACCGAGACCGTCTCTAATTCGATCAAATTCGTGCGCTAGAACCGgct                                                                                                                                             >  8:69621/4‑148 (MQ=255)  aCAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGgcg                                                                                                                                             >  5:287076/1‑148 (MQ=255)   cAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCCGGTTCCAGCGCa                                                                                                                                           <  7:286359/149‑1 (MQ=255)   cAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCCGGTTCCAGCGCa                                                                                                                                           <  6:35522/149‑1 (MQ=255)    aGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCATGCGCCGGTTCCAGCGCac                                                                                                                                          <  1:213867/149‑2 (MQ=255)      gcgcAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGtt                                                                                                                                        <  5:112207/149‑1 (MQ=255)        gcAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGa                                                                                                                                      <  5:210264/149‑1 (MQ=255)          aaTGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGc                                                                                                                                    <  7:50695/149‑1 (MQ=255)              gCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACCGCAGTTGTATGAGGGTTTTATAGGCGACCGCGCCCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATg                                                                                                                                <  4:21074/149‑1 (MQ=255)               cAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATgg                                                                                                                               >  7:173899/1‑149 (MQ=255)                aaGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGc                                                                                                                              >  3:118730/1‑149 (MQ=255)                aaGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGc                                                                                                                              <  4:38581/149‑1 (MQ=255)                       gCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGa                                                                                                                       <  7:304217/149‑1 (MQ=255)                                 aTCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCgatga                                                                                                             >  5:306486/1‑149 (MQ=255)                                   cTATGGCAAGGACGAGTTCGCCACGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACAGTCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGCCCAGAGGATGtt                                                                                                                   >  2:236521/1‑141 (MQ=255)                                   cTATGGCAAGGACGAGATCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCATGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGtt                                                                                                                   <  1:236521/141‑1 (MQ=255)                                       ggCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCcgcg                                                                                                       >  3:187998/1‑149 (MQ=255)                                          aaGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGttt                                                                                                    <  6:287076/149‑1 (MQ=255)                                            ggACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGc                                                                                                  >  3:244738/1‑149 (MQ=255)                                            ggACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGc                                                                                                  >  1:34092/1‑149 (MQ=255)                                              aCGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCg                                                                                                <  6:306486/149‑1 (MQ=255)                                                 aGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGa                                                                                                                                      <  6:125710/108‑1 (MQ=255)                                                 aGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGa                                                                                                                                      >  5:125710/1‑108 (MQ=255)                                                      gCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCtgtt                                                                                        >  4:301517/1‑149 (MQ=255)                                                         ccGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTGTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTgttgtt                                                                                     <  2:255970/149‑1 (MQ=255)                                                          cgagcgagGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTc                                                                                    <  4:138027/149‑1 (MQ=255)                                                                     ccGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCggg                                                                          >  1:264870/1‑148 (MQ=255)                                                                            agcagTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGATGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTacc                                                                  >  8:71932/1‑147 (MQ=255)                                                                                gTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTcc                                                              <  4:89648/149‑1 (MQ=255)                                                                                   gTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCt                                                            >  3:18680/1‑148 (MQ=255)                                                                                    tATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTg                                                           >  7:27238/1‑148 (MQ=255)                                                                                    tATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGa                                                          <  2:150269/149‑1 (MQ=255)                                                                                     aTGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGCGCCATGGCCAGTGGATGTTCGATGATCCGCGTTGGCCGGAGATGTTGTTCCGCTACCGGGCCGGTAGCTTccaggaa                                                         >  8:54898/1‑144 (MQ=255)                                                                                          ggTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCtt                                                    >  5:159660/1‑149 (MQ=255)                                                                                            ttttATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGa                                                  <  4:24574/149‑1 (MQ=255)                                                                                                 taGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCgg                                              <  3:301517/148‑1 (MQ=255)                                                                                                  aGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGc                                             >  8:19038/1‑148 (MQ=255)                                                                                                     tGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCg                                            >  8:189435/1‑146 (MQ=255)                                                                                                         cgcgACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCg                                     >  8:299418/1‑149 (MQ=255)                                                                                                           cgACCGGCGCTTATGCGAGCAAGGTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGc                                   <  8:173899/149‑1 (MQ=255)                                                                                                               cGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCAAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACg                               >  1:86139/1‑149 (MQ=255)                                                                                                                ggCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAAc                                                       <  3:170195/124‑1 (MQ=255)                                                                                                                ggCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAAc                                                       >  4:170195/1‑124 (MQ=255)                                                                                                                  cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCTAAGAGCGGCAACGTTg                            <  4:244738/149‑1 (MQ=255)                                                                                                                  cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                            <  7:189435/149‑1 (MQ=255)                                                                                                                  cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                            <  5:246594/149‑1 (MQ=255)                                                                                                                  cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                            <  2:34092/149‑1 (MQ=255)                                                                                                                  cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                            <  4:118730/149‑1 (MQ=255)                                                                                                                  cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGAGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTg                            <  7:299418/149‑1 (MQ=255)                                                                                                                   gcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTgg                           <  4:187998/149‑1 (MQ=255)                                                                                                                    cTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGt                          <  4:18680/149‑1 (MQ=255)                                                                                                                                 tcaactGCGCCGGTTCCAGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCg                                     <  2:192472/120‑1 (MQ=255)                                                                                                                                 tcaactGCGCCGGTTCCAGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGcc              >  2:191317/6‑148 (MQ=255)                                                                                                                                 tcaactGCGCCGGTTCCAGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGcc              >  6:135481/6‑148 (MQ=255)                                                                                                                                 gTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTcc                                                              <  2:279495/100‑1 (MQ=255)                                                                                                                                 gTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTcc                                                              >  1:279495/1‑100 (MQ=255)                                                                                                                                  caactGCGCCGGTTCCAGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCaa            >  7:230324/5‑149 (MQ=255)                                                                                                                                      tGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCt                                                            >  7:284523/1‑97 (MQ=255)                                                                                                                                      tGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCt                                                            <  8:284523/97‑1 (MQ=255)                                                                                                                                             gAACCTGCTCAGTTTAGCCATGGCCAGTGGATGTTCGATTATCCGCGTTTGCCGGAGCTGTTGTTCCGTTACGGGGGACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCtta  <  7:69621/148‑3 (MQ=255)                                                                                                                                                  |                                                                                                                                               AACAGGCGCAATGGCAAGAAAAGCTGGAGCACATCTATGGCAAGGACGAGTTCGCCCCGAGCGAGGACCCGGAACAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGA  >  NC_002947/1554923‑1555210

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