Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	239,053	C→T	5.7%	P45L (CCA→CTA) ‡	hemD →	uroporphyrinogen‑III synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	239,053	0	C	T	5.7%	273.9 / 3.4	90	P45L (CCA→CTA) ‡	hemD	uroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand):  ref base C (39/45);  new base T (3/2);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

GGCCGGTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGC  >  NC_002947/238908‑239200                                                                                                                                                  |                                                                                                                                                       ggCCGGTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCcaca                                                                                                                                                     <  4:390954/149‑1 (MQ=255)     ggTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCgg                                                                                                                                                 <  8:210588/149‑1 (MQ=255)     ggTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCgg                                                                                                                                                 >  8:131335/1‑149 (MQ=255)     ggTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCgg                                                                                                                                                 >  7:210588/1‑149 (MQ=255)         aCCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGt                                                                                                                                             <  7:107354/149‑1 (MQ=255)          cccGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑c                                                                                                                                           <  2:391011/149‑1 (MQ=255)            cGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGc                                                                                                                                         >  8:280784/1‑149 (MQ=255)                  ccAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAg                                                                                                                                   >  7:112226/1‑149 (MQ=255)                  ccAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGCCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCTTCCAGCGACAGC‑CGCTgtggg                                                                                                                                   >  3:242999/1‑145 (MQ=255)                        gcgcCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGc                                                                                                                             <  6:48903/149‑1 (MQ=255)                          gcCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGCCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCTTCCAGCGACAGC‑CGCTGTGGGGCTTGC‑ac                                                                                                                          <  4:242999/149‑2 (MQ=38)                             tGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGAc                                                                                                                       >  8:32157/1‑149 (MQ=255)                               ttGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACtt                                                                                                                     <  3:210434/149‑1 (MQ=255)                                   tGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGc                                                                                                                             <  5:103467/138‑1 (MQ=255)                                   tGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGc                                                                                                                             >  6:103467/1‑138 (MQ=255)                                    gACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAgg                                                                                                                >  8:289628/1‑149 (MQ=255)                                       ccGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGa                                                                                                             <  7:131335/149‑1 (MQ=255)                                       ccGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGa                                                                                                             <  2:263756/149‑1 (MQ=255)                                         ggCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGAt                                                                                                            >  1:110023/1‑148 (MQ=255)                                               ccGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATtgtg                                                                                                     <  7:289628/149‑1 (MQ=255)                                                cGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTgg                                                                                                    >  3:114211/1‑149 (MQ=255)                                                 gAGGACTGCGCGGGGCTGGCGCATAGCCTGGCGGCAGCGGGAGTGGGCAGCTGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGt                                                                                                   <  2:22245/149‑1 (MQ=255)                                                 gAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGt                                                                                                   >  1:155196/1‑149 (MQ=255)                                                 gAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGt                                                                                                   <  3:194770/149‑1 (MQ=255)                                                 gAGGACTGCGCGGCGCTGCCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCAAGGCGATCATTGTGGt                                                                                                   >  8:375163/1‑149 (MQ=255)                                                     aCTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGc                                                                                               <  7:379555/149‑1 (MQ=255)                                                      cTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGc                                                                                               >  3:76151/1‑148 (MQ=255)                                                       tGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACt                                                                                                                      >  2:164605/1‑124 (MQ=255)                                                       tGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACt                                                                                                                      <  1:164605/124‑1 (MQ=255)                                                                   gcgcAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCg                                                                                 >  4:269118/1‑149 (MQ=255)                                                                   gcgcAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCg                                                                                 >  2:1441/1‑149 (MQ=255)                                                                     gcAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCGACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGc                                                                               >  1:128323/1‑149 (MQ=255)                                                                     gcAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGc                                                                               >  6:405342/1‑149 (MQ=255)                                                                     gcAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGc                                                                               >  4:221752/1‑149 (MQ=255)                                                                      cAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCt                                                                              <  7:32157/149‑1 (MQ=255)                                                                      cAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCt                                                                              >  5:186157/1‑149 (MQ=255)                                                                        aaGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTTCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGt                                                                            <  3:16022/149‑1 (MQ=255)                                                                         aGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGACAGCAAGCCGGCCGCCCGGCTGtt                                                                           <  6:251559/149‑1 (MQ=255)                                                                          gCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGTACGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGtt                                                                           <  1:121380/149‑1 (MQ=255)                                                                          gCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTg                                                                          <  2:155196/149‑1 (MQ=255)                                                                              ggcggcAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGa                                                                                                                                     >  6:183445/1‑87 (MQ=255)                                                                              ggcggcAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGa                                                                                                                                     <  5:183445/87‑1 (MQ=255)                                                                               gcggcAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑Ga                                                                    <  5:271807/149‑1 (MQ=255)                                                                               gcggcAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑Ga                                                                    <  4:114211/149‑1 (MQ=255)                                                                                  gcAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCa                                                                 >  5:185278/1‑149 (MQ=255)                                                                                       gggAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTgg                                                            <  4:76151/149‑1 (MQ=255)                                                                                         gAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGcc                                                          >  4:80375/1‑149 (MQ=255)                                                                                         gAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCAATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGcc                                                          >  6:198805/1‑149 (MQ=255)                                                                                              ggCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTg                                                                          <  1:278536/129‑1 (MQ=255)                                                                                              ggCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTg                                                                          >  2:278536/1‑129 (MQ=255)                                                                                              ggCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCgaga                                                      <  3:221752/148‑1 (MQ=255)                                                                                              ggCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGAc                                                     <  6:185278/149‑1 (MQ=255)                                                                                              ggCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGAc                                                     >  6:61325/1‑149 (MQ=255)                                                                                                  gcagTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGc                                                 <  4:208996/149‑1 (MQ=255)                                                                                                      ttGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATc                                              >  3:204500/1‑148 (MQ=255)                                                                                                         cctgccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGc                                          >  5:59643/1‑149 (MQ=255)                                                                                                          ctgccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGAc                                                     <  6:28654/137‑1 (MQ=255)                                                                                                          ctgccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGcc                                         <  3:145776/149‑1 (MQ=255)                                                                                                          ctgccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGcc                                         >  1:395418/1‑149 (MQ=255)                                                                                                            gccCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCAAAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAACCGGCCGCCCGGCTGTTGCTT‑GAGCGTCTGGCCGAGACCAGCATCCAGCCTc                                       >  7:109103/1‑149 (MQ=255)                                                                                                             ccccTGCTGGAGAATGAACCCATTACCGAGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTc                                       <  8:173995/148‑1 (MQ=255)                                                                                                               ccTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTTAGAGCATCTGGCCGAGACCAGCATCCAGCCTCCg                                     <  1:259152/149‑1 (MQ=255)                                                                                                                    tgGAGATTGAAGCCATTACCGTGCTAGCGGCTCAGCGGGTGT‑CGCTGGAAGGCTTGC‑ATGACTTCGAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGa                               <  4:204500/149‑1 (MQ=255)                                                                                                                     ggAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAgg                                                                                                                <  1:109340/68‑1 (MQ=255)                                                                                                                     ggAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAgg                                                                                                                >  2:109340/1‑68 (MQ=255)                                                                                                                     ggAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAg                              >  8:185372/1‑149 (MQ=255)                                                                                                                      gagaTTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGa                               <  5:380588/147‑1 (MQ=255)                                                                                                                          ttGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTgg                         <  1:1441/149‑1 (MQ=255)                                                                                                                           tGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCTATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTgg                         <  5:61325/148‑1 (MQ=255)                                                                                                                           tGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGt                        >  6:320611/1‑149 (MQ=255)                                                                                                                              aCCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCa                     <  7:185372/149‑1 (MQ=255)                                                                                                                              aCCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCa                     >  7:339987/1‑149 (MQ=255)                                                                                                                               cccATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCAc                    <  1:187715/149‑1 (MQ=255)                                                                                                                                ccATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCAcc                   >  6:76931/1‑149 (MQ=255)                                                                                                                                  aTTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGc                                 <  7:276089/133‑1 (MQ=255)                                                                                                                                  aTTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGc                                 >  8:276089/1‑133 (MQ=255)                                                                                                                                  aTTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGt                 >  6:32285/1‑149 (MQ=255)                                                                                                                                  aTTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGACTGGTTc                      >  3:49335/1‑144 (MQ=255)                                                                                                                                  aTTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGACTGGTTc                      <  4:49335/144‑1 (MQ=255)                                                                                                                                   ttACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTg                <  3:237394/149‑1 (MQ=255)                                                                                                                                     gtcatGCAAGCCTTCCAGCGACAGC‑CGCTGTGGGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGAc                                                     >  5:93235/5‑110 (MQ=21)                                                                                                                                     gtcatGCAAGCCTTCCAGCGACAGC‑CGCTGTGGGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGAc                                                     <  6:93235/106‑1 (MQ=21)                                                                                                                                     aCCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTgg               <  6:59643/148‑1 (MQ=255)                                                                                                                                      tcatGCAAGCCTTCCAGCGACAGC‑CGCTGTGGGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGt             >  8:353268/4‑149 (MQ=37)                                                                                                                                      ccGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGt             >  2:58020/1‑149 (MQ=255)                                                                                                                                        gTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGa           <  5:320611/149‑1 (MQ=255)                                                                                                                                        gTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGa           <  2:128323/149‑1 (MQ=255)                                                                                                                                        gTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGa           <  2:395418/149‑1 (MQ=255)                                                                                                                                         tGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGCGATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGt             >  1:271607/1‑147 (MQ=255)                                                                                                                                            aaGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGc        >  2:99763/1‑148 (MQ=255)                                                                                                                                              gCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAAc     <  6:186157/149‑1 (MQ=255)                                                                                                                                                cccACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCg   >  7:2175/1‑149 (MQ=255)                                                                                                                                                 ccACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCtt                                                                       <  4:360745/81‑1 (MQ=255)                                                                                                                                                 ccACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCtt                                                                       >  3:360745/1‑81 (MQ=255)                                                                                                                                                 ccACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCCGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGc  >  8:84352/1‑149 (MQ=255)                                                                                                                                                 ccACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGc  >  4:209527/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                       GGCCGGTCACCCGTGAGCCAATGGCGCCTGTTGCTGACCCGGCCTGCCGAGGACTGCGCGGCGCTGGCGCAAAGCCTGGCGGCAGCGGGAGTGGGCAGCAGTTGCCTGCCCCTGCTGGAGATTGAACCCATTACCGTGCAAGCCCCACAGCGGCTGT‑CGCTGGAAGGCTTGC‑ATGACTTCCAGGCGATCATTGTGGTCAGCAAGCCGGCCGCCCGGCTGTTGCTT‑GAGCATCTGGCCGAGACCAGCATCCAGCCTCCGCAGCGAGGCTGGTTCACCGTGGGTGAGGCAACCGC  >  NC_002947/238908‑239200

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