Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,587,622:1	+A	7.0%	coding (458/915 nt)	catA‑II →	catechol 1,2‑dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,587,622	1	.	A	7.0%	223.1 / 10.6	57	coding (458/915 nt)	catA‑II	catechol 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base . (26/27);  new base A (2/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

ACCATCGAAGGCCCGTTGTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGC  >  NC_002947/3587474‑3587767                                                                                                                                                      |                                                                                                                                                  aCCATCGAAGGCCCGTTGTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGACAGGTGCTGGCCCCGCACGGACGCCCGCTGCCGGGTGCCAc                                                                                                                                                    >  3:286139/1‑149 (MQ=255)          ggCCCGTTGTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCTGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCt                                                                                                                                          >  4:88812/1‑149 (MQ=255)           gCCCGTTGTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCTGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCtg                                                                                                                                         <  3:88812/149‑1 (MQ=255)           gCCCGTTGTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACAGGTCGAccg                                                                                                                                          >  2:291943/1‑148 (MQ=255)           gCCCGTTGTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCtg                                                                                                                                         >  2:202486/1‑149 (MQ=255)            cccGTTGTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCtgt                                                                                                                                        <  4:371485/149‑1 (MQ=255)               gTTGTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTgg                                                                                                                                      >  4:69781/1‑148 (MQ=255)               gTTGTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGc                                                                                                                                     <  2:230549/149‑1 (MQ=255)               gTTGTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGc                                                                                                                                     <  5:167072/149‑1 (MQ=255)                  gTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACAGGTCGACC‑GTGGCac                                                                                                                                   <  1:291943/148‑2 (MQ=255)                   tACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGTCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATg                                                                                                                                  <  2:120578/148‑1 (MQ=255)                             gcgcACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACAGGTCGACC‑GTGGCacc                                                                                                                                  <  5:93852/138‑3 (MQ=255)                             gcgcACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCACAGGTCGACC‑GTGGCacc                                                                                                                                  >  6:93852/1‑136 (MQ=255)                             gcgcACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGt                                                                                                                       <  3:69781/149‑1 (MQ=255)                               gcACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCAGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTgg                                                                                                                     >  8:121661/1‑149 (MQ=255)                                   cGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTAcc                                                                                                                 >  1:29888/1‑149 (MQ=255)                                         cccAAGGTGAAGTGCGCATGGACGACGGCAGCGAAGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCg                                                                                                           >  8:46781/1‑149 (MQ=255)                                           cAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTg                                                                                                                      >  7:322829/1‑136 (MQ=255)                                           cAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTg                                                                                                                      <  8:322829/136‑1 (MQ=255)                                             aGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGttct                                                                                                       <  4:26887/149‑1 (MQ=255)                                                  aaGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGa                                                                                                   <  1:202486/148‑1 (MQ=255)                                                       cgcATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCaaa                                                                                              <  7:338141/148‑1 (MQ=255)                                                       cgcATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCaaa                                                                                              <  2:29888/148‑1 (MQ=255)                                                        gcATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGc                                                                                            >  7:139695/1‑149 (MQ=255)                                                           tGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAg                                                                                         >  8:229351/1‑149 (MQ=255)                                                                 acgGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCg                                                                                   <  1:226031/149‑1 (MQ=255)                                                                      aGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACaa                                                                              >  8:346959/1‑149 (MQ=255)                                                                      aGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACaa                                                                              <  7:229351/149‑1 (MQ=255)                                                                          aggaggGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTg                                                                          >  5:145408/1‑149 (MQ=255)                                                                          aggaggGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTg                                                                          >  5:365780/1‑149 (MQ=255)                                                                               ggCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTa                                                                                 >  4:33834/1‑137 (MQ=255)                                                                               ggCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTa                                                                                 <  3:33834/137‑1 (MQ=255)                                                                               ggCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCg                                                                     >  2:241467/1‑149 (MQ=255)                                                                               ggCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCg                                                                     <  4:73413/149‑1 (MQ=255)                                                                                 cGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCGGTCGGCGTACAACCTGCGTCGgc                                                                   <  8:139695/149‑1 (MQ=255)                                                                                 cGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGgc                                                                   <  7:149590/149‑1 (MQ=255)                                                                                     gCCACGGCGATGATCCTGGACGGCGAGGTGCTGGACCCGCACGGACGCCGGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAACGCCAGTCGGCGTACACCCTGCGTCGGCGCa                                                                <  6:21084/148‑1 (MQ=255)                                                                                        aCGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCaaa                                                                                              >  7:219410/1‑115 (MQ=255)                                                                                        aCGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCaaa                                                                                              <  8:219410/115‑1 (MQ=255)                                                                                                gTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATg                                                    >  5:133888/1‑149 (MQ=255)                                                                                                    cTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGGTAACGATGcac                                                 <  8:23539/148‑3 (MQ=255)                                                                                                    cTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGGTAACGATGcac                                                 >  7:23539/1‑146 (MQ=255)                                                                                                     tGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAg                                               >  5:236719/1‑149 (MQ=255)                                                                                                             cAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGc                                                                         >  3:334318/1‑115 (MQ=255)                                                                                                             cAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGc                                                                         <  4:334318/115‑1 (MQ=255)                                                                                                                  gCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACcgcg                                  <  7:346959/149‑1 (MQ=255)                                                                                                                   cTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACcgcgc                                 <  6:365780/149‑1 (MQ=255)                                                                                                                   cTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACcgcgc                                 <  6:236719/149‑1 (MQ=255)                                                                                                                   cTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACcgcgc                                 >  2:74939/1‑149 (MQ=255)                                                                                                                        cccGCACGGACGCCCGCTGCCGGGTGCCAC‑TGTCGACCTGTGTCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCt                            <  4:286139/149‑1 (MQ=255)                                                                                                                         ccGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTc                           >  8:173116/1‑149 (MQ=255)                                                                                                                             aCGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGCTGCCCAGGGGCGCTACCGCGCGCGCTCCATc                       >  8:222325/1‑149 (MQ=255)                                                                                                                             aCGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATc                       >  4:29044/1‑149 (MQ=255)                                                                                                                                 aCGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGc                   >  1:322151/1‑149 (MQ=255)                                                                                                                                     ccGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATc               >  3:71993/1‑149 (MQ=255)                                                                                                                                      cGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCg              >  1:78463/1‑149 (MQ=255)                                                                                                                                       gCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCgg             <  7:121661/149‑1 (MQ=255)                                                                                                                                                  ccAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGc  <  7:173116/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  ACCATCGAAGGCCCGTTGTACGTGGCTGGCGCACCGATTGCCCAAGGTGAAGTGCGCATGGACGACGGCAGCGAGGAGGGCGTGGCCACGGTGATGTTCCTGGAAGGCCAGGTGCTGGACCCGCACGGACGCCCGCTGCCGGGTGCCAC‑GGTCGACCTGTGGCATGCCAATACCCGTGGTACCTACTCGTTCTTCGACCAAAGCCAGTCGGCGTACAACCTGCGTCGGCGCATCGTTACCGATGCCCAGGGGCGCTACCGCGCGCGCTCCATCGTGCCATCGGGCTATGGCTGC  >  NC_002947/3587474‑3587767

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