Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,789,762	G→C	11.6%	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,762	0	G	C	11.6%	217.5 / 16.7	76	L289L (CTG→CTC)	atsA	arylsulfatase
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base C (5/4);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTACCTGGCGCAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTGTTCATGTCC  >  NC_002947/3789615‑3789906                                                                                                                                                     |                                                                                                                                                 gTACCTGGCGCAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGa                                                                                                                                                 <  8:41573/149‑1 (MQ=255) gTACCTGGCGCAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGa                                                                                                                                                 <  5:233959/149‑1 (MQ=255)  tACCTGGCGCAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAg                                                                                                                                                <  4:257682/149‑1 (MQ=255)    ccTGGCGCAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGc                                                                                                                                               <  1:46629/148‑1 (MQ=255)     cTGGCGCAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGa                                                                                                                                             >  1:167924/1‑149 (MQ=255)        gcgcAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATg                                                                                                                                           >  6:361808/1‑148 (MQ=255)         cgcAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTCAGCGATgcc                                                                                                                                         <  8:335746/149‑3 (MQ=255)          gcAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTCAGCGATgccc                                                                                                                                        <  8:117859/149‑4 (MQ=255)              tACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACaaa                                                                                                                                     >  5:282426/1‑148 (MQ=255)                ccGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAAc                                                                                                                                  <  4:244449/149‑1 (MQ=255)                  gAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGc                                                                                                                                <  2:167924/149‑1 (MQ=255)                      ccGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCGCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGt                                                                                                                            <  2:103629/149‑1 (MQ=255)                       cGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTc                                                                                                                           >  3:260728/1‑149 (MQ=255)                       cGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTc                                                                                                                           <  5:361808/149‑1 (MQ=255)                        gCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAg                                                                                                                                                <  1:78424/127‑1 (MQ=255)                        gCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAg                                                                                                                                                >  2:78424/1‑127 (MQ=255)                        gCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTcc                                                                                                                          <  8:136415/149‑1 (MQ=255)                        atcagaaCCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTcc                                                                                                                          <  4:346593/143‑1 (MQ=255)                         cTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCcg                                                                                                                         <  8:125482/149‑1 (MQ=255)                          tACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCcg                                                                                                                         >  4:182817/1‑148 (MQ=255)                          tACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTCAGCGATgcc                                                                                                                                         >  3:376148/1‑130 (MQ=255)                          tACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTCAGCGATgcc                                                                                                                                         <  4:376148/132‑3 (MQ=255)                           aCGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCcgca                                                                                                                       >  6:222307/1‑149 (MQ=255)                              aCCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCcgcacg                                                                                                                     >  1:378684/1‑148 (MQ=255)                                  aGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCa                                                                                                                <  7:67930/149‑1 (MQ=255)                                     cTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATgg                                                                                                             >  8:217318/1‑149 (MQ=255)                                       aCGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAAGGCCGGCCTGGTGCCCGCCGACTTCTCACTCCAAGCACACCGGGAAGGTGTGCCGTGCTGGGCAACGCTGAGCGAGGAACAAAAACGCCAGTCCGCACCCACCATggc                                                                                                            >  1:270119/1‑148 (MQ=255)                                        cGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAg                                                                                                           <  4:46194/148‑1 (MQ=255)                                                  gCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCgg                                                                                                 <  4:281216/148‑1 (MQ=255)                                                   cccACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGc                                                                                                >  1:385097/1‑148 (MQ=255)                                                      aCGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATg                                                                                            <  3:182817/149‑1 (MQ=255)                                                       cGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATgg                                                                                           <  6:282426/149‑1 (MQ=255)                                                          aCGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCa                                                                                        >  7:103912/1‑149 (MQ=255)                                                                 aaGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCt                                                                                 <  4:225330/149‑1 (MQ=255)                                                                     gcCAACGCGAA‑GCCGGCCTGGTGTCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGAc                                                                             >  4:380833/1‑149 (MQ=255)                                                                      ccAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACg                                                                            <  4:260728/149‑1 (MQ=255)                                                                       cAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGc                                                                           >  3:371759/1‑149 (MQ=255)                                                                        aaCGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGcc                                                                          >  2:173035/1‑149 (MQ=255)                                                                          cgcgAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGa                                                                        <  8:103912/149‑1 (MQ=255)                                                                           gcgAT‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGAt                                                                                                                                            <  7:281765/80‑1 (MQ=255)                                                                           gcgAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGAt                                                                                                                                            >  8:281765/1‑80 (MQ=255)                                                                               a‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTgg                                                                   <  2:385097/149‑1 (MQ=255)                                                                                 gCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTggg                                                                  <  4:113576/149‑1 (MQ=255)                                                                                     gCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGgcgc                                                              <  2:267284/149‑1 (MQ=255)                                                                                          gTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCtggt                                                         >  7:347168/1‑149 (MQ=255)                                                                                          gTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCtg                                                           <  7:217318/147‑1 (MQ=255)                                                                                            gCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCtggtgg                                                       <  5:222307/149‑1 (MQ=255)                                                                                             cccGCCGACTTCTCTCTGCAATCACAACATGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATTAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCtggttg                                                       <  8:347168/148‑3 (MQ=255)                                                                                               cgccgACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGaa                                                     >  3:34734/1‑148 (MQ=255)                                                                                                    aCTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGaa                                                     <  1:132631/143‑1 (MQ=255)                                                                                                    aCTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGaa                                                     >  2:132631/1‑143 (MQ=255)                                                                                                         tctctGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGc                                            >  1:177917/1‑147 (MQ=255)                                                                                                           tctGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGa                                        >  2:334007/1‑149 (MQ=255)                                                                                                             tGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGgtgt                                                                                                       <  5:298265/84‑1 (MQ=255)                                                                                                             tGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGgtgt                                                                                                       >  6:298265/1‑84 (MQ=255)                                                                                                             tGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAg                                         <  1:243241/146‑1 (MQ=255)                                                                                                             tGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGAcc                                      >  7:108701/1‑149 (MQ=255)                                                                                                               cAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCgg                                    >  2:44315/1‑149 (MQ=255)                                                                                                                aaGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGc                                   <  4:34734/149‑1 (MQ=255)                                                                                                                aaGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGc                                   >  4:119064/1‑149 (MQ=255)                                                                                                                    acaACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGc                               >  1:383528/1‑149 (MQ=255)                                                                                                                    acaACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGc                               >  4:142532/1‑149 (MQ=255)                                                                                                                    acaACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGACGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGc                               >  6:174115/1‑149 (MQ=255)                                                                                                                     caACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCt                              <  3:380833/149‑1 (MQ=255)                                                                                                                                        gCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATgg                                                                                           <  6:43932/69‑1 (MQ=255)                                                                                                                                        gCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATgg                                                                                           >  5:43932/1‑69 (MQ=255)                                                                                                                                         ttgttcATCGCTCAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTGt          >  6:174780/6‑149 (MQ=255)                                                                                                                                         cTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTg           >  1:26025/1‑148 (MQ=255)                                                                                                                                         cTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTg           <  2:26025/148‑1 (MQ=255)                                                                                                                                           gttcATCGCTCAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGaaa                         >  1:33363/4‑132 (MQ=255)                                                                                                                                           gttcATCGCTCAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGaaa                         <  2:33363/129‑1 (MQ=255)                                                                                                                                           gttcATCGCTCAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTGtt         >  7:60272/4‑148 (MQ=255)                                                                                                                                           gttcATCGCTCAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTGtt         >  5:191294/4‑148 (MQ=255)                                                                                                                                            ggCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTGTTc        <  1:173035/148‑1 (MQ=255)                                                                                                                                              cATCGCTGAGCGATGAACAAAAACGACACTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTCCGCAAGACCGTCGAGCTGGAAAACACGATCATCCTGTTCATg     >  4:314194/1‑149 (MQ=255)                                                                                                                                                  gCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTGTTCATGTcc  >  4:96888/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GTACCTGGCGCAGTACCGAGGCCGCTACGACCAAGGCTACGAAGTCATCGCCCACGAACGCCTGAAGCGCCAACGCGAA‑GCCGGCCTGGTGCCCGCCGACTTCTCTCTGCAAGCACAACTGGAAGGTGTGCCGAGCTGGGCATCGCTGAGCGATGAACAAAAACGCCAGTCCGCACGCACCATGGAGGTGTATGCGGCCATGGTCACGGCGCTGGACGCCGAGGTGGGGCGCCTGGTGGAACACCTACGCAAGACCGGCGAGCTGGAAAACACGATGATCCTGTTCATGTCC  >  NC_002947/3789615‑3789906

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