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.6%	208.6 / NA	70	A356V (GCC→GTC)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (28/40);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGA  >  NC_000913/3061673‑3062156                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 tcttcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCAGGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATGTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGa                                                                                                                                                                                                                                            <  2:255894/250‑1 (MQ=255)         ggcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCCGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGACCAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGACTTTGACGAAGCGCt                                                                                                                                                                                                                                   >  1:513955/1‑251 (MQ=255)           catcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTg                                                                                                                                                                                                                                 <  1:58393/251‑1 (MQ=255)               ggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGttt                                                                                                                                                                                                                             <  1:183756/251‑1 (MQ=255)                    ggATCTGTTTATGCACGTCGCCATGTTAAGTACGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCAATGGCGCGGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTAGAACAACGTTATCCGCACCACCCTTGAAGCGCTGGATGCGCCGCTGGGAGGTACTCAGTCACTGCATACCCACGCTTATGACGACACACTTGGTTTGcca                                                                                                                                                                                                                         >  2:528857/1‑249 (MQ=255)                           tttATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGAttt                                                                                                                                                                                                                 >  1:255015/1‑251 (MQ=255)                             tATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTtct                                                                                                                                                                                                               >  1:157638/1‑251 (MQ=255)                                   cGTCGCCATTGTGCGTGCGGCACGGTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCCCTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTAGCCGCGCCTCCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACTGATTTCTCAGCAc                                                                                                                                                                                                         <  2:534582/251‑1 (MQ=255)                                      cGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGGTCACAGGACCCGAAATCACTGGCGCGGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCa                                                                                                                                                                                                      <  2:408895/251‑1 (MQ=255)                                          aTGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCCTACCAACGTCTTTGCCGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGc                                                                                                                                                                                                  >  2:434056/1‑251 (MQ=255)                                               gcgtgcgGCACGTTCTTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGCAATACCTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCCTTGAAGCGCTGGCTGCGACGCTGGGAGGTACTCAGTCACTGCATCCCAAAGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCCCGCATTGCCCGcc                                                                                                                                                                                             >  2:82766/1‑250 (MQ=255)                                                   gcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAaca                                                                                                                                                                                          <  1:434056/250‑1 (MQ=255)                                                       cACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACAccc                                                                                                                                                                                       >  1:341633/1‑249 (MQ=255)                                                           ttatttatGGAGCGAAGCGGTCAGTGGATTTGGCGCAAAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACATCAGGCTGGAGCCTGACTGAACAGGATACGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGAGGCGACGCTGAGCGGTACTAGGTCACTGCATACCATCGTCTTTGACGAAACGCTTGGCTTGCCTACCGATTTCGCAGCCCGCATTGCCCGCAACACCCAgctcc                                                                                                                                                                                 >  2:434331/1‑247 (MQ=255)                                                            tatttatGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGatcat                                                                                                                                                                                >  1:278972/1‑251 (MQ=255)                                                                   ggAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGGTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCa                                                                                                                                                                          <  1:8828/250‑1 (MQ=255)                                                                    gAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGCCGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCgagatcggaagcgcgtcgtgcggg                                                                                                                                                                        >  2:490823/1‑228 (MQ=255)                                                                    gAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCg                                                                                                                                                                                               <  1:490823/228‑1 (MQ=255)                                                                      gCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGAGGCGTAACCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGgaa                                                                                                                                                                      <  2:57078/251‑1 (MQ=255)                                                                           gCGGGCAGTGGATGTGGCGCACAGGCCCCGAAATCACTTGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCAGCCAGGAAGAATc                                                                                                                                                                 <  2:142524/251‑1 (MQ=255)                                                                              gTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGAttt                                                                                                                                                                                                                 <  2:486055/200‑1 (MQ=255)                                                                              gTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGAttt                                                                                                                                                                                                                 >  1:486055/1‑200 (MQ=255)                                                                              gTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGa                                                                                                                                                              <  1:556344/251‑1 (MQ=255)                                                                              gTCAGTGGATATGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGAATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGATGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACCCCCAGATCATCATCCAGgaagaaaa                                                                                                                                                                 >  2:344008/1‑246 (MQ=255)                                                                                cAGTTGATTTGGCGCGCAGGACCCGAAATCCTTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACGCAGTCACTGCATACCAACGTCTGTGACGAAGCGCTTGGTTTGCGTACCGATGTCTAAGGACGCATTGCGCGCAACACCCAGTTCATCATCCGGGAGGAATCAGaa                                                                                                                                                             <  2:245149/250‑1 (MQ=255)                                                                                   tGGATTTGGCGCACAGAACCCGAAATAACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGCACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCAATGCATACCAACGTCTTTGACGAAGAGCTTGGTTTGCCTACCCATTTCTCAGCACGCATTGCCCGCCACACCCAGCTCATCCTCCCGGAGGAATCAGAActat                                                                                                                                                         >  2:167693/1‑249 (MQ=255)                                                                                                   gACCCGAAATCACTGGCGCTGCGTACCCACTTCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTTGGCTGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATcc                                                                                                                                         <  2:130532/251‑1 (MQ=255)                                                                                                      ccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTCCCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGACCTCTGCCGCACCGTCGATCCACt                                                                                                                                      >  2:315449/1‑251 (MQ=255)                                                                                                         aaaTCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGc                                                                                                                                   <  1:351909/251‑1 (MQ=255)                                                                                                          aaTCACTGGCGCTGCGTCCCCACGGCCAGACCTCCGGCCGGTGCCTGACTGAACAGGTTCCGGATAATAACGTTATCCGCACCGCCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGTTCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGcc                                                                                                                                  <  2:68042/251‑1 (MQ=255)                                                                                                            tCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCgg                                                                                                                                >  1:542104/1‑251 (MQ=255)                                                                                                             ccctggacccgCGTACCCACTGCCCGTCCGCAGGCTGGAGACTGACTGAACAGGATCCGGATAACAACTTGATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGTCTTGCCCGCAACACCCAGAGCATCATCCAGGAAAAATCATAACTCTGCCGCACCGTCGATCCACTGGCCGGa                                                                                                                               <  2:345068/241‑1 (MQ=255)                                                                                                               caggCGCTGCGTGCCCACGGTCATAACTAAGGCCGGAGCCTGACTGAACAGGATCCGTAGACCAACGTTATCCGCACCACCATTGAAGCGCTGGCTGTGACGCTGGGCGGTACTCAGTTACTGCAGACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCGCAGCACGCATTGCCCGCAACCCCCAGATCATCATCCAGTAAGAATCAGAACGCTGCCGCATCGTCGATCCACTGTCCGGATc                                                                                                                             <  2:377128/249‑1 (MQ=255)                                                                                                                     cTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGCTCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGACTACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCCCTGGCCGGATCatata                                                                                                                        >  2:410003/1‑249 (MQ=255)                                                                                                                       gCGTCCCCACGGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGACCCGTATAACAACGTTATCCGCACCCACATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCATGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACa                                                                                                                     <  2:58404/251‑1 (MQ=255)                                                                                                                            cccACTGCCAGACCGCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAg                                                                                                                <  1:302971/251‑1 (MQ=255)                                                                                                                             ccACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCCCTGGCCGGATCCTATTACATTGAGt                                                                                                               >  2:235317/1‑251 (MQ=255)                                                                                                                                   ccAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTg                                                                                                          >  1:550618/1‑250 (MQ=255)                                                                                                                                     aGACCTAAGGCTGGAGCCTGACTGAACAGGATCCGTATAACATCGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGAcc                                                                                                       <  1:239362/251‑1 (MQ=255)                                                                                                                                      gACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCg                                                                                                      <  1:233704/251‑1 (MQ=255)                                                                                                                                              gCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATc                                                                                              <  1:118652/251‑1 (MQ=255)                                                                                                                                                tGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGcc                                                                                                                                  >  1:131598/1‑213 (MQ=255)                                                                                                                                                   aGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCaa                                                                                         <  1:57386/251‑1 (MQ=255)                                                                                                                                                       ggactgaACAGGATCCGTATACCAACGGTATCCGCACCATCATTGAAGCGGTGGCTGCGACGCTGGGCGGTACGCAGTTCCTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCTATGTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGAACCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACaa                                                                                     <  2:57521/250‑1 (MQ=255)                                                                                                                                                        gactgaACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACa                                                                                      <  1:527950/249‑1 (MQ=255)                                                                                                                                                            gaACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGcc                                                                                  <  1:402504/249‑1 (MQ=255)                                                                                                                                                              aCAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCag                                                                                >  1:430236/1‑249 (MQ=255)                                                                                                                                                                        tataACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACCTTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATcc                                                                     >  2:146614/1‑250 (MQ=255)                                                                                                                                                                         ataACAACGTTATCCGCCCCACCATGGAAGCGCGGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCa                                                                    <  2:157638/250‑1 (MQ=255)                                                                                                                                                                          taACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGCTCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGcaatgatacaac                                                                  >  2:289687/1‑240 (MQ=255)                                                                                                                                                                           aacaacGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCCTCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCacc                                                                  >  2:150067/1‑248 (MQ=255)                                                                                                                                                                           aacaacGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAAc                                                                  <  1:471136/250‑1 (MQ=255)                                                                                                                                                                           aacaacGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAAc                                                                  >  2:127701/1‑250 (MQ=255)                                                                                                                                                                              cacGTTACCTGCACCATCATTGCAGCGCTGGCTGCGACGCTGGGTGGTCCTCAGTCTCTGCCTACCTACGTCTTTTACGAATCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGACCTCTGCCGCCCCGTCGATCCACTGGCCGAATCCTATTTCATTGAGTCGCTTACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGAt                                                              <  1:182379/250‑1 (MQ=255)                                                                                                                                                                                cGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATcg                                                            <  1:242152/251‑1 (MQ=255)                                                                                                                                                                                 gTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATcga                                                           <  1:23360/251‑1 (MQ=255)                                                                                                                                                                                             cccactGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTgg                                               <  1:454848/246‑1 (MQ=255)                                                                                                                                                                                                aTTGAAGCGCTGGCTGCGACGCGGGGAGGTACTCAGTCAATGCATACCAACGTCTTTGACGAAGCGCTTGGATTGCCTACCGATTTCTCAGCACGCATTGACCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGACGCACCGTCGATCCACTGACCGGATCCTATTCCATTGGGTCGCTGACCGATACAATCGTCAAAAAAGCCAGAGCTATTATCAAACAGATCGACGACGCCGGTggga                                             >  2:82373/1‑248 (MQ=255)                                                                                                                                                                                                        gcgtgcggcgccGCTGGGCGGTACGCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTGTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACACAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGGCGATCCACTGGCCGGATCCTATTACATTGAGGCGCTGACCGATCAAATCGGCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAg                                    <  2:73792/240‑1 (MQ=255)                                                                                                                                                                                                                     cTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCaa                                                                                         >  2:339164/1‑185 (MQ=255)                                                                                                                                                                                                                     cTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCaa                                                                                         <  1:339164/185‑1 (MQ=255)                                                                                                                                                                                                                       gggCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGt                      <  1:95639/250‑1 (MQ=255)                                                                                                                                                                                                                           ggTACTCAGTCACTGCATACCAACGTCTTTGACGAAACGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCATATTAACTTGCGTAGCTGACCGATCAAATCGACAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGGGGCAAGGCGAAATCGATCGAAGCAgggaagcc                 >  2:411860/1‑245 (MQ=255)                                                                                                                                                                                                                            gTACCCAGTCACTGCAAACTAACATCTTTGACTAAGCGCTTGGTTTGGCTACCGATTTCGCAGCACGCATCGCCCGCAACACCCAAATCATCATCCAGGAAGAATCAGAACTCTTCCGCACCGTCGATCCACTGGCCGGATCCGATTACATTGAATCGCTGACCGATCAAAGCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAATCAGGTCTGCCa                <  2:47602/251‑1 (MQ=255)                                                                                                                                                                                                                                 cAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAAc            <  1:541393/250‑1 (MQ=255)                                                                                                                                                                                                                                 cAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACg           >  1:291462/1‑251 (MQ=255)                                                                                                                                                                                                                                    tCACTGCATACACACGTCTTTGACGAAGCGATTGGTTTGCCTACCGCTTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGCCGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAACTCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGACGCCGGTGGCCTGGCGAAAGCGACCGCAGCAGGCCGGCCAATACgcat        >  2:54434/1‑248 (MQ=255)                                                                                                                                                                                                                                      aCTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGa      <  1:230791/251‑1 (MQ=255)                                                                                                                                                                                                                                        tGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATc    >  1:445921/1‑251 (MQ=255)                                                                                                                                                                                                                                          cATACCAACGTCTTTGGCGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGa  <  1:243742/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 TCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGA  >  NC_000913/3061673‑3062156

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