Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	4,053,476	A→T	8.7%	Q17L (CAG→CTG)	PP_3573 →	putative Monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	4,053,476	0	A	T	8.7%	183.7 / 6.5	57	Q17L (CAG→CTG)	PP_3573	putative Monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (23/28);  new base T (2/3);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

CCACGCCACCCACCCGGCCAACCTGCCCCACAGGCCAGAAAACTCAAGGCCTTCGCCGTTTTATCGGCAAGTACCAGCTAAATACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGCCGGCAAGGACCAGGTGGCCGACTACTTCGAGCAGTACGTACGCAAGTACAACCTG  >  AE015451/4053246‑4053711                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    ccACGCCACCCACCCGGCCAACCTGCCCCACAGGCCAGAAAACTCAAGGCCTTCGCCGTTTTATCGGCAAGTACCAGCTAAATACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGcc                                                                                                                                                                                                                                                 <  1:101831/235‑1 (MQ=255) ccACGCCACCCACCCGGCCAACCTGCCCCACAGGCCAGAAAACTCAAGGCCTTCGCCGTTTTATCGGCAAGTACCAGCTAAATACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGcc                                                                                                                                                                                                                                                 <  1:749323/235‑1 (MQ=255)              ccGGCCAACCTGCCCCACAGGCCAGAAAACTCAAGGCCTTCGCCGTTTTATCGGCAAGTACCAGCTAAATACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGa                                                                                                                                                                                                                                    <  1:397463/235‑1 (MQ=255)               cGGCCAACCTGCCCCACAGGCCAGAAAACTCAAGGCCTTCGCCGTTTTATCGGCAAGTACCAGCTAAATACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCgcc                                                                                                                                                                                                                                           <  2:772335/227‑3 (MQ=255)                                aGGCCAGAAAACTCAAGGCCTTCGCCGTTTTATCGGCAAGTACCAGCTAAATACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGc                                                                                                                                                                                                                  <  1:915682/235‑1 (MQ=255)                                               aGGCCTTCGCCGTTTTATCGGCAAGTACCAGCTAAATACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCAc                                                                                                                                                                                                    >  1:1030166/1‑234 (MQ=255)                                                         cGTTTTATCGGCAAGTACCAGCTAAATACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGa                                                                                                                                                                                                                                >  1:554702/1‑196 (MQ=255)                                                         cGTTTTATCGGCAAGTACCAGCTAAATACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGa                                                                                                                                                                                                                                <  2:554701/196‑1 (MQ=255)                                                                                aaaTACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAAc                                                                                                                                                                                  >  1:565942/1‑219 (MQ=255)                                                                                    ccTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTgg                                                                                                                                                                                                <  1:696300/200‑1 (MQ=255)                                                                                           tCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGc                                                                                                                                                                                                                                              >  1:959707/1‑148 (MQ=255)                                                                                           tCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGc                                                                                                                                                                                                                                              <  2:959706/148‑1 (MQ=255)                                                                                           tCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAgc                                                                                                                                                                                       >  2:462322/1‑203 (MQ=255)                                                                                           tCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGc                                                                                                                                                                               >  2:399720/1‑211 (MQ=255)                                                                                           tCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGc                                                                                                                                                       >  1:1111356/1‑235 (MQ=255)                                                                                             gCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGcc                                                                                                                                                                                                                                                 <  2:543393/143‑1 (MQ=255)                                                                                                        ggCAACCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAgg                                                                                                                                                                     <  1:3609/208‑1 (MQ=255)                                                                                                          cATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAACCCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGc                                                                                                                                                               <  2:996008/212‑1 (MQ=255)                                                                                                                tAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAg                                                                                                                                                                                                                       >  2:454316/1‑150 (MQ=255)                                                                                                                    cTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCTGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTgg‑g                                                                                                                                               <  2:465396/217‑1 (MQ=255)                                                                                                                     tGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAGTCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGTCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAAcc                                                                                                                                                                                 <  2:608170/183‑1 (MQ=255)                                                                                                                       gCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTg                                                                                                                                                                                                                         >  2:295907/1‑141 (MQ=255)                                                                                                                          ccAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCAt                                                                                                                                                                             >  2:992630/1‑182 (MQ=255)                                                                                                                          ccAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTgg                                                                                                                                                                >  2:652497/1‑195 (MQ=255)                                                                                                                           cAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTgg                                                                                                                                                                                                >  2:913708/1‑162 (MQ=255)                                                                                                                           cAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAAcc                                                                                                                                                                                 <  1:203478/177‑1 (MQ=255)                                                                                                                           cAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGcc                                                                                                                                                      >  2:421009/1‑204 (MQ=255)                                                                                                                                    gcggTCCTCCAAGCCGCCCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAGCCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTg                                                                                                                                                                                                          <  1:252931/143‑1 (MQ=255)                                                                                                                                          ctccAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCTTGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCAt                                                                                                                                                                             >  2:748387/1‑166 (MQ=255)                                                                                                                                                       tGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCTGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGt                                                                                                                                                                                               >  2:132206/1‑135 (MQ=255)                                                                                                                                                       tGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCTGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGt                                                                                                                                                                                               <  1:132206/135‑1 (MQ=255)                                                                                                                                                        gACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTTGt                                                                                                                                     >  2:974119/1‑191 (MQ=255)                                                                                                                                                           aaaaaCCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTgg                                                                                                                                                                                                <  2:276071/130‑1 (MQ=255)                                                                                                                                                                     tGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGTTGGCGACCAGCGAGTCCCAGCGGCCGGTACGCCAGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCaa                                                                              <  1:652498/234‑1 (MQ=255)                                                                                                                                                                                 cGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGcc                                                                                                                                                                                                        <  1:256878/100‑1 (MQ=255)                                                                                                                                                                                    aCATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGCCgg                                                      >  1:789159/1‑235 (MQ=255)                                                                                                                                                                                    aCATGACACTGAACAACCTCGAAATCGACACCCTCGGCGTCGGCGCCGGCCCGGCCGGCGTCGCCATCAGCGACCACCTCCGCAAGc                                                                                                                                                                                                                  >  2:223140/1‑87 (MQ=255)                                                                                                                                                                                         acacTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGCCGGCAAgg                                                 <  1:536713/235‑1 (MQ=255)                                                                                                                                                                                               aacaacCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGAt                                                                 <  1:1108951/213‑1 (MQ=255)                                                                                                                                                                                                  cacCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGACCACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCgg                                                                                                                                                        <  2:167889/130‑1 (MQ=255)                                                                                                                                                                                                  aacCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCTGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCa                                                                                                                                                                       >  2:533781/1‑116 (MQ=255)                                                                                                                                                                                                  aacCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCTGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCa                                                                                                                                                                       <  1:533782/116‑1 (MQ=255)                                                                                                                                                                                                  aacCACGAAATCGCCACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGCCGGCAAGGACCAGGTGGc                                       <  1:125606/236‑1 (MQ=255)                                                                                                                                                                                                          aaTCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGCCGGCAAgg                                                 <  2:673378/218‑1 (MQ=255)                                                                                                                                                                                                          aaTCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGCCGGCAAGGACCa                                             >  1:626567/1‑222 (MQ=255)                                                                                                                                                                                                             cGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGCCGGCAAGGACCAGGTGGCCGACTACTTc                             <  1:295908/235‑1 (MQ=255)                                                                                                                                                                                                              gACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGCCGGCAAGGACCAGGTGGCCGACTACTTCGa                           <  1:552302/236‑1 (MQ=255)                                                                                                                                                                                                                 acCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTCCCGGCCGCTgg                                                                                                                                                 >  2:65186/1‑123 (MQ=255)                                                                                                                                                                                                                   ccTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGCCGGCAAGGACCAGGTGGCCGACTACTTCg                            <  1:58377/230‑1 (MQ=255)                                                                                                                                                                                                                        cagtcgGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGcc                                                        <  1:266895/195‑1 (MQ=255)                                                                                                                                                                                                                           tcgGCGCAGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGCCGGCAAGGACCAGGTGGCCGACTACTTCGAGCAGTACGTACGc              <  1:707274/236‑1 (MQ=255)                                                                                                                                                                                                                           gccacgcCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGCCGGCAAGGACCAGGTGGc                                       <  2:408363/207‑1 (MQ=255)                                                                                                                                                                                                                                tcgtgCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGAGGCATTCGCCGGCAAGGACCCGGGGGCCGACTACTTCGAg                          >  1:229057/5‑219 (MQ=255)                                                                                                                                                                                                                                gcCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTg                                                                                                            >  2:316947/1‑137 (MQ=255)                                                                                                                                                                                                                                   ggCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCgggg                                                                                         >  2:602880/1‑153 (MQ=255)                                                                                                                                                                                                                                   ggCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGCCGGCCAGGACCAGGTGGCCGACTACTTCGAGCAgtacgtac                >  1:91242/1‑226 (MQ=255)                                                                                                                                                                                                                                   ggCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGAGGCATTCGCCGGCAAGGACCAGGGGGCCGACTACTTCGAGCAGTACGGACGCAAGTACaa      >  1:315083/1‑236 (MQ=255)                                                                                                                                                                                                                                       aGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGt                                                                                                                                                                                                           >  2:750561/1‑43 (MQ=255)                                                                                                                                                                                                                                       aGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGt                                                                                                                                                                                                           <  1:750562/43‑1 (MQ=255)                                                                                                                                                                                                                                       aGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGCCGGCAAGGACCAGGTGGCCGACTACTTCGAGCAGTACGTACGCAAGTACAACCTg  >  1:1016961/1‑236 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    CCACGCCACCCACCCGGCCAACCTGCCCCACAGGCCAGAAAACTCAAGGCCTTCGCCGTTTTATCGGCAAGTACCAGCTAAATACTATTTTCGCCATTTCCACGGCATCCCTAGTCTGGCTCCAGGCAACGGCGGTCCTCCAAACCGACCTGACAAAAACCGCCTGAACAAGGGCTCGGACATGACACTGAACAACCTCGAAATCGACACCCTCGTCGTCGGCGCCGGCCAGGCCGGCGTGGCCATGAGCGAACACCTGAGCAAGCTTGGCGTGCCGCACCTGGTACTGGAGCGCAACCGCATCGCGCAGGCCTGGCGTACCGGCCGCTGG‑GACTCGCTGGTC‑‑‑GCCAACG‑‑‑‑GCCCGGCCTGGCACGACCGCTTCCCGGGGCTGGAATTCAACCTCGACGCCGATGCATTCGCCGGCAAGGACCAGGTGGCCGACTACTTCGAGCAGTACGTACGCAAGTACAACCTG  >  AE015451/4053246‑4053711

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