Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,638,667	C→A	6.1%	H221Q (CAC→CAA)	PP_3205 →	fumarylacetoacetate hydrolase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,638,667	0	C	A	6.1%	210.2 / 3.5	65	H221Q (CAC→CAA)	PP_3205	fumarylacetoacetate hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base A (2/2);  total (35/30)
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.57e-01

CCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCA  >  NC_002947/3638523‑3638815                                                                                                                                                 |                                                                                                                                                      ccTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGTGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCAc                                                                                                                                                       >  5:308707/1‑145 (MQ=255) ccTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCAc                                                                                                                                                       <  6:308707/145‑1 (MQ=255)  cTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGACTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCCCTTGCCGCTGCACGTGCa                                                                                                                                                  >  3:157301/1‑149 (MQ=255)    gCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACt                                                                                                                                                >  2:103419/1‑149 (MQ=255)     cGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCAc                                                                                                                                                 <  1:194638/147‑1 (MQ=255)        gCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGaa                                                                                                                                            >  3:28585/1‑149 (MQ=255)        gCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGaa                                                                                                                                            <  7:14055/149‑1 (MQ=255)            tGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGc                                                                                                                                        >  8:138199/1‑149 (MQ=255)                  cgcgGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTg                                                                                                                                               <  3:276751/136‑1 (MQ=255)                  cgcgGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTg                                                                                                                                               >  4:276751/1‑136 (MQ=255)                   gcgGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGt                                                                                                                                 <  4:28585/149‑1 (MQ=255)                      ggAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGt                                                                                                                                    >  5:368451/1‑143 (MQ=255)                      ggAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGt                                                                                                                                    <  6:368451/143‑1 (MQ=255)                      ggAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCg                                                                                                                              <  4:165627/149‑1 (MQ=255)                       gAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCgg                                                                                                                             <  7:138199/149‑1 (MQ=255)                                    gcttcggcttTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATggc                                                                                                                >  4:131253/1‑149 (MQ=255)                                      ttcggcttTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATggcgg                                                                                                              >  3:25454/1‑149 (MQ=255)                                           cttTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATg                                                                                                                  <  3:94518/140‑1 (MQ=255)                                           cttTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATg                                                                                                                  >  4:94518/1‑140 (MQ=255)                                              tctcCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATg                                                                                                       <  2:320631/148‑1 (MQ=255)                                               ctcCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCgg                                                                                                                             >  8:210911/1‑125 (MQ=255)                                               ctcCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCgg                                                                                                                             <  7:210911/125‑1 (MQ=255)                                                tcCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCAc                                                                                                    >  5:303968/1‑149 (MQ=255)                                                 ccAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACt                                                                                                   >  7:300481/1‑149 (MQ=255)                                                  cAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACtt                                                                                                  >  3:54873/1‑149 (MQ=255)                                                   aGGCCAAGCCCGCCAGCAGTTTCGCGCCCATTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTc                                                                                                 >  1:26389/1‑149 (MQ=255)                                                         aGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTc                                                                                          <  8:300481/149‑1 (MQ=255)                                                                 agcagTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGtt                                                                                  >  6:31700/1‑149 (MQ=255)                                                                          gcgcCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCa                                                                         >  3:3432/1‑149 (MQ=255)                                                                           cgcCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGtt                                                                                                                                <  1:54774/94‑1 (MQ=255)                                                                           cgcCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGtt                                                                                                                                >  2:54774/1‑94 (MQ=255)                                                                              ccGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATgcg                                                                     <  5:31700/149‑1 (MQ=255)                                                                                gTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATgcggc                                                                   >  8:192777/1‑149 (MQ=255)                                                                                  tGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGgcgc                                                                 >  3:331763/1‑149 (MQ=255)                                                                                     ggTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCAAGCGtcgcgtc                                                              <  4:54873/149‑4 (MQ=255)                                                                                        cACCCCGGACGAACGCGGCGAGGCCTGGCGCGACGGCCGGGTGCCCTTGCCGCTGCACGTGCACTGGAACGGCGAGTGTTTCGGCCACCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCCTGCGGCGCTgacgc                                                           <  7:876/149‑1 (MQ=255)                                                                                         aCCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCAc                                                                                                                                                 <  3:14601/63‑1 (MQ=255)                                                                                         aCCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCAc                                                                                                                                                 >  4:14601/1‑63 (MQ=255)                                                                                             cGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGgcgc                                                                 <  2:65512/138‑1 (MQ=255)                                                                                             cGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGgcgc                                                                 >  1:65512/1‑138 (MQ=255)                                                                                                 cGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCtt                                                    >  3:77588/1‑147 (MQ=255)                                                                                                 cGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGa                                                  <  7:192777/149‑1 (MQ=255)                                                                                                  gAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAg                                                 >  4:185088/1‑149 (MQ=255)                                                                                                    aCTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTg                                               >  1:97760/1‑149 (MQ=255)                                                                                                    aCTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTg                                               >  1:222832/1‑149 (MQ=255)                                                                                                        ggCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCgg                                            <  3:131253/148‑1 (MQ=255)                                                                                                           gAGGCCTGGCGCAACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCAc                                        >  8:85998/1‑149 (MQ=255)                                                                                                            gggCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACg                                       >  6:320980/2‑149 (MQ=255)                                                                                                                   gggcgACGGCCGGGGGCAATTGCCGCTGAACGGCCACGGGAACGGCGAGGGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCGTTTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTg                                 <  1:103419/147‑1 (MQ=255)                                                                                                                   gcgcGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTgg                                >  6:30307/1‑149 (MQ=255)                                                                                                                    cgcgACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTggg                               >  7:314066/1‑149 (MQ=255)                                                                                                                    cgcaACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTggg                               <  7:85998/149‑1 (MQ=255)                                                                                                                       gACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTc                            >  7:106409/1‑149 (MQ=255)                                                                                                                          ggCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCgg                         <  8:106409/149‑1 (MQ=255)                                                                                                                          ggCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCgg                         <  4:331763/149‑1 (MQ=255)                                                                                                                            ccGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCGACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCa                       >  1:134536/1‑149 (MQ=255)                                                                                                                             cGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCAc                      <  2:97760/149‑1 (MQ=255)                                                                                                                              gggTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCAcc                     <  5:30307/149‑1 (MQ=255)                                                                                                                              gggTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCAcc                     <  4:157301/149‑1 (MQ=255)                                                                                                                               agTGCACGTGCAGCGGCAAGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCg                    >  2:295782/2‑149 (MQ=255)                                                                                                                                gTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCgt                   <  3:185088/149‑1 (MQ=255)                                                                                                                                  gCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCgtgt                 >  5:16705/1‑149 (MQ=255)                                                                                                                                   cACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCgtgt                 <  6:303968/148‑1 (MQ=255)                                                                                                                                    acgtGCAGCGGCAAGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTcc               <  1:295782/146‑1 (MQ=255)                                                                                                                                      gtGCAGCGGCAAGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCaa             >  5:86334/2‑149 (MQ=255)                                                                                                                                          cGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGcc         >  4:21262/1‑149 (MQ=255)                                                                                                                                             ggCAAGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAg      <  6:86334/148‑1 (MQ=255)                                                                                                                                                 cGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCAGTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCa  <  2:134536/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                      CCTGCGGGCCCTGGCACCGCGGGAAATGAAGACCGGCTTCGGCTTTCTCCAGGCCAAGCCCGCCAGCAGTTTCGCGCCCGTTGCGGTCACCCCGGACGAACTCGGCGAGGCCTGGCGCGACGGCCGGGTGCACTTGCCGCTGCACGTGCACTGGAACGGCGAGTGGTTCGGCCATCCCCATGGCGGCGCCATGCACTTCG‑GTTTCCACCAGTTGATTGCCCATGCGGCGCTGACGCGACGCTTGAGTGCGGGCACGCTACTGGGCTCCGGCACCGTGTCCAACGCCGAGCGCA  >  NC_002947/3638523‑3638815

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