Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	8,975	2 bp→TA	coding (738‑739/954 nt)	talB →	transaldolase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	8,975	0	G	T	95.0%	184.9 / ‑0.4	60	P246P (CCG→CCT)	talB	transaldolase B
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (21/36);  total (23/37)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01
*	NC_000913	8,976	0	G	A	96.6%	173.1 / ‑2.1	60	A247T (GCA→ACA)	talB	transaldolase B
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (21/36);  total (23/37)
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.94e-01

GCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTG‑ACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTC  >  NC_000913/8733‑9216                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                                 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)                                                                                                                                                                                                                                           cATCGCACCTACACTGCTGAAAGAGCTGGCTGAGAGCGAAGGGGCTATCGAACGTAAACTGTATTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCGGCTGTAATCATTCTTAGCGTGACCGGGAAGTc  <  2:314384/251‑1 (MQ=255)                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                                 GCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTG‑ACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTC  >  NC_000913/8733‑9216

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: ‑