Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,538,110	A→C	19.2%	G219G (GGT→GGG)	nanM ←	N‑acetylneuraminic acid mutarotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,110	0	A	C	19.2%	117.3 / 7.4	49	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (32/6);  new base C (2/8);  total (34/15)
Fisher's exact test for biased strand distribution p-value = 2.68e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATAT‑AAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAGCCATTGAAGATATTCTGGTTAACACCACCAGT  >  NC_000913/4537875‑4538340                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           cATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATAT‑AAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑t                                                                                                                                                                                                                                        <  2:188999/239‑1 (MQ=255) cATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATAT‑AAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑t                                                                                                                                                                                                                                        >  1:389224/1‑239 (MQ=255)       cATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATAT‑AAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑Tc                                                                                                                                                                                                                                 >  2:342377/1‑239 (MQ=255)                      tgtGGTAATTTTCTAGTGAAGCTTTGAGTCCGGCCCCTCCGGCAAATAT‑AAGAGACTCCTTGCTTATCTCCGCAAAACCGCAAGCTACGACATCTGGTGATGAGACGGGAGCAGGCTTATTCCCTTTT‑AAATTATGACGGGTCAAATCAGGTT‑C‑ACATACGGCATCCGT‑TCGCAATCGT‑GGTTTGC‑CTCCGCCATTAATAAGCCAGTCTTTATCCCC‑TTTAT‑Tc‑acccgcgcgccacca                                                                                                                                                                                                                 >  1:1065987/3‑226 (MQ=255)                      tCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATAT‑AAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTCTCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCa                                                                                                                                                                                                                >  1:723732/1‑239 (MQ=255)                      tCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATAT‑AAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑ACCCGCCGCACCa                                                                                                                                                                                                                >  2:364209/1‑239 (MQ=255)                      tCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATAT‑AAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCAcc                                                                                                                                                                                                                 >  1:235388/1‑238 (MQ=255)                      tCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATAT‑AAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCa                                                                                                                                                                                                                >  2:1084625/1‑239 (MQ=255)                      tCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATAT‑AAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCa                                                                                                                                                                                                                >  2:577521/1‑239 (MQ=255)                                                 aacccgCCCCCCCGGAAAAAAAAAAGAAAATCATTCCTTACCCCCAAAAAACCGCCAGCACCGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTA‑AAATTATTACCGGTGAAATAAAGTA‑C‑AAATACGGCATCCGT‑CCGCAATCCT‑GGTTGGG‑CTCCGCAATAAATAAGCCAGTTTTTATCCCC‑TTTAT‑CC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                                                                                                      <  2:662594/234‑1 (MQ=255)                                                                               tCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                                                                                                        >  2:683940/1‑239 (MQ=255)                                                                               tCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                                                                                                        >  1:545881/1‑239 (MQ=255)                                                                                      ttATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGAGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                                                                                                                 >  2:846804/1‑239 (MQ=255)                                                                                        aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                                                                                                                >  2:1058485/1‑238 (MQ=255)                                                                                        aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                                                                                                                >  2:64736/1‑238 (MQ=255)                                                                                             cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGTAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCCCC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAACAGCCAACAgca                                                                                                                                          >  1:102988/1‑237 (MQ=255)                                                                                               cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                                                                                                          >  2:429189/1‑237 (MQ=255)                                                                                                     cgccacatgcgtctgcTGGGGAGGATACGGGACAAAGTTTCCTCAATTTA‑AAATTCTTCACGGTAAAAACAAGTA‑A‑AAACCCGGC‑TCCGTCACGCAACCGT‑GTTTTGG‑CTTCGACATTAA‑AACCCCCGTTTTAACCCC‑CTTATCAC‑AC‑ATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                                                                                                  <  1:1058485/224‑1 (MQ=255)                                                                                                                             agaCGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                                                                                                          >  2:775058/1‑239 (MQ=255)                                                                                                                             agaCGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                                                                                                          >  1:269993/1‑239 (MQ=255)                                                                                                                                  ggAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCA‑TAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaca                                                                                                    >  2:900218/1‑237 (MQ=255)                                                                                                                                    gagaaagtttttCCATTTTAAAATTTTTCCGGGAGAAAT‑AAGTTAC‑AAATCGGGCA‑CCGT‑TCCGAAACCC‑TGTTTGG‑GTTCCCCAATAAAAAAGCCGGTTTTTTCCCCTTTTAT‑CC‑AA‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                                                                                                     <  2:289468/228‑1 (MQ=255)                                                                                                                                      cAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAg                                                                                                 >  2:400274/1‑239 (MQ=255)                                                                                                                                              ttCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAtt                                                                                         >  2:107088/1‑239 (MQ=255)                                                                                                                                                              ttaCCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCCTTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCATAATAGTGAGCATTGATTTTATCTATAGc                                                                          >  2:806865/1‑239 (MQ=255)                                                                                                                                                                     tGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGaa                                                                   >  1:1183734/1‑239 (MQ=255)                                                                                                                                                                      taaaaCAAGTT‑AAAAATACGCCACCCGT‑TCGCAATCGT‑GGTTGGG‑CTCCGCCATTAATACCCCAGTTTTTACCACC‑TTTAT‑CC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGa                                                                    <  2:770567/237‑1 (MQ=255)                                                                                                                                                                               tt‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCa                                                          >  2:17411/1‑238 (MQ=255)                                                                                                                                                                               tt‑C‑AAATACGGCATCAGT‑TCGAAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAACTAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCa                                                          >  2:704844/1‑238 (MQ=255)                                                                                                                                                                                  taaatgcggtcATCCGT‑TCGC‑ATCCTGGGTTTGG‑TTGGACCATTAATAGGCCGGTTTTTATCA‑C‑GTTAT‑TCAAC‑ACCCGCCCCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGc                                                        <  1:577521/229‑1 (MQ=255)                                                                                                                                                                                  gtaaatccgccttcttt‑tCTCAATCCT‑TTTTTTG‑TTCCGCCATTAATAGCCCAGTTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGc                                                        <  2:298367/223‑1 (MQ=255)                                                                                                                                                                                   gaaaTCCGGCATCCTT‑TCGCAATCCT‑GGTTGGG‑CTCCGCCATTAATAAGCCAGGTTTAATCCCC‑TTTAT‑CC‑AC‑ACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGcc                                                       <  2:102988/238‑1 (MQ=255)                                                                                                                                                                                           gCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTgag                                               >  1:349399/1‑239 (MQ=255)                                                                                                                                                                                            cacccGT‑CCGCAATCCT‑GGTTTGG‑CTTCGCCATTAAAAACCCAGGTTTTATCCCC‑TTTAT‑CC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTgag                                               <  2:870011/235‑1 (MQ=255)                                                                                                                                                                                                  gtatgCACTCCT‑TGTTTGGCCTTCCCCATTATAAGGCCTGTTTGTCTCCCC‑TTAAT‑TC‑CA‑ATCCGCCGCCCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCtt                                          <  1:848262/235‑1 (MQ=255)                                                                                                                                                                                                     cGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAACTAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCaaa                                      >  1:932218/1‑239 (MQ=255)                                                                                                                                                                                                     cGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCaaa                                      >  2:840486/1‑239 (MQ=255)                                                                                                                                                                                                          tCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAg                                  >  2:209735/1‑238 (MQ=255)                                                                                                                                                                                                                 tttGG‑TTTCGCTATAAATATGCCGGGTTTTATCACC‑TTTAT‑TC‑CA‑ACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAGCCATTGa                           <  2:723732/239‑1 (MQ=255)                                                                                                                                                                                                                 tttGG‑TTCCGCCATAAAGAAGCAGTTTTTTACCCCC‑TTTCT‑CC‑AC‑ACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAGCCATTGa                           <  2:804144/239‑1 (MQ=255)                                                                                                                                                                                                                   tGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAGCCATTGAAg                         >  2:436950/1‑239 (MQ=255)                                                                                                                                                                                                                   tGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAGCCATTGAAg                         >  1:211103/1‑239 (MQ=255)                                                                                                                                                                                                                         tccataaaaatacgcgtgtttgtATCCCT‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAGCCATTGAAGATAtt                    <  1:107088/217‑1 (MQ=255)                                                                                                                                                                                                                         tCGCAATAAATACGCCAGTTTTACCCACC‑TTTAT‑CC‑CC‑AGCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAGCCATTGAAGATAtt                    <  2:236115/239‑1 (MQ=255)                                                                                                                                                                                                                                          ggTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCATAATAGCCATTGAAGATATTCTGGTTAACACCACCAg   >  1:777296/1‑239 (MQ=255)                                                                                                                                                                                                                                           gTTTTTTCTCC‑TTTTA‑TC‑ACAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGTGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAGCCATTGAAGATATTCTGGTTAACACCACCAg   <  2:754179/239‑1 (MQ=255)                                                                                                                                                                                                                                           gTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTgaga                                              >  2:1045051/1‑195 (MQ=255)                                                                                                                                                                                                                                           gTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAGCCATTGAAGATATTCTGGTTAACACCACCAGt  >  2:107993/1‑239 (MQ=255)                                                                                                                                                                                                                                           gTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAGCCATTGAAGATATTCTGGTTAACACCACCAGt  >  1:971343/1‑239 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           CATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATAT‑AAGAGAATCATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACCGGTGAAATCAAGTT‑C‑AAATACGGCATCCGT‑TCGCAATCCT‑GGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAT‑TC‑AC‑AACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTTATCTATAGCGGTTGAATCTTTTCCAGCCTCGTTGAGATCTTCAAAATAGCCATTGAAGATATTCTGGTTAACACCACCAGT  >  NC_000913/4537875‑4538340

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