Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,297,668	C→T	D567N (GAC→AAC)	fdhF ←	formate dehydrogenase‑H, selenopolypeptide subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,297,668	0	C	T	100.0%	230.5 / NA	71	D567N (GAC→AAC)	fdhF	formate dehydrogenase‑H, selenopolypeptide subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/35);  total (36/35)

TGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCC  >  NC_000913/4297543‑4297798                                                                                                                              |                                                                                                                                     tgCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTAcc                                                                                                                         >  2:676378/1‑139 (MQ=255)      aGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑c                                                                                                                  >  2:1034221/1‑139 (MQ=255)        ccAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGc                                                                                                                >  2:1272133/1‑139 (MQ=255)           ggTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAgg                                                                                                             >  1:63210/1‑139 (MQ=255)            gTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGAt                                                                                                                                   >  2:1759757/1‑118 (MQ=255)            gTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGAt                                                                                                                                   <  1:1759757/118‑1 (MQ=255)                 tcaGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTgaaga                                                                                                       >  1:818754/1‑139 (MQ=255)                  cagccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTgaagaa                                                                                                      <  2:818754/139‑1 (MQ=255)                  cagccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTgaagaa                                                                                                      >  1:720952/1‑139 (MQ=255)                  cagccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTgaagaa                                                                                                      <  1:765840/139‑1 (MQ=255)                      cagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCa                                                                                                                       >  1:1856932/1‑120 (MQ=255)                      cagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCa                                                                                                                       <  2:1856932/120‑1 (MQ=255)                       agcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTgcg                                                                                                 <  2:63210/139‑1 (MQ=255)                          gccgccAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCa                                                                                              >  1:1607496/1‑139 (MQ=255)                            cgccAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGa                                                                                            <  2:1453641/139‑1 (MQ=255)                                   gCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                                     <  2:145345/139‑1 (MQ=255)                                       cacaGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGa                                                                                            <  1:827397/128‑1 (MQ=255)                                       cacaGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGa                                                                                            >  2:827397/1‑128 (MQ=255)                                        acaGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTa                                                                                >  2:715894/1‑139 (MQ=255)                                               ccGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                         <  1:1034221/139‑1 (MQ=255)                                               ccGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                         <  2:1607496/139‑1 (MQ=255)                                                    cATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCGCAGGTGAAGAACTGCGCCAGACCGTTCgggg                                                                                  >  2:1167431/1‑127 (MQ=255)                                                    cATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGTCGCGCAGGTGAAGAACTGCGCCAGACCGTTCgggg                                                                                  <  1:1167431/127‑1 (MQ=255)                                                      tCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCttt                                                                  >  2:481412/1‑139 (MQ=255)                                                         aaCGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTaaa                                                               <  2:720952/139‑1 (MQ=255)                                                                agagTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATaa                                                        >  1:1705032/1‑139 (MQ=255)                                                                 gagTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATaag                                                       <  1:715894/139‑1 (MQ=255)                                                                   gtagtGACCAACTTCACGCACCGTTGATAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATaagaa                                                     <  1:3878/139‑1 (MQ=255)                                                                       tGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                         <  1:1032903/115‑1 (MQ=255)                                                                       tGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACt                                                                         >  2:1032903/1‑115 (MQ=255)                                                                            aaCTTCACGCACCGTTGACAGTACCATCGGTTACTCGTCGGTGAGTTTGTTGATCGGCGCTACCCAGC‑‑AGCAGGTGAAGAAATGCGTCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACGGATAAGAAGTCCCCTGa                                            >  1:664430/1‑139 (MQ=255)                                                                            aaCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTgaagaa                                                                                                      <  1:628604/81‑1 (MQ=255)                                                                            aaCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTgaagaa                                                                                                      >  2:628604/1‑81 (MQ=255)                                                                            aaCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                                     >  2:1676907/1‑98 (MQ=255)                                                                            aaCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                                     <  1:256011/98‑1 (MQ=255)                                                                            aaCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                                     >  2:256011/1‑98 (MQ=255)                                                                            aaCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCg                                                                                     <  1:1676907/98‑1 (MQ=255)                                                                            aaCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGa                                            >  2:1795833/1‑139 (MQ=255)                                                                              cTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATc                                          >  2:686126/1‑139 (MQ=255)                                                                                 cACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGc                                       <  1:1795833/139‑1 (MQ=255)                                                                                 cACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGc                                       <  1:1272133/139‑1 (MQ=255)                                                                                             aCAGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATc                           >  2:643426/1‑139 (MQ=255)                                                                                               aGTACCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATcgc                         <  1:686126/139‑1 (MQ=255)                                                                                                  aCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATc                                          >  1:1806265/1‑119 (MQ=255)                                                                                                  aCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATc                                          <  2:1806265/119‑1 (MQ=255)                                                                                                  aCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGc                      >  2:460147/1‑139 (MQ=255)                                                                                                  aCCATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGc                      <  2:1786769/139‑1 (MQ=255)                                                                                                    cATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAg                                                            >  2:1274575/1‑99 (MQ=255)                                                                                                    cATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGAt                                                          <  2:448944/101‑1 (MQ=255)                                                                                                    cATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGAt                                                          >  1:448944/1‑101 (MQ=255)                                                                                                    cATCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGCGGTATCAAACTTCTCTTTAAACAg                                                            <  1:1274575/99‑1 (MQ=255)                                                                                                      tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAg                                                                                             >  1:1612858/1‑64 (MQ=255)                                                                                                      tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAg                                                                                             <  2:1612858/64‑1 (MQ=255)                                                                                                      tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAgg                  <  2:369034/139‑1 (MQ=255)                                                                                                      tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAgg                  >  2:1531730/1‑139 (MQ=255)                                                                                                      tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAgg                  >  1:685261/1‑139 (MQ=255)                                                                                                      tCGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAgg                  <  1:643426/139‑1 (MQ=255)                                                                                                       cGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATaagaa                                                     <  2:1569967/103‑1 (MQ=255)                                                                                                       cGGGTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATaagaa                                                     >  1:1569967/1‑103 (MQ=255)                                                                                                         ggTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCTAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGc                                       <  1:1241801/115‑1 (MQ=255)                                                                                                         ggTACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTGTAAACAGATAAGAAGTCCCCTGATCGGc                                       >  2:1241801/1‑115 (MQ=255)                                                                                                           tACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGa                                                                                            <  2:264860/60‑1 (MQ=255)                                                                                                           tACTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGa                                                                                            >  1:264860/1‑60 (MQ=255)                                                                                                             cTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTga                                                                                                          >  1:1278855/1‑44 (MQ=255)                                                                                                             cTCGTCGGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTga                                                                                                          <  2:1278855/44‑1 (MQ=255)                                                                                                                gtcgGTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACtgaatg        >  1:447611/1‑139 (MQ=255)                                                                                                                    gTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCgg                       >  2:1056634/1‑120 (MQ=255)                                                                                                                    gTGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCgg                       <  1:1056634/120‑1 (MQ=255)                                                                                                                     tGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTgaagaa                                                                                                      >  2:166051/1‑40 (MQ=38)                                                                                                                     tGAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTgaagaa                                                                                                      <  1:166051/40‑1 (MQ=38)                                                                                                                      gAGTTTGTTGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGccc  <  2:685261/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                     TGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGT‑‑CGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCC  >  NC_000913/4297543‑4297798

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