Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,697,627	T→C	R249R (CGT→CGC)	unknown →	Sulfate transporter, CysZ‑type

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,697,627	0	T	C	100.0%	228.4 / NA	69	R249R (CGT→CGC)	unknown	Sulfate transporter, CysZ‑type
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/46);  total (23/46)

ACCTTTTTATTCACCGCATGGATGATGGCGATTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCAGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGTTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACATATTATGACTATTTACAGTGACAATTCTTATTCTATCGGCAATACCCCGTTAGTACGTTTAAAACAT  >  1/1697385‑1697868                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  aCCTTTTTATTCACCGCATGGATGATGGCGATTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGcaac                                                                                                                                                                                                                                            <  1:415652/250‑1 (MQ=255)    tttttATTCACCGCATGGATGATGGCGATTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAgg                                                                                                                                                                                                                                        >  1:818051/1‑251 (MQ=255)    tttttATTCACCGCATGGATGATGGCGATTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAgg                                                                                                                                                                                                                                        >  1:799998/1‑251 (MQ=255)    tttttATTCACCGCATGGATGATGGCGATTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAgg                                                                                                                                                                                                                                        >  1:690127/1‑251 (MQ=255)     ttttATTCACCGCATGGATGATGGCGATTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAggg                                                                                                                                                                                                                                       >  1:478690/1‑251 (MQ=255)     ttttATTCACCGCATGGATGATGGCGATTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAggg                                                                                                                                                                                                                                       <  1:850699/251‑1 (MQ=255)           tCACCGCATGGATGATGGCGATTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAAcc                                                                                                                                                                                                                                 <  1:852097/251‑1 (MQ=255)              ccGCATGGATGATGGCGATTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCgg                                                                                                                                                                                                                               >  1:617797/1‑250 (MQ=255)               cGCATGGATGATGGCGATTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCa                                                                                                                                                                                                                             >  1:160460/1‑251 (MQ=255)                gCATGGATGATGGCGATTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCa                                                                                                                                                                                                                             >  1:475446/1‑250 (MQ=255)                       tgatgGCGATTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAgtg                                                                                                                                                                                                                     <  1:878690/251‑1 (MQ=255)                          tgGCGATTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAgtgaat                                                                                                                                                                                                                  <  1:842386/251‑1 (MQ=255)                               aTTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTGGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAgtgaatgtgaa                                                                                                                                                                                                             <  1:260309/251‑1 (MQ=255)                               aTTCAATATTGCGACTATCCTTGTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAgtgaatgtgaa                                                                                                                                                                                                             <  1:741131/251‑1 (MQ=255)                                   aaTATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAAtcgt                                                                                                                                                                                                         <  1:546911/251‑1 (MQ=255)                                    atatTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAAtcgt                                                                                                                                                                                                         >  1:678202/1‑250 (MQ=255)                                               aTCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTGACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCgg                                                                                                                                                                                              >  1:428200/1‑250 (MQ=255)                                                 ccTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAggg                                                                                                                                                                                                                                       <  1:571174/207‑1 (MQ=255)                                                  cTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGAc                                                                                                                                                                                          <  1:144773/251‑1 (MQ=255)                                                           aCCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTAt                                                                                                                                                                                  >  1:284950/1‑250 (MQ=255)                                                            ccACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCg                                                                                                                                                                                <  1:693036/251‑1 (MQ=255)                                                              acaAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTc                                                                                                                                                                              <  1:158419/251‑1 (MQ=255)                                                                     cccTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAAt                                                                                                                                                                       <  1:572762/251‑1 (MQ=255)                                                                         ttCGGCATAATGAAAAACGAACAGGCGATTAAGCGCAATGGGACGATGACCTTCGGCAGTTTTGTGACGCTTTGTTCTTTTGTGACGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCt                                                                                                                                                                   <  1:647987/251‑1 (MQ=255)                                                                             gCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACAATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCg                                                                                                                                                                                <  1:55767/234‑1 (MQ=255)                                                                                  aTGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTtat                                                                                                                                                          >  1:544058/1‑251 (MQ=255)                                                                                     aaaaaCGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCACTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTatat                                                                                                                                                        <  1:292740/250‑1 (MQ=255)                                                                                      aaaaCGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGa                                                                                                                                                      <  1:466846/251‑1 (MQ=255)                                                                                      aaaaCGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGa                                                                                                                                                      <  1:526332/251‑1 (MQ=255)                                                                                         aCGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGAttt                                                                                                                                                   >  1:323783/1‑251 (MQ=255)                                                                                                  cGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAAt                                                                                                                                          <  1:494788/251‑1 (MQ=255)                                                                                                         gcgcAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTAtt                                                                                                                                   >  1:342101/1‑251 (MQ=255)                                                                                                            cAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAgtgaatgtgaat                                                                                                                                                                                                            <  1:99153/175‑1 (MQ=255)                                                                                                              aTGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTtata                                                                                                                                                         <  1:291358/224‑1 (MQ=255)                                                                                                              aTGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGtttt                                                                                                                              >  1:273070/1‑251 (MQ=255)                                                                                                                 tgACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGa                                                                                                                                                      <  1:227040/224‑1 (MQ=255)                                                                                                                  gACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGt                                                                                                                                                             <  1:100995/216‑1 (MQ=255)                                                                                                                                      tGGTGACGCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTAt                                                                                                                            >  1:332868/1‑229 (MQ=255)                                                                                                                                             gCTTTGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTTTCAttcgctttctctttta                                                                                                >  1:695100/1‑235 (MQ=255)                                                                                                                                                ttGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGtttttattcttaaaataaaacccc                                                                                              >  1:782470/1‑227 (MQ=255)                                                                                                                                                 tGTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTc                                                                                           <  1:757248/251‑1 (MQ=255)                                                                                                                                                  gTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAgtgaat                                                                                                                                                                                                                  <  1:377389/131‑1 (MQ=255)                                                                                                                                                  gTACTTTTGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATcgg                                                                                                                         >  1:679828/1‑220 (MQ=255)                                                                                                                                                      ttttGTGCCGGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCattat                                                                                      <  1:125886/251‑1 (MQ=255)                                                                                                                                                                  ggtgAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTc                                                                                           <  1:486146/234‑1 (MQ=255)                                                                                                                                                                    tgAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGgag                                                                         <  1:726741/250‑1 (MQ=255)                                                                                                                                                                           ggTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGa                                                                                                                                                      <  1:828111/166‑1 (MQ=255)                                                                                                                                                                               attattCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACATATTATg                                                              <  1:730460/250‑1 (MQ=255)                                                                                                                                                                                ttattCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACATATTATGAc                                                            <  1:797682/251‑1 (MQ=255)                                                                                                                                                                                    tCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGCAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCg                                                                                                                                                                                >  1:790355/1‑131 (MQ=255)                                                                                                                                                                                                gtgtgGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGcc                                                                                                <  1:672938/199‑1 (MQ=255)                                                                                                                                                                                                gtgtgGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGGTTTCATTGTAAAAAAAAAACCACTTCATTATTTTTTCTCTCGAGAAAAAATTATGAATATTTatag                                                   >  1:674852/1‑241 (MQ=255)                                                                                                                                                                                                   tgGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACATATTATGACTATTTACAGTGACAATTc                                          <  1:202669/250‑1 (MQ=255)                                                                                                                                                                                                        cAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACATATTATGACTATTTACAGTGACAATTCTTATTc                                    <  1:463747/251‑1 (MQ=255)                                                                                                                                                                                                        cAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACATATTATGACTATTTACAGTGACAATTCTTATTc                                    <  1:416205/251‑1 (MQ=255)                                                                                                                                                                                                          aCGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGgag                                                                         <  1:460053/212‑1 (MQ=255)                                                                                                                                                                                                          aCGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACATATTATGACTATTTACAGTGACAATTCTTATTCTa                                  <  1:808072/251‑1 (MQ=255)                                                                                                                                                                                                             gCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACATATTATGACTATTTACAGTGACAATTCTTAtt                                     <  1:153383/245‑1 (MQ=255)                                                                                                                                                                                                              cAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACATATTATGACTATTTACAGTGACAATTCTTATTCTATCg                               <  1:773830/250‑1 (MQ=255)                                                                                                                                                                                                                   tgGGTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTCAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACATATTATGACTATTTACAGTGACAATTCTTATTCTATCGGCAAt                          <  1:113009/250‑1 (MQ=255)                                                                                                                                                                                                                    gggTAACGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATg                                                                                                                                                       >  1:509274/1‑124 (MQ=255)                                                                                                                                                                                                                     ggTAACGGAATACCGTGGTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACATATTATGACTATTTACAGTGACAATTCTTATTCTATCGGCAATAc                        <  1:747434/250‑1 (MQ=255)                                                                                                                                                                                                                         aCGGAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCattat                                                                                      <  1:334369/184‑1 (MQ=255)                                                                                                                                                                                                                            gAATACCGTGCTTTATTTAAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACATATTATGACTATTTACAg                                                   <  1:901312/216‑1 (MQ=255)                                                                                                                                                                                                                               tACCGTGCTTTATTTAAGCGCTCGCAACAAGCGCAAACCGGCACTGAAGTGAATGTGCATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACATATTATGACTATTTACAGTGACAATTCTTATTCTATCGGCAATACCCCGTTAGTAc             <  1:252353/251‑1 (MQ=255)                                                                                                                                                                                                                                        ttatttaAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATTCTCAGGAGAACATATTATGACTATTTACAGTGACAATTCTTATTCTATCGGCAATACCCCGt                   <  1:348332/236‑1 (MQ=255)                                                                                                                                                                                                                                          atttaAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTc                                                                                                                                                                              >  1:861286/1‑79 (MQ=255)                                                                                                                                                                                                                                          atttaAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTTTTGCTTAAAATAAAGCCACCTCATTATTTTTTTTCTGGAGAACATAATATTAATATTTACAGTGACAAATCTTTTTTTTTTTGCAACACCCCCCTAGTACGTTTaaaaaaa  >  1:838676/1‑248 (MQ=255)                                                                                                                                                                                                                                          atttaAGCGCTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACAtatt                                                                 <  1:656674/188‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  ACCTTTTTATTCACCGCATGGATGATGGCGATTCAATATTGCGACTATCCTTTTGATAACCACAAAATCCCTTTCGGCATAATGAAAAACGAACTGGCGATTAAGCGCAATGTGACGATGACCTTCGGCAGTTTGGTGACGCTTTGTACTTTTGTGCCAGTGGTGAATCTGGTGATTATTCCCGTGGCAGTGTGTGGCGCAACGGCAATGTGGGTAACGGAATACCGTGCTTTATTTAAGCGTTCGCAACAAGGGCAAACCGGCACTGAAGTGAATGTGAATCGTCGCTCTTCCGGCGACGGTGTTATCGTCAGAAAATAGCTTTTGGTTATATGATTTATAAAAAATAGCTATTGTTTTATCGGCGGTTCTCATGCTAAAACTAACGCCACTTCATTATTTATCCTCAGGAGAACATATTATGACTATTTACAGTGACAATTCTTATTCTATCGGCAATACCCCGTTAGTACGTTTAAAACAT  >  1/1697385‑1697868

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