Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	8,963	G→T	L242L (CTG→CTT)	talB →	transaldolase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	8,963	0	G	T	95.1%	178.4 / ‑0.4	61	L242L (CTG→CTT)	talB	transaldolase B
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (22/36);  total (23/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

TCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTG‑ACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACC  >  NC_000913/8719‑9208                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       tCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCAt                                                                                                                                                                                                                                                  <  1:112467/250‑1 (MQ=255)       atgttggTGCTTGGTCGGTAGCGGGCGTGTTCCTGATCCCGCCGTTTGTTGGCCGTCTTCTTGACTGGTACAAAGCGTATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGTCGTGGTTTCTGTATCTGAAATCTACTAGTACGACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGTCCGTCTT‑TCCATCGCAcct                                                                                                                                                                                                                                           <  2:352235/245‑2 (MQ=255)               gCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACActgct                                                                                                                                                                                                                                   <  1:157778/251‑1 (MQ=255)                cTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTAGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACActgctg                                                                                                                                                                                                                                  >  1:422864/1‑251 (MQ=255)                 ttGTGCGGAAGCGGGCGTTTTCCTGATCTCGCCGTTTGTTTGCCGTATGCTTGACTGGTACAAAGCGAATTCCGATAAGAAAGAGGACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACGGTCAGGCTGCGACCTTCTT‑ACCATCGCTCCTACActgctg                                                                                                                                                                                                                                  <  2:345158/250‑1 (MQ=255)                       ggAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTGGGCCGTAATCTTGACTGGTCCAAAGCGAAGACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCGACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGAGCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGc                                                                                                                                                                                                                           <  1:219642/251‑1 (MQ=255)                          aGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACCCTGCTGAAAGAGCTgg                                                                                                                                                                                                                        >  2:314753/1‑251 (MQ=255)                                              cGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAAGTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGc                                                                                                                                                                                                        <  1:441323/247‑1 (MQ=255)                                                     tgttgGCCGTATTCTTGGCTAGTCCAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGAAGAAGATCCAGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACACAGAGCACGGTTATGAACCCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAATTTCTGGAACTGACAGGCGGAGGCCGTCAT‑ACCATCGCACCTTAACTTCTGAAAGAGCTGGATGAGGACGCAGGGGCTATctcacgt                                                                                                                                                                                              >  2:69604/1‑244 (MQ=255)                                                       ttgGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTaa                                                                                                                                                                                            >  1:334059/1‑250 (MQ=255)                                                           ccGTATTCTTGACTGGTACAAAGCGAAGACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCCCCTACTCTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGt                                                                                                                                                                                       <  1:491214/251‑1 (MQ=255)                                                            cGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTc                                                                                                                                                                                      >  2:245930/1‑251 (MQ=255)                                                             gttttCTTGACTGCTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATGCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGTTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCt                                                                                                                                                                                     <  2:640292/248‑1 (MQ=255)                                                               aTTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCtt                                                                                                                                                                                    <  1:129773/250‑1 (MQ=255)                                                               aTTCGTGACTGGTACAAAGCTAATACCGATACGAAAGAGTACGATACTGCAGAAGATCCGGGCGTGTTTTCTGTATCTGAAATCTACCAGTACTAGACAGAGCACGGTTATGAAACCGTGGTTAGGGGCGCAAGCTTCCGTTACATCGGCGTATTTCTTGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACATACACTGCTGAAATAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCtt                                                                                                                                                                                    <  2:409709/250‑1 (MQ=255)                                                                  cTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAAAGCGAAGGGGCTCTCGAACGTAAACTGACTTacac                                                                                                                                                                                >  2:78522/1‑251 (MQ=255)                                                                    tGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTaca                                                                                                                                                                                 <  1:617074/248‑1 (MQ=255)                                                                    tGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACAcc                                                                                                                                                                               <  1:619031/250‑1 (MQ=255)                                                                       cTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCg                                                                                                                                                                           >  1:513840/1‑251 (MQ=255)                                                                       cTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTACGGCAGAAGGTCCGGGCGTGGTTTCTGTATCTAAAATCTGCCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTAAGGGCGCAAGCTTCAGTAACATCGGCGAAATTCTGGGACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCGACACTGCTGAAAGAGCTGGCTTAGATCGAAGGGGCTATCGAACGTAAACTTCCTTACACCGGCg                                                                                                                                                                           >  1:302319/1‑251 (MQ=255)                                                                              aaaGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGaa                                                                                                                                                                    <  1:552061/251‑1 (MQ=255)                                                                                 gCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAACCTGTCTTACACCGGCGAAGTGAAAgc                                                                                                                                                                 >  1:449963/1‑251 (MQ=255)                                                                                    aaTACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTAGCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGcgcg                                                                                                                                                              <  2:669653/251‑1 (MQ=255)                                                                                            tAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGgcgc                                                                                                                                                      >  2:592902/1‑251 (MQ=255)                                                                                                     gTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTg                                                                                                                                             >  1:194028/1‑251 (MQ=255)                                                                                                       ccGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACCAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAg                                                                                                                                           <  2:194028/250‑1 (MQ=255)                                                                                                             cGGCCGAAGAGCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTGATTGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAACCTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAg                                                                                                                                     <  2:264019/251‑1 (MQ=255)                                                                                                              ggCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGCGAGCGCAGGGGCTATCGAACGAAAACTGTCTTACACCGGCGAAGAGAAAGCGCGTCCGGCGCGTCTCACTGAGTCCGAGt                                                                                                                                    >  2:175862/1‑251 (MQ=255)                                                                                                              ggCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGCGTCCGAGt                                                                                                                                    >  1:5365/1‑251 (MQ=255)                                                                                                                cAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTACACTGTCTTACACCGGCGAAGGGAAAGCGCGTCCGGCGCGTAACACTGAGTCCGGGTTc                                                                                                                                  >  2:674092/1‑251 (MQ=255)                                                                                                                              gTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGcaca                                                                                                                     >  1:625837/1‑250 (MQ=255)                                                                                                                                 gTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCa                                                                                                                 <  1:492517/251‑1 (MQ=255)                                                                                                                                    tCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTG‑ACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGcgcg                                                                                                                                                     <  1:251238/212‑1 (MQ=255)                                                                                                                                    tCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTG‑ACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGcgcg                                                                                                                                                     >  2:251238/1‑212 (MQ=255)                                                                                                                                       gTATCTGAAATCTACCAGTACTACAAAGAGCATGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGAt                                                                                                             <  1:515846/249‑1 (MQ=255)                                                                                                                                       gTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGa                                                                                                              <  1:626036/248‑1 (MQ=255)                                                                                                                                        tATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGAGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATcc                                                                                                           >  1:554604/1‑250 (MQ=255)                                                                                                                                           ttGAAATCTCCCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATGCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAAt                                                                                                        <  2:486648/249‑1 (MQ=255)                                                                                                                                              aaaTCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGGCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCa                                                                                                    <  1:640019/251‑1 (MQ=255)                                                                                                                                                             acAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCg                                                                                     <  1:314753/251‑1 (MQ=255)                                                                                                                                                                aaGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATAGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGAACCTACACTGCTGAAAGAGCTGACTGAGAGCGACGGGGCTAACGAACGTAAACAGTCTTAAACCGGCGAAGTGAACGCGCGTCCGGCGCGTATCACTGCGTCCGAATTCCTGTGGCAGCCCAACCAGGATCAAATGGCCGTAGGTAAACTGGCGGa                                                                                   >  2:346714/1‑250 (MQ=255)                                                                                                                                                                     cACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTAt                                                                             <  1:78522/251‑1 (MQ=255)                                                                                                                                                                          ttATGAATCCGTGGTTATGGGCTCACGTGTCCGGCACATCGTCGAACTTCTGGAATTGTCAGTCTGTGACCGT‑TGTACCATCGCATCTACACTGCTGAACTGGCTGGCTGAGAGCAAAGGGGCTATCGCACTTAAGCTGTCTTACATCGGCGAAGTGAAAGCGCGTCCGGCGCTTACCTCTGAGGCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATccta                                                                         <  2:467216/250‑3 (MQ=255)                                                                                                                                                                            aTGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTaa                                                                       <  1:516890/250‑1 (MQ=255)                                                                                                                                                                                         tattggCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGGTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGAc                                                          <  2:109445/246‑1 (MQ=255)                                                                                                                                                                                            gggCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAgg                                                      <  1:592902/251‑1 (MQ=255)                                                                                                                                                                                              gcgcAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGa                                                     >  1:314384/1‑250 (MQ=255)                                                                                                                                                                                                gcAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGaaaa                                                  <  1:26594/251‑1 (MQ=255)                                                                                                                                                                                                gcAAGCTTCCGTAACACCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGaaaa                                                  <  2:113579/251‑1 (MQ=255)                                                                                                                                                                                                            aaCATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGaaaaaa                                       <  1:166882/250‑1 (MQ=255)                                                                                                                                                                                                            aaCATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAAt                                      <  1:346714/251‑1 (MQ=255)                                                                                                                                                                                                             aCATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGTTAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATg                                     <  2:625837/251‑1 (MQ=255)                                                                                                                                                                                                              cATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGa                                    <  1:368330/251‑1 (MQ=255)                                                                                                                                                                                                                      aaaTTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGAt                            >  1:477385/1‑251 (MQ=255)                                                                                                                                                                                                                        aTTCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAATAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATct                          <  2:291601/251‑1 (MQ=255)                                                                                                                                                                                                                         ttCTGGAACTGGCAGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACAGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGCTAAACTGGCGGAAGGTATCCGTAAGGTTGCTATTGACCAGGAAAAACTGGGAAAAATGATCGGCGATctg                         >  1:27872/1‑251 (MQ=255)                                                                                                                                                                                                                                      aGGCTGCGACCGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGcaca                                                                                                                     <  1:447661/146‑1 (MQ=255)                                                                                                                                                                                                                                      aGGCTGCGACCGTAAT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACAGAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGcaca                                                                                                                     >  2:447661/1‑146 (MQ=255)                                                                                                                                                                                                                                                cGTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGAcc  <  1:245930/251‑1 (MQ=255)                                                                                                                                                                                                                                                 gTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCCATGGCAGTAGATaa                                                                                            >  2:608804/1‑160 (MQ=255)                                                                                                                                                                                                                                                 gTCTT‑ACCATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATaa                                                                                            <  1:608804/160‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       TCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTG‑ACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACC  >  NC_000913/8719‑9208

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