Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	214,309	G→A	5.5%	D12N (GAC→AAC)	tauD_1 →	Alpha‑ketoglutarate‑dependent taurinedioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	214,309	0	G	A	5.5%	260.2 / 3.2	77	D12N (GAC→AAC)	tauD_1	Alpha‑ketoglutarate‑dependent taurinedioxygenase
Reads supporting (aligned to +/- strand):  ref base G (30/42);  new base A (2/2);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

GATGCGCCGCACGACTTGAACCCACACGGGAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCTCGACCACC  >  ACB143/214183‑214444                                                                                                                               |                                                                                                                                        gATGCGCCGCACGACTTGAACCCACACGGGAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGTGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGcc                                                                                                                                 >  8:351782/1‑135 (MQ=255) gATGCGCCGCACGACTTGAACCCACACGGGAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGccc                                                                                                                                >  1:120311/1‑136 (MQ=255) gATGCGCCGCACGACTTGAACCCACACGGGAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGcc                                                                                                                                 >  7:109548/1‑135 (MQ=255) gATGCGCCGCACGACTTGAACCCACACGGGAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGcc                                                                                                                                 >  7:162510/1‑135 (MQ=255) gATGCGCCGCACGACTTGAACCCACACGGGAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGcc                                                                                                                                 >  8:188449/1‑135 (MQ=255)     cgccgcACGACTTGAACCCACACGGGAAATCCCCGTGCGTTTGTGATCGGTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGcc                                                                                                                            <  7:183839/136‑1 (MQ=255)     cgccgcACGACTTGAACCCACACGGGAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGcc                                                                                                                            <  7:222336/136‑1 (MQ=255)     cgccgcACGACTTGAACCCACACGGGAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGcc                                                                                                                            <  4:127104/136‑1 (MQ=255)                   aaCCCACACGGGAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGaa                                                                                                               <  1:64750/135‑1 (MQ=255)                   aaCCCACACGGGAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGaa                                                                                                               <  7:294925/135‑1 (MQ=255)                     cccACACGGGAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGtt                                                                                                            <  5:180180/136‑1 (MQ=255)                     cccACACGGGAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGtt                                                                                                            <  8:162510/136‑1 (MQ=255)                           cGGGAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGcccc                                                                                                      >  2:6671/1‑136 (MQ=255)                            gggAAATCCCCGTGCGTTTCTGATCGTTTCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTCGAAGTTCGCCCCt                                                                                                     >  7:77929/1‑136 (MQ=255)                            gggAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTCAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCACGCCGCAACTAATTGAAGTTCGCCCCt                                                                                                     >  4:233574/1‑136 (MQ=255)                            gggAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCAGTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGCACGCCTTGTCGTCCATCGACCGCG‑CCAGCCGCAAGTTTTTGAAGTTGCCCCCtt                                                                                                    >  6:55225/1‑136 (MQ=255)                                        tGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCg                                                                                         >  1:372412/1‑136 (MQ=255)                                         gagTTTCTGATCGTTGCGCGCCATTTAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCATGAACGACATGTAGTCCATCGCCAGCGCCCAGCCGCAAGTTTTTTAAGTTAGCCCCTTATAGGGTTCCGt                                                                                        >  5:55765/3‑136 (MQ=255)                                           gTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCAAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCg                                                                                      >  4:149076/1‑136 (MQ=255)                                              tCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGgc                                                                                    <  1:231104/135‑1 (MQ=255)                                              tCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGgc                                                                                    <  6:158921/135‑1 (MQ=255)                                                     gTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCgaga                                                                             <  1:304309/135‑1 (MQ=255)                                                     gTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCgaga                                                                             <  5:135632/135‑1 (MQ=255)                                                     gTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCgaga                                                                             <  4:282939/135‑1 (MQ=255)                                                     gTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCgaga                                                                             <  2:156163/135‑1 (MQ=255)                                                     gTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCgaga                                                                             <  3:96469/135‑1 (MQ=255)                                                     gTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCgaga                                                                             <  7:133545/135‑1 (MQ=255)                                                     gTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCgaga                                                                             <  4:168550/135‑1 (MQ=255)                                                     gTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGAt                                                                            <  3:126090/136‑1 (MQ=255)                                                          gcgcCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGttt                                                                                                                    <  8:183905/91‑1 (MQ=255)                                                          gcgcCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGttt                                                                                                                    >  7:183905/1‑91 (MQ=255)                                                          gcgcCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCATACTCGGGTGCCGTCGGCGCCGAGATAtttg                                                                       >  8:373864/1‑134 (MQ=255)                                                          gcgcCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCAACAGCGCCCAGCCGCAAGTTTTTGAAGTTCg                                                                                                          <  6:309988/101‑1 (MQ=255)                                                          gcgcCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCAACAGCGCCCAGCCGCAAGTTTTTGAAGTTCg                                                                                                          <  6:315726/101‑1 (MQ=255)                                                          gcgcCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCAACAGCGCCCAGCCGCAAGTTTTTGAAGTTCg                                                                                                          >  5:309988/1‑101 (MQ=255)                                                          gcgcCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCAACAGCGCCCAGCCGCAAGTTTTTGAAGTTCg                                                                                                          >  5:315726/1‑101 (MQ=255)                                                               aTTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTg                                                                  <  7:319255/136‑1 (MQ=255)                                                                 tGAGGTCGCTCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTg                                                                  <  7:234317/134‑1 (MQ=255)                                                                        gcgcATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGc                                                          >  3:177166/1‑135 (MQ=255)                                                                           cATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGCCCTGGCCaac                                                      >  4:35066/1‑135 (MQ=255)                                                                           cATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCaac                                                      >  4:29684/1‑135 (MQ=255)                                                                           cATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCaaa                                                      >  7:139432/1‑136 (MQ=255)                                                                           cATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCaaa                                                      >  2:51444/1‑136 (MQ=255)                                                                           cATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCaaa                                                      >  6:361014/1‑136 (MQ=255)                                                                                   gggAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGGCCTGGCCAAACCGGTCAg                                              >  2:202527/1‑136 (MQ=255)                                                                                     gactGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAgcg                                            <  4:10516/133‑1 (MQ=255)                                                                                     gAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAgcg                                            >  7:179390/1‑136 (MQ=255)                                                                                     gAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAgcg                                            <  2:42062/136‑1 (MQ=255)                                                                                           aaaaTGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAg                                       >  3:118626/1‑135 (MQ=255)                                                                                            aaaTGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGa                                     <  7:330591/136‑1 (MQ=255)                                                                                            aaaTGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGa                                     <  1:12608/136‑1 (MQ=255)                                                                                            aaaTGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGa                                     <  5:62124/136‑1 (MQ=255)                                                                                             aaTGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGAt                                    >  2:155266/1‑136 (MQ=255)                                                                                             aaTGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGAt                                    <  3:326599/136‑1 (MQ=255)                                                                                              aTGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGAtt                                   <  6:6824/136‑1 (MQ=255)                                                                                               tGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGAttt                                  <  2:330652/136‑1 (MQ=255)                                                                                                 tCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTc                                 <  8:139432/135‑1 (MQ=255)                                                                                                     ccgccTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCg                            <  2:15092/136‑1 (MQ=255)                                                                                                         cTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATc                        <  8:91677/136‑1 (MQ=255)                                                                                                         cTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATc                        <  2:240727/136‑1 (MQ=255)                                                                                                              aCGCCTTTTCGTCCATCGACAGCGCCCAGCCGGAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATc                    <  6:45846/135‑1 (MQ=255)                                                                                                                gCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATcgcg                 <  8:124148/136‑1 (MQ=255)                                                                                                                  cTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGcc               <  5:92397/136‑1 (MQ=255)                                                                                                                     gtcgtcCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCAc            >  1:246021/1‑136 (MQ=255)                                                                                                                     gtcgtcCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCAc            >  4:359716/1‑136 (MQ=255)                                                                                                                     gtcgtcCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTAACCATCGCGCCCAc            >  6:329293/1‑136 (MQ=255)                                                                                                                     gtcgtcCATCCACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCAc            >  6:8803/1‑136 (MQ=255)                                                                                                                      tcgtcCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGtt                                                                                                            <  4:175939/39‑1 (MQ=255)                                                                                                                      tcgtcCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGtt                                                                                                            >  3:175939/1‑39 (MQ=255)                                                                                                                        gtcCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCt          <  8:133008/135‑1 (MQ=255)                                                                                                                        gtcCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCt          <  5:357199/135‑1 (MQ=255)                                                                                                                            aTCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCTCGacc     <  5:36150/136‑1 (MQ=255)                                                                                                                            aTCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCTCGacc     <  3:359716/136‑1 (MQ=255)                                                                                                                            aTCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCTCGacc     <  6:68899/136‑1 (MQ=255)                                                                                                                            aTCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCTCGacc     <  5:13422/136‑1 (MQ=255)                                                                                                                              cGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCTCGaccac   >  3:192973/1‑136 (MQ=255)                                                                                                                              cGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCTCGaccac   >  7:124239/1‑136 (MQ=255)                                                                                                                               gACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCTCGaccacc  <  3:52008/136‑1 (MQ=255)                                                                                                                               gACAGCGCCCAGCCGCAAGTTTTCGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCTCGaccacc  <  8:77929/136‑1 (MQ=255)                                                                                                                               |                                                                                                                                        GATGCGCCGCACGACTTGAACCCACACGGGAAATCCCCGTGCGTTTCTGATCGTTGCGCGCCATTGAGGTCGCGCATTGCGTGGGAGTGAAAAATGTCAGCCGCCTTGAACGCCTTGTCGTCCATCGACAGCGCCCAGCCGCAAGTTTTTGAAGTTCGCCCCTTCTCGGGTGCCGTCGGCGCCGAGATCTTAGGCCTGGACCTGGCCAAACCGGTCAGCGCCGAGGATTTCACCCGTATCCATCGCGCCCACCTCGACCACC  >  ACB143/214183‑214444

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