Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,789,805	G→C	5.4%	L289L (CTG→CTC)	atsA →	arylsulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,789,805	0	G	C	5.4%	291.2 / 4.3	93	L289L (CTG→CTC)	atsA	arylsulfatase
Reads supporting (aligned to +/- strand):  ref base G (37/50);  new base C (3/2);  total (40/53)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CGTACCTGGCGCAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTGTTCATGTCCGACAA  >  NC_002947/3789657‑3789954                                                                                                                                                     |                                                                                                                                                      cGTACCTGGCGCAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGa                                                                                                                                                      <  2:2254270/150‑1 (MQ=255)    aCCGGGCGCAGTAACGGGGCCGCTACGACCAAGGTTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCCCTCTGCAAGCACAACTGGGAGGTGTGCCGAGCTGGGCATCGCTGAGCg                                                                                                                                                   <  2:847677/150‑1 (MQ=255)      cTGGCGCAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGAt                                                                                                                                                 <  1:987196/150‑1 (MQ=255)       tGGCGCAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATg                                                                                                                                                >  2:1676223/1‑150 (MQ=255)       tGGCGCAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATg                                                                                                                                                >  2:1595979/1‑150 (MQ=255)                aCCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAAc                                                                                                                                       <  2:2042440/150‑1 (MQ=255)                  cGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGc                                                                                                                                     <  2:2143809/150‑1 (MQ=255)                  cGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGc                                                                                                                                     >  2:1923204/1‑150 (MQ=255)                  cGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGc                                                                                                                                     <  1:1830446/150‑1 (MQ=255)                  cGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGc                                                                                                                                     <  2:1648323/150‑1 (MQ=255)                    aGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCa                                                                                                                                   >  2:1228981/1‑150 (MQ=255)                    aGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCa                                                                                                                                   <  2:1968595/150‑1 (MQ=255)                    aGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCa                                                                                                                                   <  1:931124/150‑1 (MQ=255)                    aGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCa                                                                                                                                   <  2:1370035/150‑1 (MQ=255)                      gCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGCACAAAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGt                                                                                                                                 >  1:1204679/1‑150 (MQ=255)                       ccGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTc                                                                                                                                <  1:1228981/150‑1 (MQ=255)                         gCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCcg                                                                                                                              <  2:2296056/150‑1 (MQ=255)                         gCTACGACCAAGGCGACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGACCAAAAACGCCAGTCcg                                                                                                                              >  1:870334/1‑150 (MQ=255)                          cTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCcgc                                                                                                                             >  1:775627/1‑150 (MQ=255)                               caccaGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCcgcacgca                                                                                                                        >  2:1577147/5‑150 (MQ=255)                                 cAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCAcc                                                                                                                      >  2:213942/1‑150 (MQ=255)                                    ggCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATg                                                                                                                   <  2:1421376/150‑1 (MQ=255)                                       tACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCCCGGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAg                                                                                                                >  1:623738/1‑150 (MQ=255)                                       tACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCAGCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAg                                                                                                                <  1:813117/150‑1 (MQ=255)                                       tACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAAAGTGGGATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAg                                                                                                                >  2:666339/1‑150 (MQ=255)                                         cGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGgt                                                                                                              <  1:1595979/150‑1 (MQ=255)                                         cGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGATGCGTATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGgt                                                                                                              <  1:1879342/150‑1 (MQ=255)                                              tCATCGCCCACTAACGCCTGAAGCGCCAACGCGAAGCCTGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATg                                                                                                         <  1:2078860/150‑1 (MQ=255)                                              tCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATg                                                                                                         >  2:1745769/1‑150 (MQ=255)                                              tCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATg                                                                                                         >  1:618005/1‑150 (MQ=255)                                                 tCGCTCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCgg                                                                                                      >  1:1541990/1‑150 (MQ=255)                                                 tCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTCAGCGATgccc                                                                                                                                             <  2:797821/111‑4 (MQ=255)                                                 tCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTCAGCGATgccc                                                                                                                                             >  1:797821/1‑108 (MQ=255)                                                     ccACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAg                                                                                                                <  2:1360910/136‑1 (MQ=255)                                                     ccACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAg                                                                                                                >  1:1360910/1‑136 (MQ=255)                                                     ccACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCAt                                                                                                  >  1:2307425/1‑150 (MQ=255)                                                      cACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATg                                                                                                 <  2:870334/150‑1 (MQ=255)                                                      cACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATg                                                                                                 <  1:1313043/150‑1 (MQ=255)                                                        cGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGt                                                                                               <  1:1923204/150‑1 (MQ=255)                                                         gAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTc                                                                                              >  2:753775/1‑150 (MQ=255)                                                          aaCGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCa                                                                                             <  1:1734522/150‑1 (MQ=255)                                                            cGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACg                                                                                           <  2:2431444/150‑1 (MQ=255)                                                             gCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGTACCATGGAGGTGTATGCGGCCATGGTCACgg                                                                                          <  1:1554430/150‑1 (MQ=255)                                                                 gAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGGGGTGTATGCGGCCATGGTCACGGCGCt                                                                                      <  1:1676223/150‑1 (MQ=255)                                                                    gcgcCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGa                                                                                   <  1:188384/150‑1 (MQ=255)                                                                           cgcgAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAg                                                                            <  2:2307425/150‑1 (MQ=255)                                                                           cgcgAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGAGGCCATGGTCACGGCGCTGGACGCCGAg                                                                            <  2:224279/150‑1 (MQ=255)                                                                            gcgAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAgg                                                                           >  1:867613/1‑150 (MQ=255)                                                                               aaGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTgg                                                                        <  1:9995/150‑1 (MQ=255)                                                                                  ccGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGgc                                                                     <  1:1061733/150‑1 (MQ=255)                                                                                  ccGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGgc                                                                     <  2:839060/150‑1 (MQ=255)                                                                                    ggCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGgcgc                                                                   >  2:569301/1‑150 (MQ=255)                                                                                       cTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCtg                                                                >  1:355582/1‑150 (MQ=255)                                                                                        tGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCtgg                                                               >  2:290288/1‑150 (MQ=255)                                                                                        accaggccgccgACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCtgg                                                               >  1:2204936/7‑150 (MQ=255)                                                                                         ggTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCtggt                                                              <  1:958406/150‑1 (MQ=255)                                                                                          gTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCtggtg                                                             <  1:1515234/150‑1 (MQ=255)                                                                                          gTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCtggtg                                                             <  1:1577147/150‑1 (MQ=255)                                                                                              ccgccgACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAac                                                         <  2:2204936/150‑1 (MQ=255)                                                                                               cgccgACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAaca                                                        >  2:2227257/1‑150 (MQ=255)                                                                                                gccgACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAacac                                                       >  2:366344/1‑150 (MQ=255)                                                                                                gccgACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATCGGCACGGCGCTGTACGCCGAGGTGGTGCGCCTGGTGGAacac                                                       >  2:145614/1‑150 (MQ=255)                                                                                                 ccgACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACAcc                                                      <  1:753775/150‑1 (MQ=255)                                                                                                    aCTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTAc                                                   <  2:775627/150‑1 (MQ=255)                                                                                                    aCTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTAc                                                   <  1:1756983/150‑1 (MQ=255)                                                                                                      ttCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGaa                                                          >  1:693240/1‑141 (MQ=255)                                                                                                      ttCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGaa                                                          <  2:693240/141‑1 (MQ=255)                                                                                                      tgccctctCCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCTTCGCTGAGCGATAAAAAAAACCGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCCCTGGACGCCGAGGGGGGGCGCCTGGGGGAACCCCTACGc                                                 <  2:1540364/146‑1 (MQ=255)                                                                                                         tctctGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAg                                              >  2:1942031/1‑150 (MQ=255)                                                                                                          ccctGCAAGCAACAATAGAAGGAGTGCAGAGATGGGCATATCGTATCAATTTACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGa                                             <  2:1842347/148‑1 (MQ=255)                                                                                                             tGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCg                                          >  1:1111379/1‑150 (MQ=255)                                                                                                              gCAAGTACAACTGGAAGGTGTGCCGAGGTGAGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCACGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCgg                                         <  2:2345872/150‑1 (MQ=255)                                                                                                               cAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGgcgc                                                                   >  2:2219675/1‑123 (MQ=255)                                                                                                               cAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGgcgc                                                                   <  1:2219675/123‑1 (MQ=255)                                                                                                                      aaCTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTgg                                 >  1:757927/1‑150 (MQ=255)                                                                                                                      aaCTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTgg                                 >  1:758084/1‑150 (MQ=255)                                                                                                                         tGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGGGCACTGCGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGaaa                              >  1:25113/1‑150 (MQ=255)                                                                                                                         tGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCtggtg                                                             >  1:1696139/1‑119 (MQ=255)                                                                                                                         tGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCtggtg                                                             <  2:1696139/119‑1 (MQ=255)                                                                                                                          ggAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGaaaa                             >  2:659476/1‑150 (MQ=255)                                                                                                                               gtgtGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGTCGAGCGGGAAAACACga                        <  2:1541990/150‑1 (MQ=255)                                                                                                                               gtgtGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACga                        <  1:1349857/150‑1 (MQ=255)                                                                                                                               gtgtGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACga                        <  1:707910/150‑1 (MQ=255)                                                                                                                                    ccGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGa                                      >  2:223388/1‑131 (MQ=255)                                                                                                                                    ccGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGa                                      <  1:223388/131‑1 (MQ=255)                                                                                                                                         cTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTGtt              <  2:160923/150‑1 (MQ=255)                                                                                                                                            ttcATCGCTCAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCgg                                         <  1:1444628/118‑1 (MQ=255)                                                                                                                                            ttcATCGCTCAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCgg                                         >  2:1444628/3‑120 (MQ=255)                                                                                                                                            ggCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGTCGTGGATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTGTTCAt           <  1:145614/150‑1 (MQ=255)                                                                                                                                               aTCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTAc                                                   >  1:1139672/1‑107 (MQ=255)                                                                                                                                               aTCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTAc                                                   <  2:1139672/107‑1 (MQ=255)                                                                                                                                               aTCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTGTTCATGTc        >  2:637860/1‑150 (MQ=255)                                                                                                                                               aTCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTGTTCATGTc        >  1:1854546/1‑150 (MQ=255)                                                                                                                                                     gAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCTAGCTGGAAAACACGATGATCCTGTTCATGTCCGACaa  <  1:366344/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CGTACCTGGCGCAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTGTTCATGTCCGACAA  >  NC_002947/3789657‑3789954

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