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	100.0%	200.1 / NA	58	N305N (AAC→AAT)	melA	alpha‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/33);  total (25/33)

CGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAACTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGC  >  NC_000913/4342581‑4343060                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             cGGAACTGCTGGCGGCTTATGAAGCCGGGCCGGCGCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGATGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGCTTTGATTGAGCGTTATAAAGTACAGCTGGTAGAGTACCCGAAACGCTGCCTCGAGCAGCTGACGAGTTGGCAt                                                                                                                                                                                                                                        >  2:110369/1‑251 (MQ=255)  ggAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATa                                                                                                                                                                                                                                       <  2:419247/251‑1 (MQ=255)        gctgGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑Gc                                                                                                                                                                                                                                <  1:176576/251‑1 (MQ=255)                           gggcaggCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAAAGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAgt                                                                                                                                                                                                             <  2:6401/251‑2 (MQ=255)                                ggcaCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTccc                                                                                                                                                                                                        <  1:490944/251‑1 (MQ=255)                                ggcaCCGAAGCCGAATATTCATGGCAATACACGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCGGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTccc                                                                                                                                                                                                        <  1:115044/251‑1 (MQ=255)                                 gcaCCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTCCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCg                                                                                                                                                                                                       >  1:502200/1‑251 (MQ=255)                                     cGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGAtt                                                                                                                                                                                                   <  1:611948/251‑1 (MQ=255)                                           cGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGACATGTTCAAAAAGCTGGGCTATTTCGTCTCGGAATCGTCAGAACATCTTGCTGAGTACACACCGTGGTTTATTAAGCCCGGTCGTGAGTATTTGATTGAGCGTTATAAAGTCCCGCTGGATGAGTTCCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATAtt                                                                                                                                                                                             <  1:653/251‑1 (MQ=255)                                              atatTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTaaa                                                                                                                                                                                          <  1:643891/251‑1 (MQ=255)                                                atTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAAcc                                                                                                                                                                                        >  2:502709/1‑251 (MQ=255)                                                atTCATGGCAATACTCGCTGACAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGTGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGACAGGTCGTGAAGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATCAAGA‑GCTGGAGGAGTTTACAAAAGACTCCCGGAATGATATTAAAcc                                                                                                                                                                                        >  2:264696/1‑251 (MQ=255)                                                  gcATGGCAATACGCGCTGCCAGAATATTTTGCGCGACTAAATGTCCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGATGATGTGATTGAACGTTATAAAGTACCGCTGTGTGAGTACCCGAAACGCTAAGTCGAGCATCGTGCGGATTGTCATAAAGA‑GCTGGGGGAGTATAAAAAAGCCTTCCTGATTGAGATTATACCGt                                                                                                                                                                                      <  2:696134/250‑1 (MQ=255)                                                   cATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTc                                                                                                                                                                                     >  2:273582/1‑251 (MQ=255)                                                    agtgCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGGCAGAACATTTTGCTGAGTGCACACCGTGGTTTATTAAGCCAGGTCGGGAGCATTTGATTGAGCGATATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTAGAGCAGCTGGCGAATGGGCATAAAGA‑GCGGGAGGAGTATaaaaaa                                                                                                                                                                                                               <  2:333938/221‑1 (MQ=255)                                                    aTGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGACTCGTCCGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGCTTTGATTGAGCGTTATACAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGCGTATaaaaaa                                                                                                                                                                                                               >  1:333938/1‑224 (MQ=255)                                                                 cTGCCAGAATATTGTGCGCTACGACATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGCACATTTTGCTGCGTACACCCCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAACAAAGCCTCCCGGATTGATATTAAACCGTCCCGGGAATATGcc                                                                                                                                                                        >  2:357809/1‑250 (MQ=255)                                                                           atTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCAt                                                                                                                                                             <  1:62588/251‑1 (MQ=255)                                                                                      cggaaTGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCt                                                                                                                                                  <  1:576616/248‑1 (MQ=255)                                                                                      cGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCt                                                                                                                                                  <  2:380312/251‑1 (MQ=255)                                                                                               cAAAAAGCTTGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGGGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGc                                                                                                                                          <  1:254633/250‑1 (MQ=255)                                                                                                aaaaaGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTgccg                                                                                                                                         >  1:18931/1‑247 (MQ=255)                                                                                                aaaaaGCGGGGCTATTGCGTCACGGAATCGGCAGAACATTTTGCTGAGTACACACCGTGGTTTAGTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGa                                                                                                                                        <  2:1283/251‑1 (MQ=255)                                                                                                   aaGCTGGGCTATTTCGACACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGAATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGcc                                                                                                                                     >  2:727907/1‑251 (MQ=255)                                                                                                        gggCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACGCACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGGGCCGAgtg                                                                                                                                <  2:146907/251‑1 (MQ=255)                                                                                                                tCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTAc                                                                                                                        >  2:341584/1‑251 (MQ=255)                                                                                                                  gTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACgg                                                                                                                      <  1:658806/251‑1 (MQ=255)                                                                                                                    cACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGCGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGCGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGCGTGTGATTTACGGcc                                                                                                                    >  2:497694/1‑250 (MQ=255)                                                                                                                        gAATCGTCAGACCATTTTGCTGAGTACAAACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGATGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAAACGTCACGGGAATATGCCAGCACAATAATGAACGATATCTGGACTGGCGAGCCGAGAGTGATTTACGGAAACGt                                                                                                                >  2:364109/1‑251 (MQ=255)                                                                                                                         aaTCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTc                                                                                                               >  2:232583/1‑251 (MQ=255)                                                                                                                            cGTCAGAACCCTTTGCTGAGTACAAACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAGAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGCGCAGCTGGCGACTTGGCATAAAGA‑GCTGGAGGAGTCTACCAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGATTTCTGGACTGGCGAGCCGAGTGTGATTTACGGCACCGTCCgg                                                                                                            >  2:515394/1‑250 (MQ=255)                                                                                                                            cGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCATGTCGTGAGGATTTGATTGAGCGTGATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATAGTAAACCGGCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGATTGGCGAGCCGAGTGAGATTTACGGCAACGTCCGt                                                                                                            <  2:460260/251‑1 (MQ=255)                                                                                                                            cGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACAGCAACGTCCGt                                                                                                            >  1:505101/1‑251 (MQ=255)                                                                                                                             gTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACTGCAACGTCCGTa                                                                                                           <  2:18931/251‑1 (MQ=255)                                                                                                                                    cATTTTGCTGAGTTCATACCGTGGTTTATTAAGCCAGGTCGTGATGATTTGATTGAGCGTTATAAAGTACCGCTGTATGAGTACCCGAAACGCTGCGTCGAGCAGCTTGTGAATTGGCATAAAGA‑GCGGGAGGAGGATAAAAAGGCCTCCCGGATGGATATGAATCCGTCACGGGAATATGCCAGCACAATCATGAACGCGATCAGGACTGGCGAGCCGAGTGTGATTTACGGCAACTTCCGTAACGATgg                                                                                                    <  1:477288/251‑1 (MQ=255)                                                                                                                                       tttGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGtt                                                                                                  <  1:597125/250‑1 (MQ=255)                                                                                                                                                tACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATCGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATa                                                                                         >  1:163818/1‑250 (MQ=255)                                                                                                                                                    cacCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGAGTTCTAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAACAAGCCTCACGGATTGATATTAAACCGTCACGGGAATATGCAAGAACAATCATGCAAGCTATCTGGACTGGCGGGCCGAGTGTGCTTTACGGCAACGTCCGTAACGCTGGTTTGATTGATAACCTg                                                                                    >  2:673220/1‑251 (MQ=255)                                                                                                                                                    aacCGTAGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGGGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTg                                                                                    >  2:78163/2‑251 (MQ=255)                                                                                                                                                                  aaGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCGGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTg                                                                      <  1:53813/251‑1 (MQ=255)                                                                                                                                                                    gCCAGGTCGTGAGGATTTAATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCg                                                                    <  1:247725/251‑1 (MQ=255)                                                                                                                                                                      cAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTg                                                                  >  2:298637/1‑251 (MQ=255)                                                                                                                                                                                   tttGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTc                                                     >  2:338247/1‑251 (MQ=255)                                                                                                                                                                                       attgaGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGt                                                 >  1:629959/1‑251 (MQ=255)                                                                                                                                                                                       attgaGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACAGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGt                                                 <  2:505101/251‑1 (MQ=255)                                                                                                                                                                                          gaGCGTTATAAAGTACCGTTGGATGAGTACCCGAAACGCTGCGTCGAGAAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCAGGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGttt                                               <  2:408679/250‑2 (MQ=255)                                                                                                                                                                                                     aGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCAGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGGGGAAGTAGCCTGTCTGGTTGATGCTAATGGc                                    >  2:3373/1‑250 (MQ=255)                                                                                                                                                                                                       tACCGCTGGGTGATTACCCGAAACGCTGCGTCTAGCATCTGGCGAATTGGCATAAAGA‑GCTGAAGTAGTATAAAAAAGCCTCCCGCATGTATATTAAACCGTCACGGGAATATGCCAGCACATTCATGAACGCTATCTGGACTTGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATTGTTTGATTGATAACCTGCCTCAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCag                                  <  1:364109/250‑2 (MQ=255)                                                                                                                                                                                                        aCCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTc                                >  1:437170/1‑251 (MQ=255)                                                                                                                                                                                                         ccGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATGGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTAGCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCa                               <  1:357809/251‑1 (MQ=255)                                                                                                                                                                                                              ggATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGcc                            <  1:502709/249‑1 (MQ=255)                                                                                                                                                                                                                       cccGAAACGCTGCGTCGAGCAGCTGGCGAATCGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCg                 <  2:163818/251‑1 (MQ=255)                                                                                                                                                                                                                              cgcggCGTCGATCAGCTTGCGAATTGGCATCAAGA‑GCGGGAGGTGTATAAAAACGCCTCCCGGATTGATATTAAACCGGCCCGGTAATAGGCCAGCACAATCATGAACGTTCTCTGGGCTGGCGAGCCGAGGGTGATTTACGGCAACGTCCGGAACGAGGGTTTGATTGATATCCTGCCTCAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGACGCTAATGGCATTCAGCCGACCAAAGTCGTTACGCt          <  2:462335/247‑1 (MQ=255)                                                                                                                                                                                                                                  gcgacGAGCAGCTGGCGAATTTGCATAAAGA‑GCTGGAGGATTATAAAAAAGCCTCCCGGATTGATAGTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTAcc       <  2:585575/246‑1 (MQ=255)                                                                                                                                                                                                                                  gCGTCGAGCTGCTGGCGAATTGGCCTAA‑GACGCTGGCGGACTATCAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACACTCATGAACGCTATTTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGGTGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCt      >  2:363846/1‑251 (MQ=255)                                                                                                                                                                                                                                   cGTCGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCtt     <  1:526682/251‑1 (MQ=255)                                                                                                                                                                                                                                     gcGAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCg   <  1:232583/250‑1 (MQ=255)                                                                                                                                                                                                                                       gAGCAGCTGGCGAATTGGCATAAAGA‑GCTGGAGGAGTATAACAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGACAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGc  >  2:594486/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             CGGAACTGCTGGCGGCTTATGAAGCAGGGCAGGCACCGAAGCCGAATATTCATGGCAATACTCGCTGCCAGAATATTGTGCGCTACGAAATGTTCAAAAAGCTGGGCTATTTCGTCACGGAATCGTCAGAACATTTTGCTGAGTACACACCGTGGTTTATTAAGCCAGGTCGTGAGGATTTGATTGAGCGTTATAAAGTACCGCTGGATGAGTACCCGAAACGCTGCGTCGAGCAGCTGGCGAACTGGCATAAAGA‑GCTGGAGGAGTATAAAAAAGCCTCCCGGATTGATATTAAACCGTCACGGGAATATGCCAGCACAATCATGAACGCTATCTGGACTGGCGAGCCGAGTGTGATTTACGGCAACGTCCGTAACGATGGTTTGATTGATAACCTGCCACAAGGATGTTGCGTGGAAGTAGCCTGTCTGGTTGATGCTAATGGCATTCAGCCGACCAAAGTCGGTACGCTACCTTCGC  >  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: ‑