Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,111,259	Δ1 bp	coding (718/900 nt)	rfbD ←	dTDP‑4‑dehydrorhamnose reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,111,259	0	G	.	98.5%	287.9 / ‑4.3	65	coding (718/900 nt)	rfbD	dTDP‑4‑dehydrorhamnose reductase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (40/24);  total (40/25)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TCATTTCATCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTG‑AGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTATTTC  >  NC_000913/2111016‑2111501                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   tcatttcatCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTG‑AGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACCAGCGCAACGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTAATGCGTTGAGATTGTTGTGTGCCAGGG‑AATGCCTg                                                                                                                                                                                                                                            >  2:253066/1‑251 (MQ=255)     ttcatCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCACAAACTATTAAATTTCTGTAGTCGTAAATAATTCATTG‑AGCCTTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGGGGACGCCGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGCGCTTGTTGAGTGCAAGGG‑ACTGCCTGCttt                                                                                                                                                                                                                                        >  2:324018/1‑251 (MQ=255)                   tAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTG‑AGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCa                                                                                                                                                                                                                           >  1:689580/1‑250 (MQ=255)                   tAATTCATCTTGATCCACCCTCACGAACAAGATGCAAAACCTATTCACTTGCTGTAGTCGTAAATAATTCATTG‑AGCATTCGTTTAACGCCAACCTGCCAGTCAGGCAAGACACGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTTCTGCGTTGCGCTTGTTTAGTGCAAGGG‑CATGCCTGCTTTGTGCGCCTTTTCa                                                                                                                                                                                                                           >  2:750850/1‑250 (MQ=255)                    aaTTCATCTTGCTCCACCAGCACGAACAAGATGCAAAAACTATTAAATTGCTGGTGTCGTAAATAATCCATTG‑AGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTACaaa                                                                                                                                                                                                                         <  2:461694/251‑1 (MQ=255)                           cTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTG‑AGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAgc                                                                                                                                                                                                                  >  1:285862/1‑251 (MQ=255)                             tGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTG‑AGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAg                                                                                                                                                                                                                   <  1:127198/248‑1 (MQ=255)                                      aTCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTG‑AGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATc                                                                                                                                                                                                       >  1:5945/1‑251 (MQ=255)                                               aagctGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTG‑AGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTCTTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTACTTTGCGCGCCTCTTCCAAAACCAGCGCAGCATAATCGAACCAGGt                                                                                                                                                                                              >  2:445680/5‑251 (MQ=255)                                                  aTGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTG‑AGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTtgt                                                                                                                                                                                           >  1:67522/1‑251 (MQ=255)                                                    gCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTG‑AGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTgg                                                                                                                                                                                         >  2:668217/1‑251 (MQ=255)                                                    gCAAAAACTATTAAAGTGCTGGAGGCGTAAATAATTAAATG‑AGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAAGTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTTTTGAGTGCAAGGG‑AATGCCTGCTTGGCGCGCCATTTCAAAAACCAGCGCAGCATAATCGTACCATGTTGTgg                                                                                                                                                                                         <  2:307469/251‑1 (MQ=255)                                                     cAAAAACTATTAAATTGCTGTAGTCGTAACTAATTCATTG‑AGCATTCGTTTCACGACAACCTGCCAGTCAGGCAAGACAAGCGCAACGTTCTGCTGCAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTCGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGGGTGCAAGGG‑AATGCCTTCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACAAGGTTGTggg                                                                                                                                                                                        >  2:473668/1‑250 (MQ=255)                                                        aaactactatactGCTGTAGCCTTAAATAATTC‑TTGTAGTATCCGTTGCTCGCCAACTTGCCAGGCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTTTATTAAGGCGAGAGTTATGTCGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTAAAGCGTTGTGCTTGTTGAGTTCAAGGG‑AATGCCTGCTTTGCTCGCCTCTTCAAAAACCAGCGCAGCATTATCGTCCCAGGTTGTGGTatc                                                                                                                                                                                     <  2:779458/239‑3 (MQ=255)                                                                  aTTGCTGTAGTCGTAAATAATTCATTG‑AGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCCAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTAACACTGGCTAcc                                                                                                                                                                           >  1:558085/1‑251 (MQ=255)                                                                                  ataatTCATTG‑AGCATTCGTTTCACGCCAACCTACCAGTCAGGCAAGACATGCGCACAGTTCTGCTGACATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGCTAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGCGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCTGGTTGTGGTACCACTGCCTACCACATGGAACAAGCCTg                                                                                                                                                           >  2:80791/1‑251 (MQ=255)                                                                                     atTCATTG‑AGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCg                                                                                                                                                         >  1:70701/1‑250 (MQ=255)                                                                                                cATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGttt                                                                                                                                              >  1:464313/1‑251 (MQ=255)                                                                                                 aTTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGGAAAAGCCTGCGACATCCGGtttt                                                                                                                                             >  2:240450/1‑250 (MQ=255)                                                                                                       ttCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAg                                                                                                                                        <  2:5945/250‑1 (MQ=255)                                                                                                        ccACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTCGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGt                                                                                                                                       <  2:715122/249‑1 (MQ=255)                                                                                                         cACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCACAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGc                                                                                                                                     >  1:674565/1‑251 (MQ=255)                                                                                                                 ccTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCCCTGGCTACCAAATGGTACAAGCCTGCGCCATCCGGTTTATTCAGAGCGacacaa                                                                                                                              >  2:63418/1‑248 (MQ=255)                                                                                                                      cAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGACGAGAGTTATGTGGACGAAGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGGGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTAGCTACCACATGGTACAAGACTGCGACATCCGGTTTATTCAGTGCGACACGCATGGc                                                                                                                         >  2:729971/1‑250 (MQ=255)                                                                                                                        gTCAGGCAAGACAAGCGCAAAGCTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGAACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCAtg                                                                                                                      >  2:113635/1‑251 (MQ=255)                                                                                                                         tCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGAGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGCGTGCAAGAG‑ACTGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAAGTTGTGGTACCACTGGCTACCAAATGGTACAAGCATGCGACATCCGGCTTATTCAGTGCGCCACgcatggcatgg                                                                                                                     >  1:263711/1‑250 (MQ=255)                                                                                                                          cAGGCAAGACAAGCGCAAAGTTCTGCTGTAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTgg                                                                                                                                                                                         <  2:723909/182‑1 (MQ=255)                                                                                                                          cAGGCAAGACAAGCGCAAAGTTCTGCTGTAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGAG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTgg                                                                                                                                                                                         >  1:723909/1‑182 (MQ=255)                                                                                                                              cAAGACAAGCGCAAAGTTCTGCTGCAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGt                                                                                                                >  2:5530/1‑251 (MQ=255)                                                                                                                              cAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTTTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGt                                                                                                                <  1:60783/251‑1 (MQ=255)                                                                                                                                   cAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAACTGGTACAAGCCTGCGAAATCCGGTTTATTCAGTGCGACACGAATGGCCTGTGCTGTACAAt                                                                                                           >  1:645454/1‑251 (MQ=255)                                                                                                                                   cAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAAt                                                                                                           >  2:462926/1‑251 (MQ=255)                                                                                                                                   aaaGCTCAAAGTGCTGCTGAAATTTTTCTGTATTAAGGCGAGTGTTATGTGGACTACGAGCTGGGTTAGGATAGGCTGTTGTTGGTACTGCTTTGAGCTTGTTGAGTGCAAGGG‑TATGCCTGCGTTGCGCGCCTCTTCAAAAACCATCGCAGCATAATCGTCCCAGGTTGTGGTTCCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACTCGAATGGCATGTGCTGTACAAt                                                                                                           <  2:558085/250‑1 (MQ=255)                                                                                                                                            aGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATcagccagc                                                                                                   >  1:306363/1‑250 (MQ=255)                                                                                                                                            aGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCcagca                                                                                                  <  1:462926/251‑1 (MQ=255)                                                                                                                                               tCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACTAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTTTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAACCCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTATAATCAGCCAGCAGt                                                                                                <  1:324018/250‑1 (MQ=255)                                                                                                                                                     gAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGc                                                                                           <  1:292197/249‑1 (MQ=255)                                                                                                                                                             tCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAACCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTgg                                                                                  <  2:250110/250‑1 (MQ=255)                                                                                                                                                              cTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGCTAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACa                                                                                                                                                                 >  2:57713/1‑170 (MQ=255)                                                                                                                                                              cTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACa                                                                                                                                                                 <  1:57713/170‑1 (MQ=255)                                                                                                                                                                  aTTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCATGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGGTGTACAATCAGCCAGCAGTTAGGCACATGGTGGCGCAcc                                                                            >  2:359650/1‑251 (MQ=255)                                                                                                                                                                     aaGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCaaa                                                                         <  1:573206/251‑1 (MQ=255)                                                                                                                                                                         cGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGGGGTACCGCTGGCTACCAAATGGTACAAGCCTTCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGa                                                                     <  1:750850/251‑1 (MQ=255)                                                                                                                                                                          gagaGTTATGTGGACGACGAGCTGGGGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGATGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGACCCAAACTGAt                                                                    >  2:116904/1‑251 (MQ=255)                                                                                                                                                                                   gtgGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCGTTTGGCGCCCCAAACTGATCGTtaata                                                            >  2:102087/1‑250 (MQ=255)                                                                                                                                                                                         gacgaGCTGGCGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACAAGAATGGCATGTGCTGTACAATCAACCAGAAGTTCAGCACCTGTTGGCGCACCAAACTGAACGATAATAACCGCt                                                      >  2:493887/1‑250 (MQ=255)                                                                                                                                                                                          acgaGTTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCGGCATAATCGTTCCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTaa                                                    <  2:321313/251‑1 (MQ=255)                                                                                                                                                                                            gaGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAAtt                                                  >  1:514614/1‑251 (MQ=255)                                                                                                                                                                                            gaGCTGGTGTAGGATAGGCTGTTGTTGGTACGGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCGTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAAtt                                                  <  2:70701/251‑1 (MQ=255)                                                                                                                                                                                             aGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGCGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAAttc                                                 >  1:149026/1‑251 (MQ=255)                                                                                                                                                                                                tGGTGCAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCCAAAACCAGCGCAGCCTAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCTACATCCGGATTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAACCAGCAGTTCAGCACCGGTTGGCGCACCAAACTAATAGTTACTAACCGCTAAtttttc                                              >  2:37596/1‑251 (MQ=255)                                                                                                                                                                                                    gtAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGc                                          >  2:170078/1‑251 (MQ=255)                                                                                                                                                                                                    gtAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACGGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGc                                          <  2:67522/251‑1 (MQ=255)                                                                                                                                                                                                      gggATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGctc                                        <  2:259727/250‑1 (MQ=255)                                                                                                                                                                                                             gCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAACACCAGCGCAGCATAAACGTCCCAGGTTGTGGACCCACTGGCTACCAAATGGTACACGCCTGCGACATCCGGTTTATTCAGTGCGACACGACTTGCATGTTCTGTACAATAATACAGCAGTACCACACCTGTTGGCGCCCCAAACTGCGCGTTAATAACCGCTAATTCATCACTCTCttgtacc                                 >  2:479398/1‑246 (MQ=255)                                                                                                                                                                                                               tgttgtTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGcc                                 >  1:201626/1‑249 (MQ=255)                                                                                                                                                                                                               ggtttgtgGTACTGCGTTGAGCTTGTTGAGTG‑AAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAg                               <  2:345703/245‑1 (MQ=255)                                                                                                                                                                                                                     tggtactgcttggagttgggggggTGCAAGGG‑AATGCCTGCTTTGCGTGCCTCTTCAAAATCCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCGGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTTTACAATCAGCTAGCAGGTCAGCACCTGTTGGCGCATCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCGTTTTCCAGACGTa                          <  2:551609/227‑1 (MQ=255)                                                                                                                                                                                                                           tGCGTTGAGCTTGTTGAGTGCAAGGG‑CATGCCTGCTTTGCGCGCCTCTTCACAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGACTGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCTAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGt                   >  1:197812/1‑251 (MQ=255)                                                                                                                                                                                                                             cGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGttt                 >  2:260404/1‑251 (MQ=255)                                                                                                                                                                                                                                 gAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGc             <  1:351160/251‑1 (MQ=255)                                                                                                                                                                                                                                      tgttgAGTGCAAGGG‑AATGCCTGCTGTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAg         <  1:170078/250‑1 (MQ=255)                                                                                                                                                                                                                                       gttgAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGtt       <  1:63418/251‑1 (MQ=255)                                                                                                                                                                                                                                         tgAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCATGAGACATCCGGTTTCTTCAGTGCGACACGAGTGGCAAGTGCTGTACAATCAGCCATCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAAACGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTAt     >  1:772720/1‑251 (MQ=255)                                                                                                                                                                                                                                            gTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTAtttc  >  2:553736/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   TCATTTCATCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTG‑AGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTATTTC  >  NC_000913/2111016‑2111501

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