Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,715,661	C→T	D77N (GAT→AAT)	pdxY ←	pyridoxal kinase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,715,661	0	C	T	97.6%	132.9 / ‑3.7	42	D77N (GAT→AAT)	pdxY	pyridoxal kinase 2
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (23/17);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

CACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTG  >  NC_000913/1715427‑1715897                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              cacaGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTaa                                                                                                                                                                                                                                   >  1:375951/1‑247 (MQ=255)      attATTTCCAGCTCAACCAGATTTGCCGCAATTATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACCCCCGGTTCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCGGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAAt                                                                                                                                                                                                                          <  1:363693/249‑1 (MQ=255)            tCCAGCTCAACCAGATTTGCCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGCTCGCAAGAATATTTCGCCTGCGGATTCGCGGCTTTCCCCTGACGGACGATACCGAGGATCTGTTCACC‑‑CTGCTCCACCGATCCCCGATAGCCACTTAATACGGCATTACAGGGGTGTCATTTATCAATgtcggc                                                                                                                                                                                                                    >  2:332285/1‑251 (MQ=255)                       cAGATTTGGCGCAATGATACCGGTTGCAGGCAGACGGTCCCGCACAGGAAACTCGGCTACCCCCGGTGCAACGATACAGCCTTTTCCCGGATGACCCATTCCCGGTTCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCCCCTGACGGACGATACCGAGGAAATGTCCTCC‑‑CTGGTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGca                                                                                                                                                                                                         <  2:554593/251‑1 (MQ=255)                              ggcgctaggATATCCCTGTCAGGCAAACCGTGCCCTACATGAACCCCTGCGACCCCCGGTGCAATGATCCGCCCTTTTTCCGGATGCCCCTTTCCCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTGTCCCCTTACGGACGTTACTGAGG‑T‑TTTTCACCCTCTTCTCCGCCGATCCCAGAGAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTc                                                                                                                                                                                                 <  2:325145/243‑1 (MQ=255)                                    aTGATATCGCTGGCAGGCAAACCGTGCCGCCCATGAAACTCTGCGACCCCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGTATTCGCGGCTTTCCCCTGACGGACGATACCGAGGATATGTTCCCC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCTATGGCGGCAATGCCTTGCACA‑ATTTCGGTTaa                                                                                                                                                                                           <  1:182586/251‑1 (MQ=255)                                         atCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGc                                                                                                                                                                                      >  1:156302/1‑251 (MQ=255)                                                  aggcaAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGggcggc                                                                                                                                                                              >  1:557127/1‑250 (MQ=255)                                                       aaCCGTCCCCCACATGAAACCCTGCAACACCCGGGTCACCAATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCCCCTGACGGACGATACCGAGTATATGTTCCCC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACg                                                                                                                                                                        <  2:557127/251‑1 (MQ=255)                                                          cGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGCTGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGTTACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCCGATAGCCACTTAATACGGCATTACAGGAGTGTAATTTATAAATGGCGGCAATGCATTGCACA‑ATTTCGGTAAAATGGCTAGGCGGCATCACGCa                                                                                                                                                                      >  2:72029/1‑250 (MQ=255)                                                              ccGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCAAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGcc                                                                                                                                                                   >  1:460704/1‑249 (MQ=255)                                                               cGCACATGACACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGCTCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATCCCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTTC‑CACATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGt                                                                                                                                                                >  2:373144/1‑251 (MQ=255)                                                               cGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGt                                                                                                                                                                >  1:325503/1‑251 (MQ=255)                                                                              gCGACACCCGGTGCCACGATACAGCCTTTTTCCGGATGACCCATTACCGGCTCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTAt                                                                                                                                                   <  1:87289/249‑1 (MQ=255)                                                                                  cgcccgggtcAACGATACATCCTTTTTCCGGATGCCCCATTCCCGGATCGCAAAAATATTTCGCCTGCGGCGTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGgtg                                                                                                                                             <  1:455292/242‑1 (MQ=255)                                                                                    cccGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGtgtg                                                                                                                                           >  1:311590/1‑251 (MQ=255)                                                                                       ggcgaaaCGCTACAGACTTTTTCAGGATGACCAATTAACGGATCGCAAAAATATTTCGCCTGCGGATTCGCGACTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CGGCTCCGCCGATACCAGATCGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGAACA‑ATTACGGTTACATGGCTGAGCGGCATCCCGCAGCCAGTCCATTTGCCGGATTGGGTGGGAtt                                                                                                                                        >  2:557041/6‑251 (MQ=255)                                                                                            aaCGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATTCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAgga                                                                                                                                    >  2:102359/1‑248 (MQ=255)                                                                                              cGATACAGCCTTTTCCCGGATGACCCTTTACCGGTTCGCAAAAATATTTCGCCTGCGTATTCGCGTCTTTCCCCTGACGGACGATACCTAGTATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAAt                                                                                                                                 <  1:232455/251‑1 (MQ=255)                                                                                                                      ccATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACg                                                                                                         <  1:373144/251‑1 (MQ=255)                                                                                                                          tacctgaacccAAAAATATTTCGCCGGCTGATTCCCGCCTTTCCCCTGACGGACGATCCCAAGTATTTGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGACTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTc                                                                                                      <  1:338964/240‑1 (MQ=255)                                                                                                                          tACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTc                                                                                                      <  1:460913/250‑1 (MQ=255)                                                                                                                              ggATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCCTCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGGGTTCAGCGGCCAGACGTTCGCGcc                                                                                                 >  2:447107/1‑251 (MQ=255)                                                                                                                                  cGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTCTTCAGCGGCCAGACGTTCGCGCCCAgg                                                                                             >  2:395917/1‑251 (MQ=255)                                                                                                                                      aaaaTATTGCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggcggc                                                                                         <  1:311574/251‑1 (MQ=255)                                                                                                                                        aaTATTTAGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGAATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGAGTTCAGCGGCCAGACGTTCGCGCCCAGGCGgcgc                                                                                       >  1:355518/1‑251 (MQ=255)                                                                                                                                         atatTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGaa                                                                                                                                    >  1:538015/1‑205 (MQ=255)                                                                                                                                            tttCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTAACC‑‑CTGCTCCGCCGATACCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCCTCACGAAGCCAGTCCATTTGCCGTATTGGGGTTGGTTAGAAAATTGAAGGGAGTTCCGGGGCCAGACGGTCGCGCCCAAGCGGCGCATCg                                                                                   >  2:368552/1‑251 (MQ=255)                                                                                                                                                   tGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTACATGGCTGGGCGGCATCACGCAGCCAGTCCCTTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCCGCGGCCAGACGTTCGCGCCCAGGCGGCGCAGCggtacctc                                                                            >  2:181478/1‑245 (MQ=255)                                                                                                                                                    gCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAAct                                                                             <  1:395917/249‑1 (MQ=255)                                                                                                                                                     cGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTGACGCCCAGGCGGGGCATCGGAAActctt                                                                          >  2:390316/1‑250 (MQ=255)                                                                                                                                                            gcgGCTTTCAACTGACGGACGCTACAGAGGCTATGTTCACC‑‑CTGCTCCGCCGCTCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCAGTATTGGGTGTGATTAGAAACTTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCCGGCGCCACATCGGaacatgacccgcact                                                                   >  2:70988/1‑238 (MQ=255)                                                                                                                                                                             gacgATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCAt                                                  >  1:108747/1‑251 (MQ=255)                                                                                                                                                                             gacgATACCGAGGATAGGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCAt                                                  <  1:478324/251‑1 (MQ=255)                                                                                                                                                                                       aGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCCGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCCCTGTTACCCGCATGACCATAAACAACGTg                                        >  1:85106/1‑251 (MQ=255)                                                                                                                                                                                                    cccctctcccgccgATCCCAGATAGCCACTTAATACGGTATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCGGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAg                             <  2:341353/243‑1 (MQ=255)                                                                                                                                                                                                         gCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCACGGCGGCGCATAGGAAACTCTTCCACGCTGTTCCCCGCATGACCATAAACAAAGTGAGACTGGGTAGCg                           >  2:339537/1‑248 (MQ=255)                                                                                                                                                                                                                     ccAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCACTGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCCTCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAACCAACGTgag                                      >  2:263652/1‑225 (MQ=255)                                                                                                                                                                                                                     ccAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTgag                                      <  1:263652/225‑1 (MQ=255)                                                                                                                                                                                                                     ccAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTtcatc            <  1:339537/251‑1 (MQ=255)                                                                                                                                                                                                                      cAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAAAGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCAAGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCACAGGCGGCGCATCGGAAACACTGCCGCACTGTTACCAGCATGACCATACACAACGTGAGACTGGCTAGCGAGAATATTTTtcatc            >  1:389346/1‑250 (MQ=255)                                                                                                                                                                                                                                attAATACGGCATGACAGGTGTGTAAGTTATCAATGGTGGCAAGGCCTTGCACA‑ATTCCGGTTAAATGTCTGGGCGGCAGCACGCCGCCAGTCCATTTGCCGTATTGGGTGTTATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGGTCCCTg  <  2:401103/249‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              CACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACC‑‑CTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACA‑ATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTG  >  NC_000913/1715427‑1715897

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