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%	147.4 / NA	46	Q10Q (CAG→CAA)	psd	phosphatidylserine decarboxylase proenzyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/25);  total (21/25)

AACGCCATCGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCT  >  NC_000913/4390085‑4390556                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                   aaCGCCATCGGCAGGCATGACCAGTACATTCGTTTCGGTATCGATTGGGCGTACTTCGTCACTCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTATCTGGCGGTGTCCGGCTGTTGCGCCTCTTTCATGTCGCCCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCGTGCGCCCCAACCCGTTAGGCGAGTAAGCCATATTTGCGGCAGAATGTATTGTa                                                                                                                                                                                                                                <  2:300281/251‑1 (MQ=255) aaCGCCATCGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGt                                                                                                                                                                                                                                 >  1:230482/1‑250 (MQ=255)      cATCGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGt                                                                                                                                                                                                                                 >  2:571247/1‑245 (MQ=255)      cATCGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGt                                                                                                                                                                                                                                 <  1:571247/245‑1 (MQ=255)      cATCGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGa                                                                                                                                                                                                                            >  1:251893/1‑250 (MQ=255)       atgggcagcgatgcccAGTAAATGCGGATTGGGATCTATTGGGCGTGCCTCGTCACGCAGAGGTCGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCGTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTGAACGAACAGATCGATAAACAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGTGCCCTAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAAGGTATTGTAGCGaaa                                                                                                                                                                                                                          <  2:190016/237‑1 (MQ=255)             aggcaTGACCAGTACATTCGGATCGGTATCGCTTGGGCGAACTTCGTCACACCGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCAGTGTCCGACTTTTGCGCCTCTTTCATGTCGCCCTTGTAGTATTTAACGAACAGCTCGATACCCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATCGTTTCGGCAGACTGTATTGTCGCGATAGTTTa                                                                                                                                                                                                                     >  1:250250/1‑250 (MQ=255)                 aTGACCAGGACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGa                                                                                                                                                                                                                <  2:154204/251‑1 (MQ=255)                   gACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAAt                                                                                                                                                                                                              <  2:251893/251‑1 (MQ=255)                   gACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAAt                                                                                                                                                                                                              <  1:258059/251‑1 (MQ=255)                       aGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTTCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTaa                                                                                                                                                                                                          <  1:123694/251‑1 (MQ=255)                       aGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTa                                                                                                                                                                                                           >  2:341154/1‑250 (MQ=255)                       aGTAAATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTGTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTaa                                                                                                                                                                                                          <  2:304906/251‑1 (MQ=255)                            aTTCGGATCGGGATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCTACCTTGTAGTCTTTAACGAACAGTTCGATAACCCGTTTTGTCAGCCAAC‑CTGCCCGCTGGCTTGCGCCCCAACCCGCCCGGCGAGGAAGCCATAGTGTCGTAGGACTGTTTTGTAGCGCAAGTTTAAATGAATTTAACAAgg                                                                                                                                                                                                     >  2:188963/1‑251 (MQ=255)                                      gTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAg                                                                                                                                                                                           <  1:678925/251‑1 (MQ=255)                                                   gTACTTCGTAACGCATAGGACGGCCACAGACTTCGTTAAAGGTGCGGTAGCTGGCGGAGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATACACAGCTTTGTCCGCCATC‑CTGCCCGCTTGCTTACGCCCCAACCCACCAGGCGAGTAAGCACTAGTTTCGGAAGAATCTATCGTAGAGCAAGTTCACATGAATTTAACAAGGTAGCCTCCCTGCCATTGTTTTGt                                                                                                                                                                              >  2:642963/1‑251 (MQ=255)                                                             aCGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGAt                                                                                                                                                                    <  1:663528/251‑1 (MQ=255)                                                                               gAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGCTAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGCATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGAt                                                                                                                                                  >  1:439318/1‑251 (MQ=255)                                                                                                   aGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGtt                                                                                                                                                                                                                       <  2:640185/162‑1 (MQ=255)                                                                                                   aGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGtt                                                                                                                                                                                                                       >  1:640185/1‑162 (MQ=255)                                                                                                       ggCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGCTCCTGAAAAAAAGGGG‑ACGATTCt                                                                                                                          >  2:253039/1‑250 (MQ=255)                                                                                                          ggTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGCTCGATAACCAGTTTTGTCAGCCAAC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAacg                                                                                                                      >  1:93891/1‑251 (MQ=255)                                                                                                          ggTGTCCGGCGTTTGCGCCTCGTTCATGGCGACCTTGTATTACTTAAAGAACAGATCGATAACCAGTTTTGTCAGCCTTC‑CTGCCCGCTTGCTTGCGCCCCAACCCGTCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAacg                                                                                                                      <  2:530387/251‑1 (MQ=255)                                                                                                                       tGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTa                                                                                                          <  1:665812/250‑1 (MQ=255)                                                                                                                                        gACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTAt                                                                                        <  1:188963/251‑1 (MQ=255)                                                                                                                                        gACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGAAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCTTGATCCGGCCTACATGCCGGCTCCTGAAAAAAAGGGG‑ACGCTTCTAACGACGGGTAGCTTAATTGTCAGTCATACGTAt                                                                                        >  2:64614/1‑251 (MQ=255)                                                                                                                                             tGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCa                                                                                      <  1:61268/249‑1 (MQ=255)                                                                                                                                             tGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGa                                                                                    <  1:453508/250‑1 (MQ=255)                                                                                                                                             tGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTGTCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGa                                                                                    <  1:606197/250‑1 (MQ=255)                                                                                                                                                agtaTTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGaaaa                                                                                 <  1:526110/250‑1 (MQ=255)                                                                                                                                                     gtagctatcAGAGCGAGAACCAGTGTTGTCAGCC‑TCGCTGCCCGCTTGCGTGCGCCCCCACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTATCTAAAGTTTAAATGGATTTAACAAGGTAGCCTCCAGGCCATTGGTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCGGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGtttt                                                                            <  2:648890/242‑1 (MQ=255)                                                                                                                                                         cGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTCTTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGCAAAAAAGGGG‑ACGATTCTAACGCCGGTTAGCTTAATTGTCGGTCATCAGTATCATATAAGTTTTTActc                                                                       >  2:277416/1‑249 (MQ=255)                                                                                                                                                          gAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCAATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGATAAGTTTTTAcgcg                                                                      >  2:382964/1‑251 (MQ=255)                                                                                                                                                           taCAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGTTTACTTGTAAGTCATCCGCATCAGAAAAGATTTTAcgagt                                                                     >  2:526197/2‑248 (MQ=255)                                                                                                                                                                      aaCCAGTTTTGTCAGCCATC‑CTGTCCGCTGGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTAcc                                                              <  2:469915/247‑1 (MQ=255)                                                                                                                                                                      aaCCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTAcc                                                              >  1:469915/1‑247 (MQ=255)                                                                                                                                                                           gTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATg                                                     <  1:19981/251‑1 (MQ=255)                                                                                                                                                                              ttGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGaaaaaaa                                                                                                                                       <  2:537495/167‑1 (MQ=255)                                                                                                                                                                              ttGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGaaaaaaa                                                                                                                                       >  1:537495/1‑167 (MQ=255)                                                                                                                                                                                 tCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTc                                                <  1:408196/250‑1 (MQ=255)                                                                                                                                                                                                        tgcGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTTTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTc                           <  1:636821/249‑1 (MQ=255)                                                                                                                                                                                                                    cGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcg             >  1:600474/1‑251 (MQ=255)                                                                                                                                                                                                                           cGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGAtttt      >  1:353435/1‑251 (MQ=255)                                                                                                                                                                                                                           cGAGTAAGCCATAGTTTCGGCAGAATGTATTGTAGCGAAAGATTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATACGGCCTACATGACGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCCTGCTTTCCAGAAAACGGTGATAGTTTTCGAAACGGGTTTCcgagctttt      >  2:222395/1‑244 (MQ=255)                                                                                                                                                                                                                               tAAGCCATAGTTTCGGCAGAATGTGTTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCt  >  2:201088/1‑251 (MQ=255)                                                                                                                                                                                                                                aatccATAGTTTCGGCAGAATTTATGGTAGCGAAAGTTTAAATAAATTTAACAAGGTAGCCTCCAGGCCAGTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACTGTGAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCt  <  1:78037/247‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                   AACGCCATCGGCAGGCATGACCAGTACATTCGGATCGGTATCGATTGGGCGTACTTCGTCACGCAGCGGACGGACAAAGAATTCGTTAAAGGTGCGGTAGCTGGCGGTGTCCGGCTTTTGCGCCTCTTTCATGTCGACCTTGTAGTATTTAACGAACAGATCGATAACCAGTTTTGTCAGCCATC‑CTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTGAAAAAAAGGGG‑ACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGATTTTCCCT  >  NC_000913/4390085‑4390556

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