Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,061,168	G→A	D107N (GAC→AAC)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,061,168	0	G	A	100.0%	174.8 / NA	52	D107N (GAC→AAC)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/25);  total (27/25)

GGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGTTTTATATCGTCGCCGCAGAAGAGCAAGGTGTTACACCTGATAAACTGACCGGCACCATTCAAAACGATA  >  NC_000913/3060924‑3061414                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       ggTTCATCGAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTCGCGTTTAACCTTg                                                                                                                                                                                                                                                  >  2:47150/1‑251 (MQ=255)           aCCGCGGAAGGGATCGCCATCAAGCCGCTGTATAACGACGCCGCTCTCGATCATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTCCGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTCTCGCCGTAACCTGGCCGCCGGGCCAAAAGGTCTTTCCGTTGCGTTTAACCTTTCCACCCCCCg                                                                                                                                                                                                                                        >  2:393814/1‑251 (MQ=255)                    gggATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACg                                                                                                                                                                                                                                >  1:183363/1‑250 (MQ=255)                        gcGCCATCAAGCCGCTTTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGGCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGTCTACGACTcc                                                                                                                                                                                                                           <  1:283464/250‑1 (MQ=255)                         cGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGAGTTTAACCTTGCCACCCACCGTGGCTCCGACTCCg                                                                                                                                                                                                                          >  2:139089/1‑251 (MQ=255)                                 aGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCcg                                                                                                                                                                                                                  <  1:147891/251‑1 (MQ=255)                                                   ccGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAAc                                                                                                                                                                                                                     <  1:159972/230‑1 (MQ=255)                                                   ccGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAAc                                                                                                                                                                                                                     >  2:159972/1‑230 (MQ=255)                                                     gATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACCCCGATACCCCGCGCGTGGCGGGCGACGTCgg                                                                                                                                                                                              >  2:203743/1‑251 (MQ=255)                                                                   ggagtggACAGGTACCCTTCCTGGTTTGCCGCCCTACGGTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCTCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGt                                                                                                                                                                                  <  1:48915/243‑1 (MQ=255)                                                                   ggAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTc                                                                                                                                                                                 <  1:212441/250‑1 (MQ=255)                                                                         gACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACATCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGTCAAAAAGGTCTTTCCGTTGCGTTGAACCTTGCCACCCACCGTGGCTACGACGCCGATAACCCGCGCGTGGCGGGCGACGTCGTCAAAGCGTGCGTCGCTATc                                                                                                                                                                           <  1:284498/250‑1 (MQ=255)                                                                             ggTACCCTTCCTGGTTTGCCGCCCTGCGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTGCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGacac                                                                                                                                                                      <  1:13476/251‑1 (MQ=255)                                                                                      ccTGGTTTTCCGCCCTACGTGCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGTCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAg                                                                                                                                                              <  1:394717/250‑1 (MQ=255)                                                                                         ggTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGata                                                                                                                                                           <  1:140579/250‑1 (MQ=255)                                                                                          gTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGGCGCTATCGACACCGTGGAAGATATg                                                                                                                                                         >  2:303012/1‑251 (MQ=255)                                                                                              gccgccCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTGTCCGTTGCGTTTAACCTTGCCACCCACCGGGGCTACGACTCCGAGAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGaaa                                                                                                                                                      <  1:393814/250‑1 (MQ=255)                                                                                                   ccTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGCCACCCTGGAAGATAAGAAAGGCCTg                                                                                                                                                >  2:400068/1‑251 (MQ=255)                                                                                                           cGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGAcc                                                                                                                                         <  1:178881/250‑1 (MQ=255)                                                                                                            gTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTACTGTTCGAcc                                                                                                                                         >  1:18500/1‑249 (MQ=255)                                                                                                                cccGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCaagtcc                                                                                                                                   >  2:330170/1‑246 (MQ=255)                                                                                                                                 ataCCGCCCAACCGTTGACCAACCGTCAGTATGCTGGTTTTTCAGCAGCAAAAGAGTCCCCCGCTTTTTATCGCCGTAACCTGGCCTCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAAGGTc                                                                                                                   <  1:239740/250‑1 (MQ=255)                                                                                                                                    ccGCCCAACCGGGGACCATCCGTCAGTATGCTGGTTTGTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTGCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGGCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGTATAAAAGGGCGGtt                                                                                                               <  2:398334/251‑1 (MQ=255)                                                                                                                                      gCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGAGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTc                                                                                                             >  1:276925/1‑251 (MQ=255)                                                                                                                                               tGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCCCCATGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGAAGTCGGCAAAGCGTGCGTAGCTATCGGCACCGTGCAAGATCTGACAGTCCTGTTCGGCCATATCCCGCTGGATAAAATGTCGGTTTCGATGAcaac                                                                                                    >  1:425455/1‑248 (MQ=255)                                                                                                                                                 gACCATCCGTCAGTCTGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGCTCCCGCTGGATCAAATGTCGGTTTCGATGACCATg                                                                                                   >  1:323287/1‑250 (MQ=255)                                                                                                                                                    cATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATCTGAAAGTCCTGTTCGACCAGATCCCGCTTGATAAAATGTCGGTTTCGATGACCATGAATg                                                                                               >  2:80447/1‑251 (MQ=255)                                                                                                                                                     aTCCGTCAGTATGCTGGTTTTTGAACAGCAAAAGAGTTCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTATCCTTGCCACCCACCGTGGCTACGACTCCGATATCCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTCTCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATgg                                                                                              <  1:265749/251‑1 (MQ=255)                                                                                                                                                       ccGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGgc                                                                                             <  1:197586/250‑1 (MQ=255)                                                                                                                                                               aTGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCGGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGAGCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTa                                                                                    <  1:31297/251‑1 (MQ=255)                                                                                                                                                                tGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCg                                                                                                                                            <  1:110737/194‑1 (MQ=255)                                                                                                                                                                tGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCg                                                                                                                                            >  2:110737/1‑194 (MQ=255)                                                                                                                                                                    ggTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGt                                                                               >  1:140780/1‑251 (MQ=255)                                                                                                                                                                          tCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGATGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTAACAGTACTGGc                                                                         >  1:87727/1‑251 (MQ=255)                                                                                                                                                                            aaCAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGt                                                                       <  2:288990/251‑1 (MQ=255)                                                                                                                                                                            aaCAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGACACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCAGCTGGATAAAATGTCGGGTTCGATGACCATGAATGGAGCCGTGCTACCAGTACTGGCGt                                                                       >  1:45690/1‑251 (MQ=255)                                                                                                                                                                             aCAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTTGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGtt                                                                      <  1:393362/251‑1 (MQ=255)                                                                                                                                                                                gCAAAAGCGTCCAACGCTTTTTATCGCCGTAACCTGGACGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGCCCATGAATGGCGCAGTGCTACCAGTACTGGCGttttt                                                                   >  2:390458/1‑250 (MQ=255)                                                                                                                                                                                 cAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCAACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGCCCCCGTGGAAGATATGAAAGTCCAGTTCGACCAGATCCCGCTGGATAAAATGGCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGTGTTTtat                                                                  >  2:279005/1‑251 (MQ=255)                                                                                                                                                                                      gagTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGGATGGCGCAGTGCTACCAGTACTGGCGTTTTATAtcgt                                                             >  2:136382/1‑251 (MQ=255)                                                                                                                                                                                          ccAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGTTTTATATCGTcgcc                                                         <  1:409071/251‑1 (MQ=255)                                                                                                                                                                                                     tATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGTTTTATATCGTCGCCGCAGAAGAGc                                               >  2:254683/1‑250 (MQ=255)                                                                                                                                                                                                                   ggCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGTTTTATATCGTCGCCGCAGAAGAGCAAGGTGTTACACCTg                                <  2:302737/251‑1 (MQ=255)                                                                                                                                                                                                                   ggCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGTTTTATATCGTCGCCGCAGAAGAGCAAGGTGTTACACCTg                                <  1:330170/251‑1 (MQ=255)                                                                                                                                                                                                                          gggCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGTTTTATAtcgtcg                                                           <  2:381400/217‑1 (MQ=255)                                                                                                                                                                                                                          gggCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGTTTTATAtcgtcg                                                           >  1:381400/1‑217 (MQ=255)                                                                                                                                                                                                                          gggCAAAAAGGTCTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGCGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGTTTTATATCGTCGCCGCAGAAGAGCAAGGTGTTACACCTGATAAACt                         >  1:70167/1‑251 (MQ=255)                                                                                                                                                                                                                                      cTTTCCGTTGCGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGTTTTATATCGTCGCCGCAGAAGAGCAAGGTGTTACACCTGATAAACTGACCGGCACCAt             >  1:287739/1‑251 (MQ=255)                                                                                                                                                                                                                                       tttCCGTTGCGTTTAACCATGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGACAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATATCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGTTTTATATCGTCGCACCCGAAGAGCAAGGTGTTACACCTGATAAACTGAACGTCACCAtt            >  2:86514/1‑251 (MQ=255)                                                                                                                                                                                                                                            gtgacGTTTAACCATGCCACCCGCCGTGGCTACGACTCCGAGAAACCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATAGACACCGTGGAAGATATGAAAGGCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGTTTTATATCGTCGCCGCAGAAGAGCAAGGTGTTACACCTGATAAACTGACCGGCACCATTCaaaa       <  2:230336/247‑1 (MQ=255)                                                                                                                                                                                                                                                cGTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGTTTTATATCGTCGCCGCAGAAGAGCAAGGTGTTACACCTGATAAACTGACCGGCACCATTCAAAACGat   >  2:365703/1‑251 (MQ=255)                                                                                                                                                                                                                                                 gTTTAACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGGTTTATATCGTCGCCGCAGAAGAGCAAGGTGTTACACCTGATAAACTGACCGGCACCATTCAAAACGata  <  2:70167/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       GGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAACCCGCGCGTGGCGGGCGACGTCGGCAAAGCGGGCGTCGCTATCGACACCGTGGAAGATATGAAAGTCCTGTTCGACCAGATCCCGCTGGATAAAATGTCGGTTTCGATGACCATGAATGGCGCAGTGCTACCAGTACTGGCGTTTTATATCGTCGCCGCAGAAGAGCAAGGTGTTACACCTGATAAACTGACCGGCACCATTCAAAACGATA  >  NC_000913/3060924‑3061414

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