Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,342,825	C→T	N305N (AAC→AAT)	melA →	alpha‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,342,825	0	C	T	97.9%	155.9 / NA	47	N305N (AAC→AAT)	melA	alpha‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (30/16);  total (31/16)
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

CGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAACTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGC  >  NC_000913/4342581‑4343060                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            cGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCAt                                                                                                                                                                                                                                       >  2:144803/1‑251 (MQ=255)           ggcggcTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTg                                                                                                                                                                                                                              <  1:144803/250‑1 (MQ=255)           ggcggcTGATGAAGCAGGGAAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCt                                                                                                                                                                                                                               <  2:96409/249‑1 (MQ=255)              ggcTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGGCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTggagg                                                                                                                                                                                                                          >  1:59494/1‑251 (MQ=255)                   aTGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGtat                                                                                                                                                                                                                     <  1:444390/251‑1 (MQ=255)                    tGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGtat                                                                                                                                                                                                                     >  1:246349/1‑250 (MQ=255)                      tagcagaccaggcctcGAAGCCGAATATTAATGGTAATCCTCACTACCAGAATATTGTGCGCTCCGAAATGTTCAAAAAGCTGGGCTATTACGTCACGGAATCGTAAGAACATTTTGCTGAGTACAAACCATGGTTTATTAAGCCAGGCCGAGAGGATTTGATTGCGCGTTATACAGTACCGCTGGATGCGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAACTGGCATAAAGAGCTGGAGGCGTATaa                                                                                                                                                                                                                   >  2:198792/16‑250 (MQ=255)                                ggcaCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTccc                                                                                                                                                                                                        >  1:84181/1‑251 (MQ=255)                                       aaGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGa                                                                                                                                                                                                 <  1:153840/251‑1 (MQ=255)                                         gCCGACTATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGCTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATGTTGCTGAGTACACACCGTGGTTTATGAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGATCAGTTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCGCCCGGATTGata                                                                                                                                                                                               <  2:177372/251‑1 (MQ=255)                                           cGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGGAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGACTTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATAtt                                                                                                                                                                                             >  1:190872/1‑251 (MQ=255)                                                    aTGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTGCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCggca                                                                                                                                                                                    <  2:59484/251‑4 (MQ=255)                                                      ggCAATACTCGCTGCCAGAATATTGGGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACg                                                                                                                                                                                  >  2:77150/1‑251 (MQ=255)                                                          aTACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGTTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGa                                                                                                                                                                               >  2:166219/1‑250 (MQ=255)                                                          aTACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGaa                                                                                                                                                                              <  1:148598/251‑1 (MQ=255)                                                                      aGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAAAATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGcac                                                                                                                                                                   >  1:310303/1‑250 (MQ=255)                                                                        aaTATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAAt                                                                                                                                                                >  1:346193/1‑251 (MQ=255)                                                                         atatTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACATACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATc                                                                                                                                                               <  1:97758/251‑1 (MQ=255)                                                                                     aCGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGCGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATACTGAACGCTATc                                                                                                                                                   >  2:315524/1‑251 (MQ=255)                                                                                          aTGTACAAAACGCAGGGCTATTTAGTCAAGGAATCGTGAGAACATATTCATGAGTACACACCGTGGTTTATTAAGCCAGGTCGAGAGGATTTGATTGAGCGTTATAAAGCACCGCTGGATGAGCACCCGAAACGCTGCGTCGAGCAGCTAGCAAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCACCCGGAATGATATAAAACCGTCACGGGAATATGCCAGAACAATCATGCACGCAATCTGGAc                                                                                                                                              >  2:458267/1‑251 (MQ=255)                                                                                              tCAAAAAGCTGGGCTCTTTCGTCACGGAATAGTCAGAACATTTTGCTGAGTACACACCGTGTTTTATTAAGCCAGGACGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAGAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACtagc                                                                                                                                          >  2:367637/1‑248 (MQ=255)                                                                                                aaaaaGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATaaaa                                                                                                                                                                                                                 <  2:129144/178‑1 (MQ=255)                                                                                                aaaaaGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATaaaa                                                                                                                                                                                                                 >  1:129144/1‑178 (MQ=255)                                                                                                   aaGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGcc                                                                                                                                     >  1:242979/1‑251 (MQ=255)                                                                                                      cggggCTATTGCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGAGTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAg                                                                                                                                  <  1:381817/249‑1 (MQ=255)                                                                                                           cTATTTCGTCACGGAATCGTCAGAGCATTTTGCTGAGTACACACCTTGGCTTATTAAGAAAGGGCGAGAGGATTTGATTGAGCGATATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGCGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCACGAACTCTCTCTGGACTGGcacgccgaaggtg                                                                                                                              >  2:337056/1‑238 (MQ=255)                                                                                                               ttCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGCATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGGGAttt                                                                                                                          >  1:332646/1‑250 (MQ=255)                                                                                                                  gTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGAACCCGAAACGCTGCGTCGAGCACCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGCCACGGGAATATGCCACCACAATCCTGAACTCTCTCTGGACTGGccacccgcgcccgacttcccc                                                                                                                      >  1:290622/1‑231 (MQ=255)                                                                                                                        gAATCGACAAAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGAGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTCCCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTCTAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGCGTGTGATTTACGTCAACGt                                                                                                                >  2:214231/1‑251 (MQ=255)                                                                                                                         aaTCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTc                                                                                                               >  1:164997/1‑251 (MQ=255)                                                                                                                         aaTCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTc                                                                                                               >  2:105600/1‑251 (MQ=255)                                                                                                                                    cATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATgg                                                                                                    >  1:308974/1‑251 (MQ=255)                                                                                                                                            tGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTttgattg                                                                                            >  1:152077/1‑251 (MQ=255)                                                                                                                                            tGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTttgattg                                                                                            >  1:251459/1‑251 (MQ=255)                                                                                                                                             gAGTACACAACGTGGTTTATTAAGCCTGGTCGTTAGGATTTTATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGTAATATGCCAGCACAATCATGAACGCTATCTGGGCTGGCGAGCCGTGTGTGATTTACGTCAACGTCCGTAACGATGTTttgatt                                                                                             <  2:201329/249‑1 (MQ=255)                                                                                                                                                     acCGTGGTTTATTAAGCCAGGTCGCGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGCGTACCCGAAACGCTGCGTCGAGCAGCTGGCGACTTGGCATACAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGCTATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTCTCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGGAACGATGGTTTTATTGATAACCTGc                                                                                   >  2:174298/1‑251 (MQ=255)                                                                                                                                                       cGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGCGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTttgatt                                                                                             >  1:48258/1‑239 (MQ=255)                                                                                                                                                       cGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTttgatt                                                                                             <  2:48258/239‑1 (MQ=255)                                                                                                                                                       cGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGcc                                                                                  >  1:113658/1‑250 (MQ=255)                                                                                                                                                                   ggccagctgtggAGGATTGGATTGAGCGTTATAAAGTACCTCTGGATGATGACCCGAAACGCTGCGTGGAGCAGCTGGCGAATGGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGAGTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCGGAGTGGGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAggtttttgc                                                                     <  2:40780/240‑8 (MQ=255)                                                                                                                                                                               aGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACCATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGcc                                                         >  1:298879/1‑251 (MQ=255)                                                                                                                                                                                    ttgattgaGCTTTAGAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGTTTTATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTc                                                     <  2:408581/250‑1 (MQ=255)                                                                                                                                                                                     tgattgaGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACATGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTg                                                   >  2:27801/1‑251 (MQ=255)                                                                                                                                                                                       attgaGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCGGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCATGGAAGTAGCCTGTCTGGt                                                 <  2:346193/251‑1 (MQ=255)                                                                                                                                                                                                              ggaggAGTACCAGAAACGCTGCGTCGAGCAGCTGTCGAATTGAAATAAAGAGCTGGCGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGACAGCACAAACATGAACGCTATCTGGACTGGCGAGCCGCGTGTGATTTACGGCCACGTCCGTAACGATGGGTTGATTGATAACCTGCCACAATGATGTTACGTGGAAGTAGCCTGTCTGGTTGTGGCTAATGGCATTCAGCCg                           >  2:402157/5‑250 (MQ=255)                                                                                                                                                                                                                                  gCGTCGAGCAGCTGGCGAATTGGCATAAAGAGCAGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTACCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCt      <  1:77150/251‑1 (MQ=255)                                                                                                                                                                                                                                      cGAGCAGCTGGCGAATTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGACGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGc  >  2:154847/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            CGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAACTGGCATAAAGAGCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGC  >  NC_000913/4342581‑4343060

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