Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,061,916	C→T	A356V (GCC→GTC)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,061,916	0	C	T	98.8%	261.8 / NA	84	A356V (GCC→GTC)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (42/41);  minor base A (1/0);  total (43/41)
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

TCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGC  >  NC_000913/3061676‑3062162                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        tcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGAc                                                                                                                                                                                                                                                  <  1:83550/248‑1 (MQ=255)      ggcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAACTCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCt                                                                                                                                                                                                                                          >  1:523890/1‑251 (MQ=255)      ggcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCt                                                                                                                                                                                                                                          >  1:133336/1‑251 (MQ=255)              catGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGc                                                                                                                                                                                                                                  >  1:525059/1‑251 (MQ=255)                    tCTGTTTATGAACGTCGCCATGTTGCGTGCGGAACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACAATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCg                                                                                                                                                                                                                            >  2:587986/1‑251 (MQ=255)                           aTGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTtctc                                                                                                                                                                                                                     >  1:226830/1‑251 (MQ=255)                           aTGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTGGCCTACCGATTtct                                                                                                                                                                                                                      <  2:244688/250‑1 (MQ=255)                                 gTCGCCAGGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCTAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACg                                                                                                                                                                                                               <  2:439831/251‑1 (MQ=255)                                    gCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCAt                                                                                                                                                                                                            >  1:540380/1‑251 (MQ=255)                                       aTGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGAACTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTTACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGc                                                                                                                                                                                                         <  1:239324/251‑1 (MQ=255)                                              gtgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAac                                                                                                                                                                                                  >  1:141522/1‑251 (MQ=255)                                                gcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAacac                                                                                                                                                                                                <  1:95531/251‑1 (MQ=255)                                                          atttatGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGatca‑t                                                                                                                                                                                      <  1:61852/250‑1 (MQ=255)                                                               tGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATcc                                                                                                                                                                                 >  1:127501/1‑250 (MQ=255)                                                                  aGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAgg                                                                                                                                                                              >  1:113037/1‑250 (MQ=255)                                                                      aaGCGGTCAGTGGAGTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCAGTGCCCGCAACACCCAGATCA‑TCATCCAGgaagaa                                                                                                                                                                         <  2:266598/251‑1 (MQ=255)                                                                         cgttCAGTGGATGTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑GCATCCAGGAAGAATCa                                                                                                                                                                      <  2:441673/248‑1 (MQ=255)                                                                          ggTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACATGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAg                                                                                                                                                                     <  2:64789/251‑1 (MQ=255)                                                                           gTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCa                                                                                                                                                                      <  1:238664/249‑1 (MQ=255)                                                                             cAGTGGCTTTCGCGCACCGGACCCGACATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACCACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAAc                                                                                                                                                                  >  2:78940/1‑251 (MQ=255)                                                                             cAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGaa                                                                                                                                                                   <  2:549913/250‑1 (MQ=255)                                                                                 ggATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGCAGAATCAGAACTCTg                                                                                                                                                              >  1:209929/1‑251 (MQ=255)                                                                                      cggCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCAAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGc                                                                                                                                                          >  2:278974/2‑250 (MQ=255)                                                                                             caGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGa                                                                                                                                                  <  1:230928/251‑1 (MQ=255)                                                                                              aGGACCCGAAATCACGGTCGCTGTGGACCCACTCCCAGACCTCAGGCTGGAGCCTGACGGAACAGGATCCGTATCACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGCGCGGTACTCAGCCACTGCATACCGACGACTTTGACGAAGCGCTCGGTTTGCCTACCGATTTCTCAGCACGCATTGCAACCACCACCCAGA‑CACTCATCCAGCAAGAATCAGAACCCTGCCGccccccata                                                                                                                                                  >  2:350063/1‑242 (MQ=255)                                                                                               ggACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGCAGAATCAGAACTCTGCCGCACCGTCGATc                                                                                                                                                >  1:276681/1‑251 (MQ=255)                                                                                                  cccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCAc                                                                                                                                             >  2:127115/1‑251 (MQ=255)                                                                                                  cccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGACTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCAc                                                                                                                                             <  1:78285/251‑1 (MQ=255)                                                                                                    cGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTg                                                                                                                                           <  2:372666/251‑1 (MQ=255)                                                                                                    cGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTg                                                                                                                                           <  1:241749/251‑1 (MQ=255)                                                                                                                cgcTGCATACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGAACCACTGGCCGGCTCCTAt                                                                                                                               >  2:568128/1‑251 (MQ=255)                                                                                                                 gcTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTAtt                                                                                                                              <  2:340440/251‑1 (MQ=255)                                                                                                                      gTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGCAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACAtt                                                                                                                         >  2:365345/1‑251 (MQ=255)                                                                                                                      gTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCGGCCGCACCGGCGATCCACTGGCCGGATCCTATTACAtt                                                                                                                         <  2:113037/251‑1 (MQ=255)                                                                                                                      gTACCCACTGACAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACTTTATCCGCACCACCAGTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCATTCTCTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCATACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGTATCCTATTACAtt                                                                                                                         <  2:431079/251‑1 (MQ=255)                                                                                                                                 cAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGAc                                                                                                              <  1:245646/251‑1 (MQ=255)                                                                                                                                  aGACCTCAGGCTGGAGCCTGACTTAACAGGAGCCGGATAACAACGTTTGCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCGCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTTGTACATTTAGCCGCTGAcc                                                                                                             <  2:525059/251‑1 (MQ=255)                                                                                                                                   gACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCGTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCg                                                                                                            <  1:586400/251‑1 (MQ=255)                                                                                                                                           gCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATc                                                                                                    <  1:10903/251‑1 (MQ=255)                                                                                                                                           gCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGGCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGTCCGATCAAATc                                                                                                    <  1:415169/251‑1 (MQ=255)                                                                                                                                                ggCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCaa                                                                                               <  2:193216/250‑1 (MQ=255)                                                                                                                                                aGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCa                                                                                                >  1:516557/1‑250 (MQ=255)                                                                                                                                                 gCCTGACTGAGCAGGAGCCGTATAACAACGTTAGCCGCACCACCATTGAAGCGCTGGCTGCTACGCTGGGCGGTACTCAGTCACAGCATACCAACGTCTTGGACGAAGCGCTTGGTTTGCCTACAGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAATAATCAGAACTCTGCCGCACCGTCTGTCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAAACGTCaaa                                                                                              <  2:278162/251‑1 (MQ=255)                                                                                                                                                     gactgaACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACaa                                                                                           <  1:440019/250‑1 (MQ=255)                                                                                                                                                      actgaAGAGGATCCGTATCCCAACGTTAACCCCAACACCCCTGGCGCGCTAGTTGCGACGCTGGGCGAAACTCAGTCACTGCCTACCACCGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTAGCCGGATCCTATTACATTGAGTAGCTGCCCGATCAAATCGTCaacctagc                                                                                         >  2:68383/1‑245 (MQ=255)                                                                                                                                                          aaCAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGcc                                                                                        >  1:137609/1‑248 (MQ=255)                                                                                                                                                           aCAGGATCCGTATAACAACGTTATCCGCACCACAATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCag                                                                                      >  1:454176/1‑249 (MQ=255)                                                                                                                                                                  ccGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCGTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACGGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCtatta                                                                              <  2:273115/250‑1 (MQ=255)                                                                                                                                                                  ccGTATAACAACGTTATCCGCACCAACATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAACCAAGCCAGAGCtattat                                                                             >  2:285904/1‑251 (MQ=255)                                                                                                                                                                     tataACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACTCTAGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGACGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCC‑TCATCCAGGAAGAATCAGAACTCTGCCACACCGTCGATCCACGGGCCGGATCCTATTACATTGAGACGCTGACCGAACACCTCGTCAAACAAGCCAAAGCTATTATcc                                                                           >  2:132095/1‑250 (MQ=255)                                                                                                                                                                     tataACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGCCGAAGCGCTTGGATTGCATACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGACAGAGCTATTATcc                                                                           >  1:363415/1‑250 (MQ=255)                                                                                                                                                                     tataACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGCTTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATcc                                                                           >  1:521609/1‑250 (MQ=255)                                                                                                                                                                         ccaacGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCGGGGCGGTACTCAGTCACTGCATTCCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTGCTCAGCACGCATTGCCCGCAACGCCCAGATCA‑TCATCCAGGAAGAATCAGAACTCGGCCGCACCGTCGATCCACGGGCCGGAACCTATTACAGTGAATCGCTTACCGGTCAGATCTTCAAACAAGCCAGAGCTATGATCCAACAg                                                                      <  2:552737/250‑1 (MQ=255)                                                                                                                                                                           aacGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACGCCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGTCCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGAt                                                                    <  1:42255/251‑1 (MQ=255)                                                                                                                                                                            acGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATc                                                                   <  2:141522/251‑1 (MQ=255)                                                                                                                                                                            acGTTAGCCGCACCGCCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCACCGTCTTTGACGAAGCACTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGACGCATCAGAACTCTGCCGCACCGCCGATCCACAGGCCGGATCATATAACATTGAGTCGCGGACCGATCAACTCGTCAAACGCGCAAGAGCTATTATCCAAc                                                                        >  2:94892/1‑246 (MQ=255)                                                                                                                                                                             cGTTATCCGCACCACCATTGAAGAGCTGGCTGCGACGCTGGGAGATACTCAGGCACTGCATACCAACGTCTTTGACGAAACGCTTGGTTTGCCTACCGATTTATCAGCACGCATTGCCCGCAAAACCCAGATCA‑TCATCCAGGAAGACTCAGAACTCTGCCGCACCGTCAATCCACTTGCAGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAAAAAACCAGAGCTATTATCCAACCGATcg                                                                  >  1:501450/1‑251 (MQ=255)                                                                                                                                                                                 aTCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATcgacga                                                              <  1:103746/251‑1 (MQ=255)                                                                                                                                                                                   ccGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATcgacga                                                              >  1:165556/1‑249 (MQ=255)                                                                                                                                                                                   ccGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCAGCACCGTCGATCCACTGGCCGGATCCTATTAAATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATcgacga                                                              >  1:544032/1‑249 (MQ=255)                                                                                                                                                                                    cGCACCCCCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTGTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGc                                                           <  2:331143/251‑1 (MQ=255)                                                                                                                                                                                    cGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGc                                                           >  1:475331/1‑251 (MQ=255)                                                                                                                                                                                          accaTTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTgg                                                     >  2:406158/1‑251 (MQ=255)                                                                                                                                                                                            caTTGAAGCGCTGGCTGCGACGCTGGGCGGAACTCAGTCACTGCATACCAAAGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGc                                                    >  2:378090/1‑250 (MQ=255)                                                                                                                                                                                               tGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACCTCTTTGACGAAGCGCTTGGCATGACTAACGACTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCAACGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGAAGAAGCCGGCGGCATgg                                                >  1:538451/1‑251 (MQ=255)                                                                                                                                                                                                 aaGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCAGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCg                                              >  2:198283/1‑251 (MQ=255)                                                                                                                                                                                                   gcgcTGGCTGCGACGCTGGGAGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCAAACAGATCGACGAAGCCGGTGGCATGGCGa                                             >  2:565712/1‑250 (MQ=255)                                                                                                                                                                                                           tgCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCCGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCg                                    >  1:408181/1‑251 (MQ=255)                                                                                                                                                                                                            gCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATc                                     <  2:226830/249‑1 (MQ=255)                                                                                                                                                                                                            gCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCg                                    <  1:565712/250‑1 (MQ=255)                                                                                                                                                                                                                    gggCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAgg                             <  2:408181/249‑1 (MQ=255)                                                                                                                                                                                                                     ggCGGTACTCAGTCACTGCCTACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGAACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATGGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATGATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCt                          <  2:257365/251‑1 (MQ=255)                                                                                                                                                                                                                        ggTACTCAGTCACTGCATACTAACGTCTTTGACGAATCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACGCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATGGAGGCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCTGGTGGCATGGCGAAATCGATCGAAGCAGGTCTGcc                       <  2:363512/251‑1 (MQ=255)                                                                                                                                                                                                                               aGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGCAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACg                 >  1:42549/1‑250 (MQ=255)                                                                                                                                                                                                                                  cACTGCATACCAACGTCTTTGACGAAGCTCTTGGTTTGCCTACCGATTGCTCAGCACGCATTGCCCGCAACACCCAGTTCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTAGTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATg             <  1:354287/251‑1 (MQ=255)                                                                                                                                                                                                                                  cACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAAt              <  1:283682/250‑1 (MQ=255)                                                                                                                                                                                                                                     tGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGACTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGACGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCACACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATc          >  2:567473/1‑251 (MQ=255)                                                                                                                                                                                                                                      gCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGAt           <  2:73451/249‑1 (MQ=255)                                                                                                                                                                                                                                        aTACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGa        >  1:273760/1‑250 (MQ=255)                                                                                                                                                                                                                                        aTACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCACGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGAGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGa        >  1:565432/1‑250 (MQ=255)                                                                                                                                                                                                                                         tACCAACGTCTTTGACGAAGCGCTTGGTTTACCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGTAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCTGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCACACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGGGACATGGCGCAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGa        >  1:532500/1‑249 (MQ=255)                                                                                                                                                                                                                                            cAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATATCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCATGTCTGCCAAAACGAATGATCGAAGAgg   <  2:348605/251‑1 (MQ=255)                                                                                                                                                                                                                                             aaCGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGACGAATCAGACCTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGc  >  1:339617/1‑251 (MQ=255)                                                                                                                                                                                                                                             aaCGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGc  >  1:228660/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                        TCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCA‑TCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGC  >  NC_000913/3061676‑3062162

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