Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,061,743	C→T	S298S (AGC→AGT)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,061,743	0	C	T	98.9%	291.9 / NA	88	S298S (AGC→AGT)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (47/40);  minor base A (1/0);  total (48/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTTCCGGCAACATGCCGCGATTTAATACCATCAGTATCAGCGGTTACCACATGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCA  >  NC_000913/3061497‑3061988                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      gtgtTCCGGCAACATGCCGCGATTTAATACCATCAGTATCAGCGGTTACCACATGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAg                                                                                                                                                                                                                                                   <  2:502208/251‑1 (MQ=255)  tgtTCCGGCAACATGCCGCGATTTAATACCATCAGTATCAGCGGTTACCACATGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGc                                                                                                                                                                                                                                                  >  1:118631/1‑251 (MQ=255)    ttCCAGCAACATGCCGCGATTTAATACCATCAGTATCAGCGGTTACCACATGGGGGAAGTGGGGGCCAACTGCGTGCAGCAGGTAGCATTTACGCCCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGAATGAATCTGTTTCTGAACGTCGCCA‑TGTTGGGTGCGGCACGTTATTGATGGAGTGAAACgg                                                                                                                                                                                                                                                >  2:23806/1‑251 (MQ=255)      ccGGCAACATGCCGCGATTTAATACCATCAGTATCAGCGGTTACCACATGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCCTCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTCGGGCACGTTATTTATGGAGTGAAGCGGTc                                                                                                                                                                                                                                              >  2:69869/1‑251 (MQ=255)       cGGCAACATGCCGCGATTTAATACCATCAGTATCAGCGGTTACCACATGGGTGAAGCGGGTTCCAACTGCGTGCAGCAGGTAGCATTTACGCTCTCTGATGGGATTGAGTACATCAAAGCAGCAATCTCAGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGAGCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCa                                                                                                                                                                                                                                             <  2:283434/251‑1 (MQ=255)       cGGCAACATGCCGCGATTTAATACCATCAGTATCAGCGGTTACCACATGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCa                                                                                                                                                                                                                                             <  2:224910/251‑1 (MQ=255)       cGGCAACATGCCGCGATTTAATACCATCAGTAGCAGCGTTTACCACATGGGTGAAGCGGGTGCCAACTGCGTTCAGCAGGTAGCATTTCCGCTCGCTGATGGGATTGAGTACCTCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCa                                                                                                                                                                                                                                             <  2:184026/251‑1 (MQ=255)            aCATGCAGCGATTAAATACCACCAGTATCAGCGGTTACCACCTGGGTGAAGCGCGCGCCAACTGCGTGCAGCAGGAACCCTTTACGCTCGATGATCGGATTGCGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGACGCCA‑TGTTCCCTGCGGCACGTTCTTTATCGAGTGCAGCGGTCAGNgg                                                                                                                                                                                                                                         >  2:606902/1‑250 (MQ=255)                           aTACCATCAGTATCAGCGGTTACCACATGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTTCGTGCCGCACGTTATTTATTGAGTTAAGTGGTCCGGGGGTTTGGCGCACCGGAc                                                                                                                                                                                                                         >  2:332587/1‑251 (MQ=255)                             aCCATCAGTTGCAGCGGTTACCACATGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCGCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAATCGGTCAGTGGATTTGGCGCACAGGAcc                                                                                                                                                                                                                        <  2:425791/250‑1 (MQ=255)                              ccATCAGTATCAGCGGTTACCACATGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCg                                                                                                                                                                                                                      <  2:522124/251‑1 (MQ=255)                                   aGNATCAGCGTTTACCACATGGGTGAAGCGGTTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGATTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGTACGTCGCCA‑TGTTTCGTGCGGCCCGTTTGTTATGGATTGAAGCGGTCAGAGTATTTGGCGCACAGCACCCGAAAt                                                                                                                                                                                                                  <  2:590362/250‑1 (MQ=255)                                             gTTACCACATGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTg                                                                                                                                                                                                       <  2:224371/251‑1 (MQ=255)                                                   acaTGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACCGGACCCGAAATCACTGGCGCTGCGTAcc                                                                                                                                                                                                 >  2:217588/1‑251 (MQ=255)                                                    caTGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAACTCACTGGCGCTGCGTAccc                                                                                                                                                                                                >  1:61134/1‑251 (MQ=255)                                                    caTGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTAccc                                                                                                                                                                                                <  1:95308/251‑1 (MQ=255)                                                         gTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATATGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGCGTGAAGCGGTCAGTGGATTTGGCGCACCGGACCCGAAATCACGGGCGCTGCGTACCCACTGc                                                                                                                                                                                           >  2:397219/1‑251 (MQ=255)                                                                               gcagcaGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctg                                                                                                                                                                      <  2:317005/250‑1 (MQ=255)                                                                                cagcagGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTTGCGCACCGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctgac                                                                                                                                                                    >  2:522748/1‑251 (MQ=255)                                                                                cagcagGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctgac                                                                                                                                                                    >  1:554741/1‑251 (MQ=255)                                                                                     ggTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGaa                                                                                                                                                                >  2:36747/1‑250 (MQ=255)                                                                                      gTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAAc                                                                                                                                                               >  2:374770/1‑250 (MQ=255)                                                                                          cATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTTGAGCCTGACTGAACAGGAt                                                                                                                                                          >  1:437844/1‑251 (MQ=255)                                                                                                 gctcgctGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGtata                                                                                                                                                   <  2:551925/251‑1 (MQ=255)                                                                                                  ctcgctGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGtata                                                                                                                                                   >  2:529100/1‑250 (MQ=255)                                                                                                     gctGATTGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaacaa                                                                                                                                               <  2:103672/251‑1 (MQ=255)                                                                                                        gATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaacaac                                                                                                                                              <  1:465583/249‑1 (MQ=255)                                                                                                         aTGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGTATGGATCTGTTTATGAACGTCGCCA‑TGTTGCTTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGTAGCCTGCCTGAACAGGATCCGCATTACAACGtt                                                                                                                                           >  2:464625/1‑251 (MQ=255)                                                                                                          ggggATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTa                                                                                                                                          <  2:25854/250‑1 (MQ=255)                                                                                                                    tACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGCTCCGTATAACAACGTTATCCGcaccac                                                                                                                                >  1:433267/1‑251 (MQ=255)                                                                                                                     aCATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCaccacc                                                                                                                               <  2:573814/251‑1 (MQ=255)                                                                                                                     aCATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCaccacc                                                                                                                               >  1:204658/1‑251 (MQ=255)                                                                                                                      cATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCTCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCATACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCAcctcca                                                                                                                              <  2:45816/251‑1 (MQ=255)                                                                                                                       aTCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCAt                                                                                                                             <  1:22349/251‑1 (MQ=255)                                                                                                                         cAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCAtt                                                                                                                            <  1:42649/250‑1 (MQ=255)                                                                                                                         cAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGGTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCAtt                                                                                                                            <  1:397219/250‑1 (MQ=255)                                                                                                                               agcaATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAgcgc                                                                                                                     >  2:444350/1‑251 (MQ=255)                                                                                                                                gcaATCTCTGCCGGACTGAAAATTGATGACTTCGCACCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAACTCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCCCCACCATTGACGCgct                                                                                                                    >  1:656834/1‑251 (MQ=255)                                                                                                                                 caATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGCACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGAAACGTATAAAAACGTTATCCGCACCACCATTGAAGCgctg                                                                                                                   >  1:439553/1‑251 (MQ=255)                                                                                                                                    tctcTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGAGCCGTATAACAACGTTATCCGCACCACCATTGAAGCgctggct                                                                                                                >  1:66693/1‑251 (MQ=255)                                                                                                                                          ccGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACg                                                                                                          >  1:301767/1‑251 (MQ=255)                                                                                                                                            ggACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGTTCAGTGGATTTGGCGCACAGGACCCGAAAGCACTGgc                                                                                                                                                                                                           <  2:94359/152‑1 (MQ=255)                                                                                                                                            ggACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGgc                                                                                                                                                                                                           >  1:94359/1‑152 (MQ=255)                                                                                                                                            cgACTGAAAATTCATGACTTCGCTGCTCGCCGGTCGTTCTTCTTCGGCATCGGCATTGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGTCACGTTATTTATTGAGTGAAGTGGTCAGTTGATGTTGCGCACAGGACCCGTTACCACTGGCGCTCCGTCCCCACTCCTAGACCTCAGGCTGGTGCCTGACTGAGCAGGATCCGTTTCACAACATAATCCGCACCGCGACTAAAGCGCTGTCTGCGACGCt                                                                                                        <  2:95923/250‑1 (MQ=255)                                                                                                                                                 gAAACTTGATTACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGAGTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGc                                                                                                    <  1:164539/250‑1 (MQ=255)                                                                                                                                                         atgaCTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCTCTGGCTGCGACGCTGGGCGGTACTCAg                                                                                           <  2:254734/251‑1 (MQ=255)                                                                                                                                                         atgaCTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAg                                                                                           <  1:175496/251‑1 (MQ=255)                                                                                                                                                           gaCTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTc                                                                                         >  1:71531/1‑251 (MQ=255)                                                                                                                                                            aCTTCGCCCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGAGGTACTCAGTc                                                                                         >  1:180961/1‑250 (MQ=255)                                                                                                                                                              ttCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACt                                                                                      >  2:240499/1‑251 (MQ=255)                                                                                                                                                                 gCTCCTCGCCTGTGATTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTC‑CCACTGTTGCGTGCGGCACGTTCTTTATGGCGAGCAGCGGTCAGTGGATTTGGCGCACAGGACCAGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATCACAACGTTATCCGCACCACCCTTGAAACGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGcc                                                                                   >  2:292938/1‑250 (MQ=255)                                                                                                                                                                  ctcctcGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCAt                                                                                  >  1:238676/1‑251 (MQ=255)                                                                                                                                                                    cctcGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGACCGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCt                                                                                                        <  2:326504/227‑1 (MQ=255)                                                                                                                                                                    cctcGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTAGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCt                                                                                                        >  1:326504/1‑227 (MQ=255)                                                                                                                                                                       cGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATAcc                                                                               <  1:628318/249‑1 (MQ=255)                                                                                                                                                                        gCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAAc                                                                            <  1:124848/251‑1 (MQ=255)                                                                                                                                                                                 tcttcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGa                                                                    <  1:217588/250‑1 (MQ=255)                                                                                                                                                                                       tcggcatcggcaTGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAgc                                                              <  1:190486/250‑1 (MQ=255)                                                                                                                                                                                         ggcatcggcatGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCt                                                           >  1:495091/1‑251 (MQ=255)                                                                                                                                                                                         ggcatcggcatGGATATGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCt                                                           >  2:54561/1‑251 (MQ=255)                                                                                                                                                                                           catcggcatGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCCTACCAACGCCTTTGACGAAGCGCTTg                                                         >  2:443228/1‑251 (MQ=255)                                                                                                                                                                                             tcggcatGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGt                                                       >  1:340074/1‑251 (MQ=255)                                                                                                                                                                                               ggcatGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGCCCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTgggtt                                                     >  2:415391/1‑248 (MQ=255)                                                                                                                                                                                                  atGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGcc                                                  >  1:238949/1‑251 (MQ=255)                                                                                                                                                                                                  atGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGACGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGcc                                                  >  1:445260/1‑251 (MQ=255)                                                                                                                                                                                                    ggATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCt                                                 >  1:222918/1‑250 (MQ=255)                                                                                                                                                                                                    ggATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAACCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCt                                                 <  2:373582/250‑1 (MQ=255)                                                                                                                                                                                                         tGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGAt                                           >  1:206548/1‑251 (MQ=255)                                                                                                                                                                                                          gTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCg                                             <  1:402279/248‑1 (MQ=255)                                                                                                                                                                                                            ttATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGAGCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGAttt                                         <  1:292938/250‑1 (MQ=255)                                                                                                                                                                                                             tATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGGAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATNtct                                       >  2:606093/1‑251 (MQ=255)                                                                                                                                                                                                               tGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGTTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTtctc                                      >  2:467973/1‑250 (MQ=255)                                                                                                                                                                                                                gAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCa                                     <  2:475593/250‑1 (MQ=255)                                                                                                                                                                                                                 aaCGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTTCTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTGGCCTACCGATTGCTCAGc                                   <  2:280284/251‑1 (MQ=255)                                                                                                                                                                                                                  aCGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGc                                   >  1:526117/1‑250 (MQ=255)                                                                                                                                                                                                                    gTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACg                                <  1:524479/251‑1 (MQ=255)                                                                                                                                                                                                                         cA‑TGTTGCGTGCGGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTg                           >  2:184360/1‑251 (MQ=255)                                                                                                                                                                                                                                 cgtgcgGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCa                     >  2:622977/1‑250 (MQ=255)                                                                                                                                                                                                                                    gcgGCACGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGAATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAacac                 >  1:18105/1‑251 (MQ=255)                                                                                                                                                                                                                                          cGTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACACGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTAACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACAGATTTCTCAGCACGCATTGCCCGCAACAccc               >  1:646229/1‑247 (MQ=255)                                                                                                                                                                                                                                           gTTATTTATGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACAccc               >  1:119771/1‑246 (MQ=255)                                                                                                                                                                                                                                                ttatGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGatcatca      <  2:301767/250‑1 (MQ=255)                                                                                                                                                                                                                                                 tatGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGatcatcat     <  1:37077/250‑1 (MQ=255)                                                                                                                                                                                                                                                 tatGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGatcatcat     <  1:332587/250‑1 (MQ=255)                                                                                                                                                                                                                                                 tatGGAGTGAAGCGGCCAGCGGCTTTGGCGCACAGGACCCGAAATCACTGCCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCGTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGatcatcatc    >  2:70862/1‑251 (MQ=255)                                                                                                                                                                                                                                                  atGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCAtt                                                                                                                            <  2:178076/130‑1 (MQ=255)                                                                                                                                                                                                                                                  atGGAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCAtt                                                                                                                            >  1:178076/1‑130 (MQ=255)                                                                                                                                                                                                                                                    ggAGTGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCa  <  1:241074/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      GTGTTCCGGCAACATGCCGCGATTTAATACCATCAGTATCAGCGGTTACCACATGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCA‑TGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCA  >  NC_000913/3061497‑3061988

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