Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,390,331	C→T	Q10Q (CAG→CAA)	psd ←	phosphatidylserine decarboxylase proenzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,390,331	0	C	T	100.0%	218.9 / NA	71	Q10Q (CAG→CAA)	psd	phosphatidylserine decarboxylase proenzyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/31);  total (40/31)

TAACGCCATCGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTCCTCAACCGCTTCCC  >  NC_000913/4390084‑4390572                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   tAACGCCCTCGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGt                                                                                                                                                                                                                                                 >  1:182061/1‑251 (MQ=255)     gCCATCGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCg                                                                                                                                                                                                                                             >  1:219162/1‑251 (MQ=255)           ggcaggcaTGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGATCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCCGAATGTATTGTAGCGAAAGtt                                                                                                                                                                                                                                       >  2:214249/1‑251 (MQ=255)           ggcaggcaTGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGtt                                                                                                                                                                                                                                       >  2:619838/1‑251 (MQ=255)            gcaggcaTGACCAGTACATGCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGGGTCCGGCTTTTGCGCCTCTGTCATGTCGACCTTGTAGTATTTAACGAACAGAG‑CGATAACCAGGTTTGTCAGCCAGCCTGCCCGCTTGCTTGCGCCCCAACCCGCGAGGCGAGTAAGCCGGAGTTTCGGCAGAATGTATTGTAGCGAAAGttt                                                                                                                                                                                                                                      <  2:614425/251‑1 (MQ=255)            gcaggcaTGACCAGAACATTCGGATCGGTATCGAGTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTAGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGTCCTTGTGGGATTTAACGAACAGAT‑CGGTAACCAGTGTTGTCAGCCATCCTGGCCGCTTGCTTGCGCCCCAACCCGCCAGACGAGTAAGCCATAGTTTCGGCAGGATGTATTGTAGCGAAAGttt                                                                                                                                                                                                                                      <  2:218261/251‑1 (MQ=255)             caggcaTGACCAGTGCATGTGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAATGTGCGGTATCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGCCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTa                                                                                                                                                                                                                                     <  1:605918/251‑1 (MQ=255)               ggcaTGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTACCGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTaa                                                                                                                                                                                                                                    >  2:414350/1‑250 (MQ=255)               ggcaTGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTaa                                                                                                                                                                                                                                    >  1:365311/1‑250 (MQ=255)                        agaaCATTCGGATCGGTATCGCTTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTa                                                                                                                                                                                                                           >  2:126550/4‑250 (MQ=255)                           aCATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACa                                                                                                                                                                                                                        <  1:206771/250‑1 (MQ=255)                           aCATTCGGAGCAGTATCGAATGGGCGTACTTCGTCCCGCCACGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACCGAT‑CGATAACCATTTTTGTCAGCCATCCTGCCCGCTTGCTTGAGCCCCAACCCGCCAGGCGAGTAAGCCATAGCTTCGGCGGAATGGATTGGAGCGAAAGTTTAAATGAATTTAACa                                                                                                                                                                                                                        >  2:462197/1‑250 (MQ=255)                                   aTCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGTTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGCCCTTGTAGTATTTAACGAACAGCT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAAACCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGGAAGTTTAAATGAATTTAACAAGGTAGCCt                                                                                                                                                                                                               >  2:446932/1‑251 (MQ=255)                                            gATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTCGCCTCCAGGCCAt                                                                                                                                                                                                      >  2:165561/1‑251 (MQ=255)                                            gATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTGTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCAt                                                                                                                                                                                                      <  2:383057/251‑1 (MQ=255)                                               tGGGCTTACTTCGTCACGCAGCGGTCGGACAACGAATTCGTTAAAGGTGCGGTAGCTGGCGGTTGCCTGCTGTTGCGCCTCTTTCATGTCGTCCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCGTGCGGCCCAACCCGCCAGGCGAGTAAGCCAGAGTTGCGGCAGAATGTATTGTAGCGAAAGTTGAAATTAATTTAACAATGTAGCCTCCAGGCCATTg                                                                                                                                                                                                    <  2:584581/250‑1 (MQ=255)                                               tGGGCGTACTTCGTCACGCAGCGGACCGACAAAGAATTCGTTAAAGGTGCGGGAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTGTTTCAGCCATCCTTCCCGCTTGCGTTCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAAGGTATTGTAGCGAAAGTTTAAATGAATGTAACAAGGTAGCTTCCAGGCCAttt                                                                                                                                                                                                    <  2:49948/250‑2 (MQ=255)                                                           gTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGCGTAAGCCATAGTTTCGGCAGACTGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCt                                                                                                                                                                                       >  1:402641/1‑251 (MQ=255)                                                              aCGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGAt                                                                                                                                                                                    <  1:215584/251‑1 (MQ=255)                                                                              aaGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACAAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCgg                                                                                                                                                                    >  1:304763/1‑251 (MQ=255)                                                                                 aaTTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCGTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGACCAGAG‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTTCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGAt                                                                                                                                                                  <  2:191494/250‑1 (MQ=255)                                                                                           aGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGaaaaaa                                                                                                                                                        >  1:568456/1‑250 (MQ=255)                                                                                                 ggTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGGATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGAc                                                                                                                                                >  2:150137/1‑251 (MQ=255)                                                                                                      cTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGCTTc                                                                                                                                           >  2:549594/1‑251 (MQ=255)                                                                                                        ggCGGTGTCCGGCTTTTGCGCTTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAG‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTTCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTTATTCGGCCTACATGCCGGTTCCTGAAAAAAA‑GGTGACGATTCt                                                                                                                                          <  2:365176/250‑1 (MQ=255)                                                                                                        ggCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACGTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTa                                                                                                                                         <  1:621525/251‑1 (MQ=255)                                                                                                          cGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGCACAGAT‑CGCTAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGATTTGTCGTTCCTGATCCGGACTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAac                                                                                                                                       >  2:273622/1‑251 (MQ=255)                                                                                                              ggccGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGGCAGCCATCCTGCCCGCTTTCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTGTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGAGCCGGCCTACATGCCGGATCCGGAAAAAAA‑GGGTACGAGTCTAacgac                                                                                                                                    <  2:407929/248‑1 (MQ=255)                                                                                                                 cGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGcc                                                                                                                                                                      <  1:522141/214‑1 (MQ=255)                                                                                                                 cGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGcc                                                                                                                                                                      >  2:522141/1‑214 (MQ=255)                                                                                                                             ctctTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCCTTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTgtc                                                                                                                    >  1:591397/1‑251 (MQ=255)                                                                                                                                      gTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCg                                                                                                           <  2:440968/251‑1 (MQ=255)                                                                                                                                          aCCTTGTAGTATTTAACGAACAGAT‑CGATATCCAGTTTTGTCAGCCATCCTGCCCTCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATc                                                                                                       <  1:150137/251‑1 (MQ=255)                                                                                                                                               gtagtaTTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAACAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGaa                                                                                                   >  1:614008/1‑250 (MQ=255)                                                                                                                                                tagtaTTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCTCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCAGCCGTATCAGaaa                                                                                                  <  2:457112/250‑1 (MQ=255)                                                                                                                                                 agtaTTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGaaaa                                                                                                 <  2:33178/250‑1 (MQ=255)                                                                                                                                                 agtaTTTAACGAACA‑ATGCCAGAACAAGGTTTGGCAGCCATACTGCCTGCGTGCTAGCGCCCCAACCCTCCAGGCGAGTAAGCCTTAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAATTGAGCCTCCAGGCCATTGGTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGAGTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAg                                                                                                <  1:355274/251‑1 (MQ=255)                                                                                                                                                  gtaTTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGt                                                                                               >  1:478170/1‑251 (MQ=255)                                                                                                                                                  gtaTTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTGGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGt                                                                                               <  2:365311/251‑1 (MQ=255)                                                                                                                                                    aTTTAACGAACAGAT‑CGATACCCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAacgacg                                                                                                                                   >  1:434817/1‑213 (MQ=255)                                                                                                                                                    aTTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAacgacg                                                                                                                                   <  2:434817/213‑1 (MQ=255)                                                                                                                                                        aaCGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCtt                                                                                                                           <  2:469016/217‑1 (MQ=255)                                                                                                                                                        aaCGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCtt                                                                                                                           >  1:469016/1‑217 (MQ=255)                                                                                                                                                          cGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTAcgc                                                                                       <  1:585077/251‑1 (MQ=255)                                                                                                                                                                gAT‑CGATAACCAGTTTTGTCATCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGttttt                                                                                 <  1:118185/251‑1 (MQ=255)                                                                                                                                                                   ccGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTAc                                                                               <  2:6419/250‑1 (MQ=255)                                                                                                                                                                     gATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCt                                                                             >  2:220867/1‑251 (MQ=255)                                                                                                                                                                     gATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCt                                                                             >  1:380357/1‑251 (MQ=255)                                                                                                                                                                     gATAACCAGTTTTGTCAGCAATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCATACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTaact                                                                             >  1:536692/1‑248 (MQ=255)                                                                                                                                                                      aTAACCAGTTTTGTCAGCCATCCTGACCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGATTTTACGCGTTTTTACCTg                                                                            >  2:81428/1‑251 (MQ=255)                                                                                                                                                                          ccAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTaaa                                                                                                                                                                                                                                   <  2:440320/97‑1 (MQ=255)                                                                                                                                                                          ccAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTaaa                                                                                                                                                                                                                                   >  1:440320/1‑97 (MQ=255)                                                                                                                                                                          ccAGTATTGTCAGCTATACTGCCCGATTGATTGCGCCACAACCCGCAAGGCGTGTAAGCCCTAGTTTCGGCAGACTGTATTGAAGCGAAAGTTTAAATGAATTTAACAAGATAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGCCGCTTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGCACAGATTTTACGCGTTTTTCCCTgcga                                                                         >  2:289479/1‑249 (MQ=255)                                                                                                                                                                             gTTTTGTCAGCCATCCTGCCCGCTTGCTTTCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAC‑GGGGACGATTCTAACGACAGTTAGCTTAATAGTCAGTCATCCGGATCAGAAAAGTTTTTACGCGTTTTTACCTGAGCCATg                                                                     >  2:482232/1‑251 (MQ=255)                                                                                                                                                                                  gTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTGTTTACCTGCGCCATGCTTTc                                                                <  2:123759/251‑1 (MQ=255)                                                                                                                                                                                   tCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTcc                                                               >  1:92047/1‑251 (MQ=255)                                                                                                                                                                                          tCCTGCGCGCTTGCTTGCGCCCCAACCCCCCAAGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCCAGAGTATAACTAAATTTCACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAACAA‑GGGCACGATTCTAACGACGCCTAGCTTACTTGTCAGTCATCCGTATCAGAAACGTTTTTACGCGTTTTTACCTGCGACATGCTTTCCAGAATAc                                                        >  2:303839/1‑251 (MQ=255)                                                                                                                                                                                           ccTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACg                                                       >  1:455157/1‑251 (MQ=255)                                                                                                                                                                                           ccTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGGTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACg                                                       <  2:511271/251‑1 (MQ=255)                                                                                                                                                                                                ccGCTTGCTGGCGCCCCAACCCGCCAGGCGTGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTGTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGAt                                                  <  2:11345/251‑1 (MQ=255)                                                                                                                                                                                                    ttgcttgcGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGtt                                              >  1:571595/1‑251 (MQ=255)                                                                                                                                                                                                       cttgcGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTc                                           >  1:280150/1‑251 (MQ=255)                                                                                                                                                                                                          gcgcCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGa                                         <  2:414417/250‑1 (MQ=255)                                                                                                                                                                                                                 aaCCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGtt                                 <  1:210878/251‑1 (MQ=255)                                                                                                                                                                                                                     cGCCAGGCGAGTAAGCCATAGTTGAGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGAGCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcg                             <  1:478125/251‑1 (MQ=255)                                                                                                                                                                                                                        cAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcgcg                           >  2:229832/1‑250 (MQ=255)                                                                                                                                                                                                                                 aagccctcgTTTCGGCAGAATGTATTGTAGCGACAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAACAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCtt                 >  2:623836/9‑251 (MQ=255)                                                                                                                                                                                                                                      aTAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTCCTc             >  2:105999/1‑250 (MQ=255)                                                                                                                                                                                                                                                cAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTCCTCAACCGCTTccc  >  2:458293/1‑251 (MQ=255)                                                                                                                                                                                                                                                cAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTCCTCAACCGCTTccc  >  1:209631/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCAAGAATACGGTGATAGTTTTCGAACCGGGTTTCCGCGATTTTCCCTTCCTCAACCGCTTccc  >  1:447496/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   TAACGCCATCGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGAT‑CGATAACCAGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAA‑GGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCTTCCTCAACCGCTTCCC  >  NC_000913/4390084‑4390572

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