Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,148,900	G→A	G335D (GGC→GAC)	gpr →	L‑glyceraldehyde 3‑phosphate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,148,900	0	G	A	100.0%	216.5 / NA	59	G335D (GGC→GAC)	gpr	L‑glyceraldehyde 3‑phosphate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/30);  total (29/30)

GACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGTGAATGGAAGAAATTGCCACAATCGACGCTGCT  >  NC_000913/3148658‑3149145                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       gACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTCGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGa                                                                                                                                                                                                                                                <  2:202538/251‑1 (MQ=255)    acCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATc                                                                                                                                                                                                                                             <  2:35134/251‑1 (MQ=255)          aaTGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCCAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGc                                                                                                                                                                                                                                       >  1:516012/1‑251 (MQ=255)             gCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAgg                                                                                                                                                                                                                                    <  1:111079/251‑1 (MQ=255)             gCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAg                                                                                                                                                                                                                                     <  2:506102/250‑1 (MQ=255)                tACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCg                                                                                                                                                                                                                                  <  1:311225/250‑1 (MQ=255)                tACCGAAGCCAAACGCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGGGCAGGCGCTGAATAATCGGACATTTAGCACCAAGGAGCTGGCGCAGAGTGATCAGCATAGCGCCGATGACGGGCTGAGTCTGTGGCAGGc                                                                                                                                                                                                                                   <  2:243217/249‑1 (MQ=255)                 aCCGAAGCCAACCTCAACAGCCGGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGAGTGAGCAGCAGATCGCCGATGACGTGCTGAATCTGTGGCAGGCGTc                                                                                                                                                                                                                                <  2:330849/251‑1 (MQ=255)                      aGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACCAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCg                                                                                                                                                                                                                           >  1:274689/1‑251 (MQ=255)                                 aCAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTAGGTCTTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCCAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGATGAATCTGTGGAAGGCGTCTTCCGATAAATGACCt                                                                                                                                                                                                                >  1:544445/1‑251 (MQ=255)                                          gcTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGATCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAAc                                                                                                                                                                                                       <  1:280365/251‑1 (MQ=255)                                             tATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTTACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACg                                                                                                                                                                                                      <  1:354177/249‑1 (MQ=255)                                                tgaatgaaATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACAGTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGgcg                                                                                                                                                                                                  <  2:353074/250‑1 (MQ=255)                                                 gaatgaaATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGGTTGGAGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGAGACTTCCGATAAATGACCTGTTAATAACGGgcgcg                                                                                                                                                                                                >  1:90977/1‑251 (MQ=255)                                                       aaTGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAAcc                                                                                                                                                                                          >  2:29862/1‑251 (MQ=255)                                                        aTGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACcg                                                                                                                                                                                         <  2:13624/251‑1 (MQ=255)                                                                  agcGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAatc                                                                                                                                                                               >  2:258971/1‑251 (MQ=255)                                                                                   gcgcAAATGGCGTTAAGCTGGTTGCTGAAAGATGAGCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGGGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCTATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGGGATTATGAc                                                                                                                                                              <  2:488025/251‑1 (MQ=255)                                                                                    cgcAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTGAGCACCAAGGAGCTGGCGCAGATTGAGCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACg                                                                                                                                                             <  1:481692/251‑1 (MQ=255)                                                                                          tGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCg                                                                                                                                                                                   >  1:73896/1‑223 (MQ=255)                                                                                          tGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCg                                                                                                                                                                                   <  2:73896/223‑1 (MQ=255)                                                                                           ggCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGagtcag                                                                                                                                                       <  1:236579/250‑1 (MQ=255)                                                                                                          gCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGcc                                                                                                                                        <  1:473448/250‑1 (MQ=255)                                                                                                            tGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCAt                                                                                                                                     >  1:524657/1‑251 (MQ=255)                                                                                                               aaGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATcac                                                                                                                                  >  1:293882/1‑251 (MQ=255)                                                                                                                    gatCGCGTGACGTCGGTCTTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATg                                                                                                                             >  1:98766/1‑251 (MQ=255)                                                                                                                    gatCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATg                                                                                                                             >  1:120634/1‑251 (MQ=255)                                                                                                                       cgcgTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAg                                                                                                                          <  1:250273/251‑1 (MQ=255)                                                                                                                          gTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAACAGCAGaa                                                                                                                        >  2:483158/1‑250 (MQ=255)                                                                                                                            gACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATGTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGc                                                                                                                     <  2:457423/251‑1 (MQ=255)                                                                                                                                       tgattgGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGGCAAGATAGCCCCTCACAAATGCAGAAAGCACAAACGTc                                                                                                            >  2:106300/1‑249 (MQ=255)                                                                                                                                       tgattgGTACCAGCCGAGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAATGATTATGACGAGTCAGTAGGTCAAGATAGCCCATCACATATTCAGAAAGAACAAAAGTCa                                                                                                           >  1:455968/1‑250 (MQ=255)                                                                                                                                           tgGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATg                                                                                                       <  1:323712/250‑1 (MQ=255)                                                                                                                                                   gCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAc                                                                                              <  1:483158/251‑1 (MQ=255)                                                                                                                                                          ggAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCGGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAAGAACGTGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACcac                                                                                      <  1:355738/251‑1 (MQ=255)                                                                                                                                                              cAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCa                                                                                  <  1:205199/251‑1 (MQ=255)                                                                                                                                                                   agagGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTTTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACTTCAGa                                                                                                         <  2:516390/224‑1 (MQ=255)                                                                                                                                                                   agagGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGa                                                                                                         >  1:516390/1‑224 (MQ=255)                                                                                                                                                                     aggagAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCATTTTGtt                                                                           <  2:257217/251‑1 (MQ=255)                                                                                                                                                                      ggagAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCGGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTGtt                                                                           >  1:124011/1‑250 (MQ=255)                                                                                                                                                                            cGTGCATGCGCTGAATAATCTGACATGTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGGGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTGTTATCCggg                                                                    <  1:394593/251‑1 (MQ=255)                                                                                                                                                                                        gAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTGTTATCCGGGACATTTTTCGc                                                         >  1:350477/1‑250 (MQ=255)                                                                                                                                                                                            aatCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATATGTGGCAGGCGTCTTCCGATAAATGAACTGTTAATAACGGACGCGCGAACCGCGCCCGAATCAGCTCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGCTAATGAC‑GTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATcc                                                     >  1:237383/1‑250 (MQ=255)                                                                                                                                                                                             atCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGCTCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGCCTTCCGATAAATGAACTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGCGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCAAAAACGCCAGATGGCTAAT‑ACACAACCACATCAGTATGTTATCCGGGCCATTTTTCGCATcaat                                                   >  2:45452/1‑248 (MQ=255)                                                                                                                                                                                               cTGACGTTTAGCACCAAGGAGCTGGCGCAGATTGCTCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTACGATAAATGACCTGTTAATAACAGGCGCGCGAAACGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGGCAAGATAACCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATa                                                  >  2:74881/1‑250 (MQ=255)                                                                                                                                                                                                   cATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGAc                                             >  1:471192/1‑251 (MQ=255)                                                                                                                                                                                                   cATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGGTTGTTATCCGGGACATTTTGCGCATCCATAAAGAc                                             <  2:280961/251‑1 (MQ=255)                                                                                                                                                                                                        aGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCa                                        <  2:392750/251‑1 (MQ=255)                                                                                                                                                                                                             cAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGCAAGCACAAACGTCAGATgg                                                                                                      >  1:381846/1‑185 (MQ=255)                                                                                                                                                                                                             cAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATgg                                                                                                      <  2:381846/185‑1 (MQ=255)                                                                                                                                                                                                                        cgcAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GAACCACAACAGTTTGTTATCCGGGACATTTTTCGCATCCATCAAGACGCGCAGTAAGTGAATGga                           >  1:420609/1‑248 (MQ=255)                                                                                                                                                                                                                        cgcAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTAAGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGTGAATGga                           >  1:67038/1‑248 (MQ=255)                                                                                                                                                                                                                           aGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTTTTAGCCGGGACATTTTTCGCATCCATAAAGGCGCGCAGTAAGTGAATGGAAGAAAt                     <  2:516012/251‑1 (MQ=255)                                                                                                                                                                                                                              ttGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGTGAATGGAAGAAATTGc                  >  1:211840/1‑251 (MQ=255)                                                                                                                                                                                                                                 aTCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTAGGACAAGATAGCCCATCACAAATGCAGAAAAACCAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTGTTATCCGGGACATTTTTCTCATCCATAAAGACGCGCAGTAAGTGAATGGAAGAAATTGCcac               >  1:528132/1‑251 (MQ=255)                                                                                                                                                                                                                                         atCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGTGCATGGCAGAAATTGCCACAATCGAcc       >  2:125251/1‑250 (MQ=255)                                                                                                                                                                                                                                            gCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGTGAATGGAAGAAAGTGCCACAATCGACgctg    <  2:319588/251‑1 (MQ=255)                                                                                                                                                                                                                                              cGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGTGAATGGAAGAAATTGCCACAATCGACgctgct  >  2:403769/1‑251 (MQ=255)                                                                                                                                                                                                                                               gATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGTCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGACCATCACAAATACAGAAAGCACAAACGTCAGATGGCTAATGAC‑GTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGTGAATGGCAGAAATTGCCACAATCGACgcttct  >  2:226187/1‑247 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       GACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGAC‑GTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGTGAATGGAAGAAATTGCCACAATCGACGCTGCT  >  NC_000913/3148658‑3149145

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