Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,067	A→G	Q155Q (CAA→CAG)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,067	0	A	G	98.4%	205.5 / ‑4.4	64	Q155Q (CAA→CAG)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (37/25);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-01

GGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCG  >  NC_000913/4619825‑4620313                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        ggcaTGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATATCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCCTTCTGGATCAGCTGGGCgg                                                                                                                                                                                                                                                 >  1:265870/4‑250 (MQ=255)   aaTGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAAGGTCATCCGGTAAAGATACCCCGTCTGTTCACTGGGAAATCGCCGGTGTCCCTGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAag                                                                                                                                                                                                                                               <  1:450264/251‑1 (MQ=255)       gACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCAATGGGAAATCGCCGGTGTCCCGGTTCTGTTTGGGTGGGGGTATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGcc                                                                                                                                                                                                                                           >  1:292921/1‑250 (MQ=255)          ggCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACAt                                                                                                                                                                                                                                        <  1:213379/251‑1 (MQ=255)                gCTGAAGTTATCGGCGCGTACGCATGGGCGCACGGAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTACTCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGCAGAGCACATGaaacc                                                                                                                                                                                                                                  >  2:130413/1‑249 (MQ=255)                         aTCGGCGCGTACGCATGGGCGCCCGAACTGTCCTCCGGTAAAGATACCCCGGCTGGGAACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGACGAACGCGCTAATATGCCGGGGTACATCGGTAAC‑TGCCACTCTTCCGGGACGGTCATTCTGGATCAGATGGGCGCAGTGCACATGAAAACCgggcagcc                                                                                                                                                                                                                         >  2:301163/1‑245 (MQ=255)                              ggcgTACGCATGGGCGCCCGAAATGTCATCCGGTAAAGATCCCCCGTCGGGTCTCTGGAAAATCGCCGGTGTCCGGTTTCGGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGCGTTGCTGGATAACCTGGTCGAACGCGCTAGTCTGCCGGGTTACACCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGAttt                                                                                                                                                                                                                    <  2:662486/250‑1 (MQ=255)                                 gTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCt                                                                                                                                                                                                                 >  1:345062/1‑251 (MQ=255)                                   aCGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGTAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGTTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCta                                                                                                                                                                                                                <  1:486503/250‑1 (MQ=255)                                    cGCATGGGCGCACGAAAGGCCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTTTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCtata                                                                                                                                                                                                              <  1:351673/251‑1 (MQ=255)                                                   aaTGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCg                                                                                                                                                                                               >  1:350011/1‑251 (MQ=255)                                                     tGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCgtg                                                                                                                                                                                             >  1:508139/1‑251 (MQ=255)                                                       tCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGtt                                                                                                                                                                                           >  1:648744/1‑251 (MQ=255)                                                        cATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGtt                                                                                                                                                                                           <  1:36325/250‑1 (MQ=255)                                                        cATCCGGTAAAGATACCCCGTCTGGTCACTGG‑AAATCGCCGGTGTCCCGGTTCTGTTTGCGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGATAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGTCGAAGAGCACATGAAAACCGGCAAGCCGATTTTATATACCTCCGCTGACTCCGTGTTcc                                                                                                                                                                                         >  2:25639/1‑251 (MQ=255)                                                               taaagctgccccgacaggTCACTGGCAAATCGCCGCTGGCCCGGTTCTGTTTGAGAGGGGCCATTTCTCCGATCACGAAAACAGCCTCCCGCAAGAGCTCCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCAATCTGGATCAGCTGGGCGAAGCGCACATAAAACCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCAGTGTTCCAGAtt                                                                                                                                                                                    >  2:277375/17‑250 (MQ=255)                                                                 aaGATACCCCGTGTGGTCACTGGGAAATCGCCGGTGTCCCGGTTCTGGTTGAGTGGGGATATTTCCCCGATCACGAAAACAGCTTCCCGCAAGCGCTGCTGGATAAACTGGTCGAACGCGTTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATtgcc                                                                                                                                                                                 <  2:153044/251‑1 (MQ=255)                                                                       ccccGTCTGGTCACTGGGAAATCGCAGGTGTCCCGGTTCTGTTTGAGTGGGGATCTTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAA‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAACTTGCCTGCCac                                                                                                                                                                           >  2:501122/1‑250 (MQ=255)                                                                               ggccACTGGGAAATCGCTGGTGTCCCGTTTTTGTTGGAGTGTGGATATTTCGCCGATCACTAAAACGGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCGCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACgaagg                                                                                                                                                                      <  2:236002/245‑2 (MQ=255)                                                                                                         ggTTCTGTTTGAGTGGGTATAGTTCTCCGATCACGAAAACAGCTTCCCGCAAGTGCTGCTGGATAAACTGTTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTCCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAAc                                                                                                                                         <  1:61272/251‑1 (MQ=255)                                                                                                                     ggggggATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCGGCCGGGTTACCTGGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCg                                                                                                                             <  1:606614/249‑1 (MQ=255)                                                                                                                       ggggATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACgaag                                                                                                                                                                       >  1:394962/1‑207 (MQ=255)                                                                                                                          gATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGCTCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTAGATAAACTCTACGAACTGTGAGAAATCGCCCgg                                                                                                                        >  2:169618/1‑250 (MQ=255)                                                                                                                          gATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGt                                                                                                                        >  1:363562/1‑251 (MQ=255)                                                                                                                              tttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGCTTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGa                                                                                                                      >  2:319516/1‑249 (MQ=255)                                                                                                                                        cACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCa                                                                                                           >  1:172842/1‑250 (MQ=255)                                                                                                                                         aCGAAACCAGTTTCCCCCAAGAGCTGCGGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAc                                                                                                         <  1:25639/251‑1 (MQ=255)                                                                                                                                         aCGAAAACAGCTTCCCGCAAGGGCTGGTGGATAAACTGTTCGATCGCGCTAATCTGCCGGTTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGTTCATTCTGGATCATCTGGGCGAAGATCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTTGATAAACTCTACGATCTGTGCGAAATCGCCCGTGAAGAGCTGACCAAc                                                                                                         <  1:102605/251‑1 (MQ=255)                                                                                                                                           gAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcgg                                                                                                       >  1:224817/1‑251 (MQ=255)                                                                                                                                           gAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcgg                                                                                                       >  1:21994/1‑251 (MQ=255)                                                                                                                                            aaaaCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAACCCGGCAAGCCGATTTTCTCTACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGCGCTGACCAAaggc                                                                                                      >  1:50065/1‑251 (MQ=255)                                                                                                                                                 aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCt                                                                                                  >  1:345111/1‑250 (MQ=255)                                                                                                                                                 aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCt                                                                                                  >  1:388344/1‑250 (MQ=255)                                                                                                                                                 aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCt                                                                                                  >  1:46657/1‑250 (MQ=255)                                                                                                                                                  gCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGCAGGAACCTTTGGCCTGGATAAAATCTACGAACTGTGCGAAATCGCCCGTGCAGAGCTGACCAACGGCGGATAc                                                                                                >  2:48938/1‑251 (MQ=255)                                                                                                                                                   cTTCCCGCAAGAGCTGCTGGCTAAACTGGCCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGACACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTACGAAATCGCCCGTGAAGCGCTGACCAACGGCGGCtcc                                                                                                >  2:651832/1‑248 (MQ=255)                                                                                                                                                   cTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTAc                                                                                                >  1:106666/1‑250 (MQ=255)                                                                                                                                                      cccGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACaacc                                                                                            >  1:172989/1‑249 (MQ=255)                                                                                                                                                           aaggtttgctgGATAACCGGGCCGACCGCGCTAATGTGCCGG‑TGTCCCCCGTAACCTGCCACTCTTCCGTTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACGCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGt                                                                                       <  2:21994/245‑1 (MQ=255)                                                                                                                                                                ctaatagATAAACTGGTCGAACGAGCTAATCTGCCGGGTTACCTCGGTAAA‑TGCCACTCTTCCGGTACGGGCATTCTGGATCAGCTGGGCGCAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGGTCAAGATTGCCTGCCACGAAGGAACCTTGGGCCTGAATAAACTCTACGAACTGTGCGGAGTCGCCCGTGAAGAGCGGACCAACTGCTGCTaacccatcggtagggt                                                                                  >  2:160310/7‑236 (MQ=255)                                                                                                                                                                   gtgGATAAACGGGCCGAACGCGTCAATCTGCCGGTTTGCCTCGGAAAC‑TGACACTCTTCCGGTACGGTCATTCTGGATCAGCGGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TAt                                                                              <  2:508139/250‑1 (MQ=255)                                                                                                                                                                   ctgGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcggc                                                                                                      >  1:146440/1‑228 (MQ=255)                                                                                                                                                                      gATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCCAAAACATCGGTCGGGT‑TATCGc                                                                           >  2:294656/1‑251 (MQ=255)                                                                                                                                                                            cgggTCGAACGCGCTAATCGGCCGGGTTACCTCGGTAAC‑TGCAACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTcc                                                                     <  1:508257/249‑1 (MQ=255)                                                                                                                                                                            cTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCa                                                                                                                                                                                        >  1:267370/1‑137 (MQ=255)                                                                                                                                                                             tGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCg                                                                    >  1:564585/1‑251 (MQ=255)                                                                                                                                                                                 caaccgcgcagttcttccGGGTTCCCCCTGTAAC‑GGCCCCTCTTCCGGGACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCGATACCTCCGCTGCCCCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTGTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCTTCCGTTTa                                                                <  2:350011/235‑1 (MQ=255)                                                                                                                                                                                          tAATCTACCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGCGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAACCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCCACACAATCGGTCGGGTGTA‑CGCCCATCCGTTTATCGGTGAc                                                        >  2:587606/1‑250 (MQ=255)                                                                                                                                                                                           aaTCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGT‑TATAGCTCGTCCGTTTATCGGTGACa                                                       >  1:78078/1‑250 (MQ=255)                                                                                                                                                                                            atcggCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATGCGGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACaaa                                                     <  1:607517/247‑1 (MQ=255)                                                                                                                                                                                                    ggTTACCTCAGTAAC‑GGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTc                                                                         <  1:27740/223‑1 (MQ=255)                                                                                                                                                                                                            cGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATa                                                                                                                                                     <  1:127671/140‑1 (MQ=255)                                                                                                                                                                                                                aaccggCCACTCTGCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTa                                  <  2:316872/246‑1 (MQ=255)                                                                                                                                                                                                                 aC‑TGCCACTCTTCCGGTACGGTCATTCTGGAACAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTGCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCTGTAACTTCCAGCGTAcc                                <  2:526310/251‑1 (MQ=255)                                                                                                                                                                                                                     gCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCgg                              <  1:319516/250‑1 (MQ=255)                                                                                                                                                                                                                      cccgacgctcttccgaactaTTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTa                            <  2:29512/232‑1 (MQ=255)                                                                                                                                                                                                                       cACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTa                            >  1:301911/1‑250 (MQ=255)                                                                                                                                                                                                                                 gTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGAc                 <  1:26666/251‑1 (MQ=255)                                                                                                                                                                                                                                   aCGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCt               <  1:301163/251‑1 (MQ=255)                                                                                                                                                                                                                                     ggTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAACTCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACAACCGTCGTGC‑TATCGCTCGTCCGGTTATCGGTGACCAAGCCGGTAACATCCCGCGTACCGGTCACCGTCACGAACTgg             >  2:251500/1‑251 (MQ=255)                                                                                                                                                                                                                                        cATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTg          <  1:130413/251‑1 (MQ=255)                                                                                                                                                                                                                                         aTTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTa                            >  1:29512/1‑232 (MQ=255)                                                                                                                                                                                                                                            cTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTgg      >  1:94773/1‑250 (MQ=255)                                                                                                                                                                                                                                                aTCAGCTGGGCGAAGAGCACATGAAAAACGGCAAGCCGATTTTCTATACCTCCGCTGACTACGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAAATCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCCGTCGGTT‑TATCGCTCGTCCCTTATTCGGTGACAAAACCGGTCACTTCCCACGTACCGGTAACCGtaccgaccgggaagtggcgacg  >  2:675156/1‑230 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        GGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCG  >  NC_000913/4619825‑4620313

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