Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,537,715	A→C	16.0%	V351G (GTG→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,537,715	0	A	C	16.0%	180.5 / 2.8	67	V351G (GTG→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (45/9);  new base C (2/11);  total (47/20)
Fisher's exact test for biased strand distribution p-value = 7.39e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAG  >  NC_000913/4537487‑4537947                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           gcgcATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTg‑c                                                                                                                                                                                                                                  <  2:1066249/238‑1 (MQ=255)    ccttttCCCTCGCCGGTAAGGGCAGTACTTCCCCAAACCCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAAAGATGGGGCGA‑TATTAGTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCA‑CTGAAGCGGTAG‑TTTTGT‑ACTGTGACTTTATTATCCTGCACAGTGATCAAAACTGAATCCGTCACCGCTTTg‑ccggc                                                                                                                                                                                                                              <  2:688215/235‑3 (MQ=255)                         gCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑ATTTGT‑ACTGTGACTTTATTATCCTTAACAGTGATCAAACCTGAATCCGTCACCGCTTTG‑CCGCCTTCAGTTTCACCGCCaataat                                                                                                                                                                                                         >  1:540755/1‑239 (MQ=255)                              ccTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCt                                                                                                                                                                                                                                       <  2:125984/205‑1 (MQ=255)                              ccTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCt                                                                                                                                                                                                                                       >  1:125984/1‑205 (MQ=255)                              ccTTCGCCAGAGGCCATGGCATTAGACTGACCAGCAACGGTAAAAACATAGTAATAATCGGGCAA‑TATTATTGG‑TTTCATATAAATTCCGCGTAATGGAAATGACAAA‑TTG‑ATAGAAGGGCCCCTTTACCGGG‑GCC‑TTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTAAACCCGTCCCCGTTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATa                                                                                                                                                                                                    <  1:1099312/239‑1 (MQ=255)                               cTTCGCAAAACGCCATGGCATTAGACTGACCAGAAACGGTAAAAACATAGTAATAATCGGGCGA‑TTTTATTGC‑GTTCATATAAATTTCGGGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTCCGGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTAATCAAAACTGAATCCGTCCCCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAg                                                                                                                                                                                                   <  1:352974/239‑1 (MQ=255)                               cTTCCCCATACCCCATGGCATTAGACTGACCAGCAACGGTAAGAACATGGTAA‑AAACGGGGGATTATTATTGG‑GTTCATAAAAATTCCGCGTAATTGAAATGAAAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCTTTAG‑TTTTGT‑ACTGTGATTTTATTATCCTTCACAGTGATAAAAACTGAATCCGTCCCCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAg                                                                                                                                                                                                   <  1:93527/239‑1 (MQ=255)                                   gCCATACGCCATGGCTTTAGACTGACCAGCACCGGTAAGAACATAGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTATCAGTGATCAAAACTGAATCCGTCACGGCTTTC‑ACGCCGGCAGTTTCACCCCAAAAAATCAAAAGACta                                                                                                                                                                                               >  2:1066687/1‑239 (MQ=255)                                               ggCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGc                                                                                                                                                                                   >  1:543746/1‑239 (MQ=255)                                                gCATGAGACTGACAAGCAACGGAAAAAAAATAGTAAAAATCGGGCAA‑TATT‑TTGCGGTTCATAAAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑CCC‑CTGAAGCGTTAG‑TTTTGT‑CCTGTGACTTTATTATCCTTCACAGTGATCAAAACTAAATCCGTCCCCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAggg                                                                                                                                                                                    <  1:709018/237‑1 (MQ=255)                                                           aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCCCCGCTTTG‑CCCCCTGCAGTTTCACCCCCAATAATCAATAGACTATTATTCCCGGGCAATGTTaac                                                                                                                                                                          >  2:709018/1‑234 (MQ=255)                                                           aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTg‑ccgcc                                                                                                                                                                                                                              <  2:768535/184‑1 (MQ=255)                                                           aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTg‑ccgcc                                                                                                                                                                                                                              >  1:768535/1‑184 (MQ=255)                                                           aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGCCTATTATTCCAGGGCACTGATACTCCg                                                                                                                                                                       >  1:204705/1‑239 (MQ=255)                                                           aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGCCTATTATTCCAGGGCAATGATACCCCg                                                                                                                                                                       >  1:706139/1‑239 (MQ=255)                                                           aCCAGCAACGGTAAGAACATAGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGGGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAACCAATAAACAATTATTCCAGGGCAATGAAACCCCg                                                                                                                                                                       >  2:13898/1‑239 (MQ=255)                                                                           aCATAATAAAAAGCGGGCAA‑TTTTATTGC‑TTTCATATAAATTTCGCGTAATTAAAATGACAAA‑TTG‑AAAGCAGGCCCCTTTTCCCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑CCTGTGACTTTATTATCCTTCACAGTGATAAAAACTGAATCCGTCACCGTTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCg                                                                                                                                                       <  1:696325/239‑1 (MQ=255)                                                                               aGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATCGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATaa                                                                                                                                                   >  1:467579/1‑239 (MQ=255)                                                                                    taatCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATACTTCACAGTGATCAAAACTGAATCCGCCACCGCTTTG‑CCGCCTGCAGTTTTACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACATTCCGATAAttaa                                                                                                                                               >  1:862287/1‑236 (MQ=255)                                                                                         gggCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TCTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAAATGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGTGATAATTCACCCg                                                                                                                                           >  1:231804/1‑237 (MQ=255)                                                                                           gCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATCCTCCGTAGGCCCGACCTTGCGATAATTAACCCGAttt                                                                                                                                       >  1:685037/1‑239 (MQ=255)                                                                                           gCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGa                                                                                                                                                      >  1:362652/1‑224 (MQ=255)                                                                                                           ttCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACtatta                                                                                                                                                                                            >  1:10409/1‑172 (MQ=255)                                                                                                           ttCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGCTACTCCGTAGGCCCGACCTTCCGATAATTCACCCGATTTATCCCATTTCCCGt                                                                                                                         >  1:25779/1‑239 (MQ=255)                                                                                                           ttCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGCTACCCCCTAGGCCCGACCTTGCGATAATTAACCCGATCTATACCATTTCCCGt                                                                                                                         >  2:79317/1‑239 (MQ=255)                                                                                                           ttCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGCTTTATCCCATTTCCCGt                                                                                                                         >  1:313084/1‑239 (MQ=255)                                                                                                           ttCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGt                                                                                                                         >  1:325168/1‑239 (MQ=255)                                                                                                           ttCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCCCAGGCCCAACCTTGCGATAATTCACCCGATTTATCCCATTTCCGGt                                                                                                                         >  2:646809/1‑239 (MQ=255)                                                                                                           ttCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATAAATAGACTATTATTCCAGGGCAATGATACCCCTTAGCCCGGACCTTGCGTAAATTACCCCGCTTTATCCCATTTccctt                                                                                                                         >  1:148773/1‑237 (MQ=255)                                                                                                                ataAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACATTTCGGAAATTCACCCGCTTTATACCATTTCCCGTTAtac                                                                                                                    >  2:627456/1‑237 (MQ=255)                                                                                                                    attgtgcgTAATTGGAATGTCAAA‑TTGTATACATGTGCCCCTTTA‑CGGG‑GCC‑CTGGA‑CGTTTGATTTTTTGTCTGGGCTTTTATTACCTTACTGATTGAAAAAACCTTAACCCCACCCCCTTTTC‑CCGCTTGACGTTTCCCCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCaa                                                                                                                 <  2:25779/234‑1 (MQ=255)                                                                                                                        cgcgTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATa                                                                                                                                                                            >  1:148626/1‑175 (MQ=255)                                                                                                                        cgcgTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTAATCCCATTTCCCGTTATCCCAAAAAta                                                                                                            >  2:293286/1‑238 (MQ=255)                                                                                                                        cgcgTAATCGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCACTGATCCTCCGTAGCCCCGACCTTGCGTTAATTACCCCGGTTTATCCCATTTCCCTTTAGTCCaaaccta                                                                                                            >  2:681226/1‑235 (MQ=255)                                                                                                                        cgcgTAA‑TGGAATGAAAAATTTG‑ATAGCAGGCCCCCTTTCCCGGG‑CCC‑CGGAAGCGTTAG‑TTTTGT‑ACGGTGATTTTATTACCCTTCACAGTGATCAAAACTGAATCCGTCCCCGCTTGG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATg                                                                                                            <  1:611345/239‑1 (MQ=255)                                                                                                                            aaaTTGAAATGACAAA‑GTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGGTTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTTATAAAAACTGAATCCGTCACCGTTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAata                                                                                                        <  1:293286/238‑1 (MQ=255)                                                                                                                              attttaatttaaaatttG‑ATAGGGGGGCCTCTTT‑CCGGG‑GCCTCGGAAGTTTTAT‑TTTTGT‑ACTGTGATTTTATTATCCTTCACAGTAAACAAAACTAAATCCATCCCTGCTTGC‑CCCACCGAAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAatat                                                                                                       <  2:325168/224‑1 (MQ=255)                                                                                                                                tttaaTGACTAA‑TTT‑ATAGGAGGGCCCCTTTCCCGGGCGCA‑CTAAAGGTTTTT‑TTTTTC‑TCTGTCTCTTTATTTCCCTTCACAGTAACAAAAATTAACTCCG‑CCCCGATGTGCCCGCCTGCATTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCa                                                                                                     <  2:313084/236‑1 (MQ=255)                                                                                                                                    aTGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTACTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTg                                                                                                  >  1:219086/1‑237 (MQ=255)                                                                                                                                       aCAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTACCCCGATTTATCCCATTTCCCGTTATGCCACAGATGAATATCAGTC‑Ttata                                                                                            >  1:603214/1‑239 (MQ=255)                                                                                                                                        cAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCg                                                                                                                                                                       >  1:624125/1‑164 (MQ=255)                                                                                                                                                     cAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGa                                                                                                                                                                                                  >  1:267826/1‑126 (MQ=255)                                                                                                                                                        ggCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCATTAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCTCCCGATTTATCCCATTTCCCGTTATGCCACAGCTTAATATCAGCG‑GTAAAAGGTTTTTACAGGcc                                                                             >  2:723067/1‑239 (MQ=255)                                                                                                                                                        ggCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAgg                                                                               >  1:311387/1‑237 (MQ=255)                                                                                                                                                        ggCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAAAATCAGTG‑CTAAATGATTTTTTAAgg                                                                               >  2:408768/1‑237 (MQ=255)                                                                                                                                                               ttACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTT‑CCGCCTGCCGTTTCCCCGCCAATAATCAATAGACTCTTATTCCCGGGCAATGATACTCCGTAGGCCCGACCTTGCGTTCATTCACCCGGTTTATCCCATTTCCCCTTTTGCCAAAAGTGAAATTCAATGCCTAT‑TTATTTTTTCCCGCCTTCAtcc                                                                      >  2:257013/1‑237 (MQ=255)                                                                                                                                                                  cgggg‑gCC‑TGGAAGCGTTAG‑TTTGTC‑ACTGTGTCTTAATTACCCTCCCCGGTAACCAAACCTAAATCCGGCCCCGTTTTG‑CCGCCGGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATGCGCa                                                                   <  1:646809/237‑1 (MQ=255)                                                                                                                                                                  ccGGG‑GCC‑CTGAAGCGTTTG‑TTTTGT‑ACTTTGATTTTATTATCCTTCACAGTGATAAAAACTGAATCCGTCACGGCTTTG‑CCGCCGGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATGCGCa                                                                   <  1:257013/239‑1 (MQ=255)                                                                                                                                                                  ccGGG‑CCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGAATTTATTATCCTCCACATTGATAAAAACTGAATCCGTCCCCGCTTTG‑CCGCCTGCAGTTCCCCCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATGCGCa                                                                   <  1:408768/239‑1 (MQ=255)                                                                                                                                                                      g‑gCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATGCGCATAGt                                                               >  1:344568/1‑239 (MQ=255)                                                                                                                                                                        gCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAatat                                                                                                       >  1:783583/1‑199 (MQ=255)                                                                                                                                                                         cc‑cTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATACTTCAAAGTGACAAAAACTGAATCAATACCCGATTAG‑CCGACTGCACATACACCGACACAACTCCACAGACAAATATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTCTCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATGCGCATAGt                                                               >  1:594693/1‑237 (MQ=255)                                                                                                                                                                                    ttAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGCTAATTCACCCGATTTATCCCATTTCCCGTTAGGCCAAAGATGAATATCAGTC‑CTAATTGATTTTTTCAGTCCTTCATGCGCCTATTTCTTACCGt                                                      >  2:1060333/1‑236 (MQ=255)                                                                                                                                                                                              t‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATGCGCATAGtt                                                              >  1:1085912/1‑219 (MQ=255)                                                                                                                                                                                                    tgACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTccc                                                                                                                           >  1:120077/1‑154 (MQ=255)                                                                                                                                                                                                    tgACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGa                                                                                         >  1:1039625/1‑187 (MQ=255)                                                                                                                                                                                                          tattatCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCt                               >  1:556505/1‑239 (MQ=255)                                                                                                                                                                                                              atCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGa                           >  1:650548/1‑239 (MQ=255)                                                                                                                                                                                                                  ttCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGCTTTATCCCATTTCCCGTTATGCCAAAGACGAATATCAGTG‑CTATATGATTTTTTCCGGCCTTCATACGCATAGTTCTTACCGTTTTGGAAATTTTTTCGTGAACCTTGGAAAcc                       >  2:541591/1‑239 (MQ=255)                                                                                                                                                                                                                    cacaGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATCCGCa                                                                   >  2:907484/1‑193 (MQ=255)                                                                                                                                                                                                                    cacaGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTAATCCCATTTCCCGTTATGCCAAAGATGAATATCAGCG‑CTATATGATTTTTTCAGGCCTTCATGCGCATATTTCTTACCGTTCTGGTAATTTTCTCTGGAAACtttaacacga                      >  2:607220/1‑231 (MQ=255)                                                                                                                                                                                                                    cacaGTGATCAAAACTAAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATGCGCa                                                                   <  1:907484/193‑1 (MQ=255)                                                                                                                                                                                                                         ttacaaaaaCTGAATCCGTCACCGTTTTC‑GCCTCTGCATTTTCCCCGCCAATAATCAATAGACTATTATACCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCt                <  1:175038/234‑1 (MQ=255)                                                                                                                                                                                                                          aaacAAAACGGAATCCGGCCCCGTTTTC‑CCCCCTGGATTTCAACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCt                <  1:1060333/235‑1 (MQ=255)                                                                                                                                                                                                                                     aaTCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAAtata    >  1:139830/1‑239 (MQ=255)                                                                                                                                                                                                                                       tCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTTCGATACTTCAACCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CCATATGATTTTTTCAGGCCTTCATGCCAATATTCATTGCCGTTCTGGAAATTTTCCTGTTGACCTTTGATTCCGGCCCCTCCCCCaaaacaaaa  >  2:1036328/1‑233 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           GCGCATCTTCCCTGTCCGGTAATGGCAGTCCTTCGCCATACGCCATGGCATTAGACTGACCAGCAACGGTAAGAACATAGTAATAATCGGGCGA‑TATTATTGC‑GTTCATATAAATTTCGCGTAATTGAAATGACAAA‑TTG‑ATAGCAGGGCCCCTTTACCGGG‑GCC‑CTGAAGCGTTAG‑TTTTGT‑ACTGTGACTTTATTATCCTTCACAGTGATCAAAACTGAATCCGTCACCGCTTTG‑CCGCCTGCAGTTTCACCGCCAATAATCAATAGACTATTATTCCAGGGCAATGATACTCCGTAGGCCCGACCTTGCGATAATTCACCCGATTTATCCCATTTCCCGTTATGCCAAAGATGAATATCAGTG‑CTATATGATTTTTTCAGGCCTTCATGCGCATAGTTCTTACCGTTCTGGTAATTTTCTCGTGAACCTTTGAATCCGGCCCCTCCGGCAAATATAAG  >  NC_000913/4537487‑4537947

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