Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,163,845	A→G	D38D (GAT→GAC)	unknown ←	tRNA 4‑thiouridine synthase (EC 2.8.1.4) / Rhodanese‑like domain required for thiamine synthesis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,163,845	0	A	G	100.0%	182.2 / NA	67	D38D (GAT→GAC)	unknown	tRNA 4‑thiouridine synthase (EC 2.8.1.4) / Rhodanese‑like domain required for thiamine synthesis
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/38);  total (29/38)

AATTAAATTCATGTTTGCCTTTACGTTTCACCCGCACGCAAAAACTCTTGTTTTCCAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCATCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAACCCTCTTAATATAAAAATATATATAA  >  1/1163598‑1164070                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                  aaTTAAATTCATGTTTGCCTTTACGTTTCACCCGCACGCAAAAACTCTTGTTTTCCAGTTGCGGAGCCACATCAGACAGGGGTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGtcgtcg                                                                                                                                                                                                                                <  1:428310/251‑1 (MQ=255)     aaaTTCATGTTTGCCTTTACGTTTCACCCGCACGCAAAAACTCTTGTTTTCCAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGtcgtcg                                                                                                                                                                                                                                <  1:877602/247‑1 (MQ=255)               ttGCCTTTACGTTTCACCCGCACGCAAAAACTCTTGTTTTCCAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATAc                                                                                                                                                                                                                  <  1:379092/251‑1 (MQ=255)               ttGCCTTTACGTTTCACCCGCACGCAAAAACTCTTGTTTTCCAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATAc                                                                                                                                                                                                                  <  1:381426/251‑1 (MQ=255)                tGCCTTTACGTTTCACCCGCACGCAAAAACTCTTGTTTTCCAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACg                                                                                                                                                                                                                 <  1:484326/251‑1 (MQ=255)                        cGTTTCACCCGCACGCAAAAACTCTTGTTTTCCAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAat                                                                                                                                                                                                         <  1:579090/251‑1 (MQ=255)                        cGTTTCACCCGCACGCAAAAACTCTTGTTTTCCAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAat                                                                                                                                                                                                         <  1:620491/251‑1 (MQ=255)                        cGTTTCACCCGCACGCAAAAACTCTTGTTTTCCAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAat                                                                                                                                                                                                         >  1:841354/1‑251 (MQ=255)                         gTTTCACCCGCACGCAAAAACTCTTGTTTTCCAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAat                                                                                                                                                                                                         <  1:718307/250‑1 (MQ=255)                             cACCCGCACGCAAAAACTCTTGTTTTCCAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTGATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGc                                                                                                                                                                                                    <  1:739063/251‑1 (MQ=255)                                            aCTCTTGTTTTCCAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGaa                                                                                                                                                                                                          >  1:392443/1‑230 (MQ=255)                                                   ttttCCAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTc                                                                                                                                                                                               <  1:217325/234‑1 (MQ=255)                                                   ttttCCAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAAcc                                                                                                                                                                               <  1:560046/250‑1 (MQ=255)                                                       ccAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTc                                                                                                                                                                          <  1:324482/251‑1 (MQ=255)                                                                  gCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGAtt                                                                                                                                                                >  1:185964/1‑250 (MQ=255)                                                                     acaTCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCg                                                                                                                                                                        <  1:796094/239‑1 (MQ=255)                                                                     acaTCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCg                                                                                                                                                              >  1:20559/1‑249 (MQ=255)                                                                         cAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCtt                                                                                                                                                           <  1:885136/248‑1 (MQ=255)                                                                           gACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAAt                                                                                                                                                      >  1:405621/1‑251 (MQ=255)                                                                                ggTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATc                                                                                                                                                     <  1:855084/247‑1 (MQ=255)                                                                                             atatCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCtt                                                                                                                                                           >  1:170558/1‑228 (MQ=255)                                                                                               atCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAAtt                                                                                                                                  >  1:643776/1‑251 (MQ=255)                                                                                                 cATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTg                                                                                                                                >  1:302184/1‑251 (MQ=255)                                                                                                         ctgctAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATaaa                                                                                                                        >  1:886234/1‑251 (MQ=255)                                                                                                           gctAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACg                                                                                                                                                                                                                 >  1:45361/1‑160 (MQ=255)                                                                                                           gctAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATc                                                                                                                                                     >  1:296652/1‑220 (MQ=255)                                                                                                           gctAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAAtt                                                                                                                      >  1:2084/1‑251 (MQ=255)                                                                                                             tAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTc                                                                                                                    <  1:146122/251‑1 (MQ=255)                                                                                                                                 tGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAAc                                                                                                <  1:401252/251‑1 (MQ=255)                                                                                                                                 tGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAAc                                                                                                <  1:222038/251‑1 (MQ=255)                                                                                                                                                 tCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAAccc                                                                                              >  1:440317/1‑237 (MQ=255)                                                                                                                                                  ccccGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACg                                                                                                                             <  1:860044/205‑1 (MQ=255)                                                                                                                                                           tcGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAg                                                                      >  1:676359/2‑251 (MQ=255)                                                                                                                                                           gCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGaa                                                                                                                                                                                                          <  1:222717/119‑1 (MQ=255)                                                                                                                                                           gCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAg                                                                      >  1:231046/1‑251 (MQ=255)                                                                                                                                                            cGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAg                                                                      >  1:816227/1‑250 (MQ=255)                                                                                                                                                            cGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGc                                                                     <  1:149778/251‑1 (MQ=255)                                                                                                                                                            cGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGc                                                                     <  1:353286/251‑1 (MQ=255)                                                                                                                                                                      gCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACg                                                                                                                                                                                                                 <  1:580027/101‑1 (MQ=255)                                                                                                                                                                      gCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCAc                                                                                                                                                                     >  1:777941/1‑145 (MQ=255)                                                                                                                                                                      gCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGaa                                                                                        <  1:531945/222‑1 (MQ=255)                                                                                                                                                                      gCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGcac                                                            >  1:70682/1‑250 (MQ=255)                                                                                                                                                                         cAATGAGCATCGGTAAATTCTCTTTGATTTTTGACCGCACTTCGATGTAATCCCAA‑GTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACaaaaa                                                       <  1:103152/251‑1 (MQ=255)                                                                                                                                                                                  tCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCAtt                                               >  1:72168/1‑251 (MQ=255)                                                                                                                                                                                    ggTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTgg                                             >  1:819567/1‑251 (MQ=255)                                                                                                                                                                                                   aTTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCtt                                                                                                                                                           >  1:360146/1‑126 (MQ=255)                                                                                                                                                                                                   aTTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTa                              >  1:277565/1‑251 (MQ=255)                                                                                                                                                                                                      tttGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAAc                           <  1:821486/251‑1 (MQ=255)                                                                                                                                                                                                              cACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAAc                           <  1:84212/243‑1 (MQ=255)                                                                                                                                                                                                                cTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTc                                                                                                                    >  1:632612/1‑152 (MQ=255)                                                                                                                                                                                                                cTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCTTAGAAAAATgg                                                                                  <  1:412048/186‑1 (MQ=255)                                                                                                                                                                                                                cTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAAt                                                                                    >  1:410257/1‑184 (MQ=255)                                                                                                                                                                                                                cTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAACCCTCTTAata                 <  1:402712/251‑1 (MQ=255)                                                                                                                                                                                                                cTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAACCCTCTTAata                 >  1:20577/1‑251 (MQ=255)                                                                                                                                                                                                                     aTGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCAc                                                                                                                                                                     <  1:355216/98‑1 (MQ=255)                                                                                                                                                                                                                     aTGTAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAAccc                                                                                              >  1:54966/1‑169 (MQ=255)                                                                                                                                                                                                                        tAATCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAACCCTCTTAATATAAAAat          <  1:333416/250‑1 (MQ=255)                                                                                                                                                                                                                         aaTCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCGCGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAACCCTCTTAATATAAAAatat        <  1:55694/251‑1 (MQ=255)                                                                                                                                                                                                                          aTCCCAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACaaaaa                                                       >  1:653869/1‑203 (MQ=255)                                                                                                                                                                                                                             ccAATGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAACCCTCTTAATATAAAAatatatat    <  1:869924/251‑1 (MQ=255)                                                                                                                                                                                                                               aaTGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAACCCTCTTAATATAAAAATATATATaa  <  1:661987/251‑1 (MQ=255)                                                                                                                                                                                                                                aTGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAACCCTCTTAatat                >  1:397246/1‑236 (MQ=255)                                                                                                                                                                                                                                 tGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGc                                            <  1:721403/207‑1 (MQ=255)                                                                                                                                                                                                                                 tGTTTCACCACGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAAc                                      <  1:321783/213‑1 (MQ=255)                                                                                                                                                                                                                                           cGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATATCACAAAAATTACCATTGGCTCAAACTCCTTTTAAACCCTCTTAATATAAAAatatatat    <  1:741743/237‑1 (MQ=255)                                                                                                                                                                                                                                           cGGCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAAc                           <  1:531921/214‑1 (MQ=255)                                                                                                                                                                                                                                             gCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAAc                           >  1:290352/1‑212 (MQ=255)                                                                                                                                                                                                                                             gCGATGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAACCCTCTTAATATAAAAata         <  1:329478/230‑1 (MQ=255)                                                                                                                                                                                                                                                aTGGTGTCGTCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAAc                           <  1:402762/209‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                  AATTAAATTCATGTTTGCCTTTACGTTTCACCCGCACGCAAAAACTCTTGTTTTCCAGTTGCGGAGCCACATCAGACAGGGTTTGTTCAAAAATATCATGGAGGCTGCTAAAGGGTTTTTCTTCCACTTGTAGAAAATGATGAATCCCCGGGATGCGTTGCAACTGCTCAATGAGCATCGGTAAATTCTCTTCGATTTTTGACCGCACTTCGATGTAATCCCAATGTTTCACCACGGCGATGGTGTCATCGTGCTTTGCCAATACGTTACGAATATTGCCGGTCAGAATTTTCACGAACCGTTTCCGCACGGATTCGCTTTTAATCATAATTTCAGGAAATAATTTGACGATAAATTTCATGTTAATTCACTTTTTAACCCATAGAAAAATGGGGCGTATTTTAGCAGATAGCACAAAAATTACCATTGGCTCAAACTCCTTTTAAACCCTCTTAATATAAAAATATATATAA  >  1/1163598‑1164070

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