Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,786,704	A→C	E482A (GAA→GCA)	gpmM →	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,786,704	0	A	C	100.0%	170.6 / NA	49	E482A (GAA→GCA)	gpmM	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (23/26);  total (23/26)

AGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATG  >  NC_000913/3786474‑3786943                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                aGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGGTGAAGAACTGGATCACTGCGTGGGAGAAGGCGCGGAAGCGGTTGAGTCCGTGGGTGGACAAATGCTGATCACCCCTGACCACGGTAACCCTGAACAGATGCGCGGTCCGGCAACGGGGCAGGCACACAAGGCACACACCCAACTGCCAGGTCCCCTGATTTACGTTGGGGGTAAGAACGAGAAAGCGGGTGCAGGCGGGAAACTTTCttaa                                                                                                                                                                                                                              >  1:433377/1‑247 (MQ=255)   gTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCACTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATc                                                                                                                                                                                                                           <  1:424103/251‑1 (MQ=255)        acCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCACTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGcc                                                                                                                                                                                                                      <  1:221403/251‑1 (MQ=255)                gATGGAAGCGGCGGTTAAAGCGGTTGAAGCACTGGATCACTGCGTGGGAGAAATCGCGAAAGGGGTTGAGTCCGTGGGGGGACAACTGCTGGTCACCCCCGGCCACCGGAAAGCTGGACAGATGCGGGATCCGGCAACGGGGCAGGCACACACGGCCCACACCAACCCGCCCGTTCCGCCGATTTACGGTGGGGGTAAGAACCTGAAAGCGGTTGCAGGCGGGCAACTTTTTGACATccgccccgccctgg                                                                                                                                                                                                              >  2:353576/1‑238 (MQ=255)                     aaGCGGCGGTTAAAGCGGTTGAAGCACTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCg                                                                                                                                                                                                         <  1:353576/251‑1 (MQ=255)                       gcggcgGTTAAAGCGGTTGAAGCACTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCt                                                                                                                                                                                                       <  1:350135/251‑1 (MQ=255)                         ggcggTTACAGCGGTTGCAGCACTGGATCACTGCGTGGCAGAAGTCGCGAAAGCCGTTGAGTCCGTGGGTGGACCAATGCTGATCACCGCTGACCCCGGTAACGCTGAACAGATGCGCGATACGGCCACGGGTCAGGCCCACACCGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCCAACTTTCTGACATCGCGCCGGCCAtgtttgcgctg                                                                                                                                                                                                      >  2:202009/1‑243 (MQ=255)                               tAAAGCGGTTGAAGCACTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGt                                                                                                                                                                                                <  1:75628/250‑1 (MQ=255)                                 aaGCGGTTGAATCACTGTATCACTGCGTGGAAGAAGTCGCGAAAGCGGTGGAGTCCGTGGGTGGACATCTGCTGATCACCGCTGACCCCGGTAACGCTGAACATATGCGCGATCCGGCAACGGGTCAGGCACACACGTCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTGTCTTACATCGCGCCGACCATGTTGTCGCTTATGGGTATg                                                                                                                                                                                             <  2:265061/251‑1 (MQ=255)                                        tGAAGCACTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATcc                                                                                                                                                                                      >  2:324395/1‑251 (MQ=255)                                         gAAGCACTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATcc                                                                                                                                                                                      <  1:202009/250‑1 (MQ=255)                                           agcacTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGc                                                                                                                                                                                   >  1:407307/5‑251 (MQ=255)                                             cacTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCa                                                                                                                                                                                  >  1:9408/3‑250 (MQ=255)                                                 ggATCCCTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACCCCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAagcg                                                                                                                                                                              >  2:403331/1‑248 (MQ=255)                                                  gATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGAt                                                                                                                                                                            <  1:146863/251‑1 (MQ=255)                                                          cGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCGGATGGGTATGGAAATCCCGCAAGAGCTGACTGGt                                                                                                                                                                     >  2:298598/1‑250 (MQ=255)                                                           gTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTa                                                                                                                                                                    >  2:170884/1‑250 (MQ=255)                                                              gaagaaGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGcc                                                                                                                                                                <  2:203933/251‑1 (MQ=255)                                                                 gaagTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCt                                                                                                                                                             <  2:431999/251‑1 (MQ=255)                                                                    gTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGtt                                                                                                                                                          >  1:315424/1‑251 (MQ=255)                                                                       gcgAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCGCGGTAACGCTGATCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTTATGGGTATTGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCAt                                                                                                                                                       <  2:127420/251‑1 (MQ=255)                                                                       gcgAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCGGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCGGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCAt                                                                                                                                                       <  2:411484/251‑1 (MQ=255)                                                                                          tgggtggACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTcccc                                                                                                                                    >  1:395725/1‑251 (MQ=255)                                                                                          tgggtagccAACTGCTGATCACCGCTGACAACGGTAACTCTGACCAGATGCGCGCTCCGGCAAAGGGTCAGGCACACCCGGAACACACCACACTCCCAGTTCCGCTGCTTTACGTTGGTGATAAGAACGTGAAAGCGGATGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCACGAGCTGACTGGTAACCCGCTGTTCATCGTGGCATAAtcactccac                                                                                                                                    >  2:183375/9‑249 (MQ=255)                                                                                            ggtggACAACTGCTGATCCCCGCTGACCACGGGAACGCTGAACAGATGTGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCAt                                                                                                                                  <  1:109/251‑1 (MQ=255)                                                                                                gACAACTGCTGATCACCGCTGACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAAGCCCGCAAGAGATGACGGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAgg                                                                                                                              <  2:3415/251‑1 (MQ=255)                                                                                                       gctgATCACCGCTAACCACGGTAACGCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGGCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGaaa                                                                                                                         <  1:270332/249‑1 (MQ=255)                                                                                                                      ccACGGTAACGCTGAACAGATGCGCGATCCGGCAACGTGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGa                                                                                                        <  2:238199/251‑1 (MQ=255)                                                                                                                                gCTGAACAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGt                                                                                              >  1:7296/1‑251 (MQ=255)                                                                                                                                     acAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAAc                                                                                         >  1:416327/2‑251 (MQ=255)                                                                                                                                       aGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACcg                                                                                       >  1:138432/1‑251 (MQ=255)                                                                                                                                         aTGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACcgcg                                                                                     >  1:216212/1‑251 (MQ=255)                                                                                                                                                  ccGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTg                                                                             <  1:118169/250‑1 (MQ=255)                                                                                                                                                                    cacacggcccccacCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATAAGGCCCATCATCt                                                           >  2:374059/1‑250 (MQ=255)                                                                                                                                                                     acacggcacacacCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGAAATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTAc                                                         >  1:292046/1‑251 (MQ=255)                                                                                                                                                                            acccaccacCCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCg                                                  <  2:103653/248‑1 (MQ=255)                                                                                                                                                                               cacCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTc                                               <  2:9408/251‑1 (MQ=255)                                                                                                                                                                                   aaCCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCCGCGTTCtt                                             >  2:332185/1‑249 (MQ=255)                                                                                                                                                                                    aCCTGCCAGTTCCGCTGATTTAAGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAgc                                          <  2:270320/251‑1 (MQ=255)                                                                                                                                                                                     ccTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGCCATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCTGGTTTGCACTAAGGCCCATCATCTACGCCAGCGTTCTTAgcg                                         <  1:127262/251‑1 (MQ=255)                                                                                                                                                                                       tGCCAGTTCCGATGATTTACGTTGGTGCTAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACCTCGAGCCGACCATGTTGTCGCTGATGGGTATGGCAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCt                                       >  2:376519/1‑251 (MQ=255)                                                                                                                                                                                                           gTTGGTGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCtt                   <  1:14192/251‑1 (MQ=255)                                                                                                                                                                                                                tGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCa                                                     >  2:162673/1‑212 (MQ=255)                                                                                                                                                                                                                tGATAAGAACGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCa                                                     <  1:162673/212‑1 (MQ=255)                                                                                                                                                                                                                         cgggAAAGCGGTTGCAGGCGGCAAACTTTCTGACAGCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCgg     <  1:211660/248‑1 (MQ=255)                                                                                                                                                                                                                         cGTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCAGGTTGGCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCatgctg     <  2:1602/251‑6 (MQ=255)                                                                                                                                                                                                                          gTGAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCgg     >  1:217138/1‑250 (MQ=255)                                                                                                                                                                                                                            gAAAGCGGTTGCAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCATGCGGATg  >  1:81551/1‑251 (MQ=255)                                                                                                                                                                                                                            gAAAGCGGGTGCAGGCGGCAACCTTTCTGACATCGCGCCTACCTTGTTTTCGCTGATGGGTATGGAAATCCCGCAAGAGATTACTGGTAAGTCGCGGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATCCCATGACACGGGCCGTGAAACCGCGCAGGTTTGTAATCAGGCCCATCATCTACGCCAGCGTGCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGTCCATGCGGATg  <  2:276338/251‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                AGGTCACACCGGGGTGATGGAAGCGGCGGTTAAAGCGGTTGAAGCGCTGGATCACTGCGTGGAAGAAGTCGCGAAAGCGGTTGAGTCCGTGGGTGGACAACTGCTGATCACCGCTGACCACGGTAACGCTGAGCAGATGCGCGATCCGGCAACGGGTCAGGCACACACGGCACACACCAACCTGCCAGTTCCGCTGATTTACGTTGGTGATAAGAACGTGAAAGCGGTTGAAGGCGGCAAACTTTCTGACATCGCGCCGACCATGTTGTCGCTGATGGGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATG  >  NC_000913/3786474‑3786943

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