Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	79,227	T→A	5.9%	D175V (GAC→GTC)	betC ←	Choline‑sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	79,227	0	T	A	5.9%	242.1 / 2.6	77	D175V (GAC→GTC)	betC	Choline‑sulfatase
Reads supporting (aligned to +/- strand):  ref base T (40/31);  new base A (3/2);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGT  >  ACB143/79097‑79362                                                                                                                                   |                                                                                                                                        ggCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTAc                                                                                                                                    <  1:314864/136‑1 (MQ=255)      ggggCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTa                                                                                                                               >  7:364466/1‑136 (MQ=255)        ggCATGGTGATATCCACAGCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACTGTCAGGCAAAATGGCCGGCCATCGTTTTCGAGCACATGGTCGTACAGGTACt                                                                                                                             <  3:68423/136‑1 (MQ=255)        ggCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACt                                                                                                                             <  4:164960/136‑1 (MQ=255)           aTGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGgc                                                                                                                          >  1:89564/1‑136 (MQ=255)           aTGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGgc                                                                                                                          >  7:81931/1‑136 (MQ=255)           aTGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGCTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGgc                                                                                                                          >  7:323253/1‑136 (MQ=255)                  tatCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCtt                                                                                                                   <  4:337090/136‑1 (MQ=255)                     ccACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGaa                                                                                                                >  1:123954/1‑136 (MQ=255)                     ccACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGaa                                                                                                                >  2:71816/1‑136 (MQ=255)                     ccACACCTTCGTAGCGGTCCCAGTAGAGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGaa                                                                                                                >  8:46531/1‑136 (MQ=255)                       acacCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAaca                                                                                                              >  7:119012/1‑136 (MQ=255)                            ttCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACaccacc                                                                                                         >  7:318981/1‑136 (MQ=255)                            ttCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACaccacc                                                                                                         >  1:40358/1‑136 (MQ=255)                                tAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCAcctcct                                                                                                     >  8:296551/1‑136 (MQ=255)                                 aGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTCCAGGTACTGGCGCGCCTTGAACACCACctcctc                                                                                                    >  7:356304/1‑136 (MQ=255)                                 aGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACctcctc                                                                                                    >  5:281620/1‑136 (MQ=255)                                 aGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACctcctc                                                                                                    >  3:385784/1‑136 (MQ=255)                                  gCGGTCCCAGTAGCGTTTGGCAATGGTGGAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCtcg                                                                                                   <  4:80380/136‑1 (MQ=255)                                    ggTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGACGTACAGGTACTGGCGCGCCTTGAACACCACCTCCtcgtc                                                                                                 >  3:65729/1‑136 (MQ=255)                                    ggTCCCAGTAGCGTTTGGCAATGGTGTAGGCGTCATGCGGGTGGGGCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGACGTACAGGTACTGGCGCGCCTTGAACACCACCTCCtcgtc                                                                                                 >  3:65735/1‑136 (MQ=255)                                        ccAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCGTCGTCGata                                                                                             >  8:365404/1‑134 (MQ=255)                                        ccAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGaaa                                                                                             >  4:338927/1‑136 (MQ=255)                                        ccAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGaa                                                                                              <  3:136461/135‑1 (MQ=255)                                        ccAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGaa                                                                                              >  5:209526/1‑135 (MQ=255)                                        ccAGTAGCGTTTGGCAATGGTGTAGGGGTAATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACATGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGaaa                                                                                             >  8:335653/1‑136 (MQ=255)                                          aGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTACTCGTCGAAATc                                                                                           <  3:63790/136‑1 (MQ=255)                                                  ttGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGAt                                                                                   >  4:86214/1‑136 (MQ=255)                                                   tGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTAGTCGAAATCCAGCTGAtt                                                                                  >  4:1484/1‑136 (MQ=255)                                                       aaTGGTGTATGGGTCATGCTGGTGGGTCATCGAAAAGGTCAGGCAAAATGGCCGGCCATAGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGAATGGTg                                                                              >  8:64394/1‑136 (MQ=255)                                                       aaTGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTAGAGGTACTGGCGCGCCTTGAACACCACCTCATCGTCGAAATCCAGCTGATTGGTg                                                                              >  2:55562/1‑136 (MQ=255)                                                         tGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTgcg                                                                            >  1:135064/1‑136 (MQ=255)                                                             gtAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGcgcacg                                                                        >  5:64750/1‑136 (MQ=255)                                                              tAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGcgcacgc                                                                       <  4:206272/136‑1 (MQ=255)                                                                 gggTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGcacgcacg                                                                    <  6:125016/136‑1 (MQ=255)                                                                 gggTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGcacgcacg                                                                    <  7:116588/136‑1 (MQ=255)                                                                   gTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGa                                                                  <  3:170086/136‑1 (MQ=255)                                                                        gCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGATTACCCGCTGATTGGTGCGCCCGCCCGGAcccac                                                             >  6:255823/1‑134 (MQ=255)                                                                        gCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGc                                                             >  6:3095/1‑136 (MQ=255)                                                                        gCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCATCGTCTAAATCCAGCTGATTGGTGCTCACGCACGGACCCGc                                                             >  6:251264/1‑136 (MQ=255)                                                                         cgtgtgggtCATCGAAACTGTCAGGCAAAATGGACGGCCATCGTTTTCGCGCACATTTACGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCACCTGATTGGTGCGCACGCACGGACCCGcc                                                            >  7:33658/4‑136 (MQ=255)                                                                         cGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGaa                                                                                              <  3:280873/102‑1 (MQ=255)                                                                         cGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGaa                                                                                              >  4:280873/1‑102 (MQ=255)                                                                         cGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGa                                                                  >  2:331253/1‑130 (MQ=255)                                                                         cGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGa                                                                  <  1:331253/130‑1 (MQ=255)                                                                         cGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGcc                                                            >  5:385029/1‑136 (MQ=255)                                                                         cGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGcc                                                            >  2:372587/1‑136 (MQ=255)                                                                         cGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGcc                                                            >  2:155092/1‑136 (MQ=255)                                                                         cGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTCATTGGTGCGCACGCACGGACCCGcc                                                            >  4:61742/1‑136 (MQ=255)                                                                         cGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACAACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGcc                                                            >  7:106625/1‑136 (MQ=255)                                                                         cGGGTGGGTCATCGAAAATTTCATTCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGcc                                                            >  3:162073/1‑136 (MQ=255)                                                                          gggtgtgtCATCGAAACGGTCAGGCATAATGGCCGGCCATCGTTTTCGAGCACATGGGCTTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCTACCTGATTGGTGCGCACGCACGGACCCGCCt                                                           <  5:371042/136‑1 (MQ=255)                                                                          gggtgggtCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCt                                                           <  4:220844/136‑1 (MQ=255)                                                                             tgggtCATCGAAACGGTCAGGCAAAATGGGCGGCGATCGTTTTCGCGCACTTGGTCGTACTGGTACTGGGGCGGCTTGAACACCAGCTGCTCGTCGAAATGCAGCTGATTGGTGCGGACGCTCGGAGCCGCCTgca                                                        <  5:300591/136‑1 (MQ=255)                                                                                gtCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTgcagca                                                     <  6:209526/136‑1 (MQ=255)                                                                                gtCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTgcagca                                                     >  6:383976/1‑136 (MQ=255)                                                                                 tCATCGAAACGGTCAGGCAAAATGGCCGGCCATCTTTTTCTCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACCCCACCTCCTCGTCGAAATACAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCAc                                                    >  7:339780/1‑136 (MQ=255)                                                                                   aTCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCg                                                  <  6:201661/136‑1 (MQ=255)                                                                                         aCGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACa                                            >  5:182685/1‑136 (MQ=255)                                                                                         aCGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACa                                            <  8:119012/136‑1 (MQ=255)                                                                                          cGGTCAGGCAAAATGGCCGGCCATCGGTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACAt                                           <  7:73010/136‑1 (MQ=255)                                                                                                ggCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTtgtg                                     <  7:17647/136‑1 (MQ=255)                                                                                                ggCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTtgtg                                     <  6:193099/136‑1 (MQ=255)                                                                                                ggCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTtgtg                                     <  1:290642/136‑1 (MQ=255)                                                                                                 gCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTgg                                    <  6:64750/136‑1 (MQ=255)                                                                                                          ccggccATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCTCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCtt                           <  1:162522/136‑1 (MQ=255)                                                                                                          ccggccATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCtt                           <  8:39432/136‑1 (MQ=255)                                                                                                          ccggccATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCtt                           <  1:372587/136‑1 (MQ=255)                                                                                                            ggccATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTgg                         <  7:335653/136‑1 (MQ=255)                                                                                                                aTCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGc                     <  8:364466/136‑1 (MQ=255)                                                                                                                        cgcgCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCgg              <  5:47141/135‑1 (MQ=255)                                                                                                                        cgcgCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCgg              <  5:170625/135‑1 (MQ=255)                                                                                                                        cgcgCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCgg              <  6:182685/135‑1 (MQ=255)                                                                                                                        cgcgCACATGGACGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCgg              <  4:65729/135‑1 (MQ=255)                                                                                                                        cgcgCACATGGACGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCgg              <  4:65735/135‑1 (MQ=255)                                                                                                                           gcACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTc          >  8:19366/1‑136 (MQ=255)                                                                                                                                 ggTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTcc      >  4:28721/1‑134 (MQ=255)                                                                                                                                   tCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGt  <  8:234425/136‑1 (MQ=255)                                                                                                                                   |                                                                                                                                        GGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGGTACTGGCGCGCCTTGAACACCACCTCCTCGTCGAAATCCAGCTGATTGGTGCGCACGCACGGACCCGCCTGCAGCACCGAGGACATGTTGTGGTACCAGCTTGGACGCACATCGGGTTCGTCCCAGT  >  ACB143/79097‑79362

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