Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	242,283	A→G	S39G (AGT→GGT)	unknown →	Ferrous‑iron efflux pump FieF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	242,283	0	A	G	100.0%	217.9 / NA	78	S39G (AGT→GGT)	unknown	Ferrous‑iron efflux pump FieF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/41);  total (37/41)

TGAAAATTCCGCGCACCCATAAATTACGCAATATGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCAGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCCGCGATTTTCCTGTTATTGCAGGGCATTCATCGTTTTAATTCACCGCAGAGGCGCTGACTAACACGTCGCTCGGCAT  >  1/242036‑242526                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    tGAAAATTCCGCGCACCCATAAATTACGCAATATGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTa                                                                                                                                                                                                                                                  >  1:723025/1‑251 (MQ=255)   aaaaTTCCGCGCACCCATAAATTACGCAATATGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTAt                                                                                                                                                                                                                                                 <  1:119130/250‑1 (MQ=255)   aaaaTTCCGCGCACCCATAAATTACGCAATATGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATg                                                                                                                                                                                                                                                <  1:890774/251‑1 (MQ=255)                cccATAAATTACGCAATATGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATaa                                                                                                                                                                                                                                   >  1:773067/1‑251 (MQ=255)                cccATAAATTACGCAATATGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATaa                                                                                                                                                                                                                                   >  1:479655/1‑251 (MQ=255)                cccATAAATTACGCAATATGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATaa                                                                                                                                                                                                                                   >  1:363173/1‑251 (MQ=255)                   aTAAATTACGCAATATGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAATGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACgg                                                                                                                                                                                                                                <  1:454758/251‑1 (MQ=255)                    gaaaTTACGCAATATGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACgg                                                                                                                                                                                                                                <  1:636563/249‑1 (MQ=255)                     aaaTTACGCAATATGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGAt                                                                                                                                                                                                                              >  1:356765/1‑251 (MQ=255)                     aaaTTACGCAATATGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGAt                                                                                                                                                                                                                              <  1:355427/251‑1 (MQ=255)                          aCGCAATATGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCAt                                                                                                                                                                                                                         >  1:725786/1‑251 (MQ=255)                             cAATATGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATg                                                                                                                                                                                                                        >  1:712081/1‑249 (MQ=255)                               atatGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTgg                                                                                                                                                                                                                    >  1:686917/1‑251 (MQ=255)                                 atGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGAc                                                                                                                                                                                                                  <  1:765464/251‑1 (MQ=255)                                        aTAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTgg                                                                                                                                                                                                           >  1:232849/1‑251 (MQ=255)                                             aaaaTCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTcc                                                                                                                                                                                                      <  1:582725/251‑1 (MQ=255)                                                 tcatcaGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTa                                                                                                                                                                                                  <  1:622768/251‑1 (MQ=255)                                                   atcaGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATg                                                                                                                                                                                                <  1:703961/251‑1 (MQ=255)                                                        ccACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAg                                                                                                                                                                                              <  1:759119/248‑1 (MQ=255)                                                        ccACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAg                                                                                                                                                                                              <  1:899416/248‑1 (MQ=255)                                                        ccACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGc                                                                                                                                                                                             <  1:202890/249‑1 (MQ=255)                                                                  aaaTTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGAttt                                                                                                                                                                                   <  1:368936/249‑1 (MQ=255)                                                                               aTAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATg                                                                                                                                                                                                >  1:780134/1‑223 (MQ=255)                                                                                 aaaGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATg                                                                                                                                                                  <  1:538110/251‑1 (MQ=255)                                                                                 aaaGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTACTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATg                                                                                                                                                                  <  1:791370/251‑1 (MQ=255)                                                                                  aaGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGTTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGc                                                                                                                                                                 <  1:873503/251‑1 (MQ=255)                                                                                         gAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATg                                                                                                                                                                                                >  1:546070/1‑213 (MQ=255)                                                                                                cTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATaa                                                                                                                                                    <  1:835252/250‑1 (MQ=255)                                                                                                cTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAAt                                                                                                                                                   <  1:647225/251‑1 (MQ=255)                                                                                                  aaCGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCa                                                                                                                                                 <  1:509587/251‑1 (MQ=255)                                                                                                          cGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACgg                                                                                                                                                                                                                                >  1:751664/1‑164 (MQ=255)                                                                                                                   cGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGc                                                                                                                                 <  1:151055/250‑1 (MQ=255)                                                                                                                     ttGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAg                                                                                                                              >  1:130709/1‑251 (MQ=255)                                                                                                                      tGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCC‑TTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGc                                                                                                                            <  1:683739/251‑1 (MQ=255)                                                                                                                       gggATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAg                                                                                                                                                                                              <  1:284905/185‑1 (MQ=255)                                                                                                                       gggATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAgg                                                                                                                             >  1:858286/1‑250 (MQ=255)                                                                                                                            gTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGc                                                                                                                                                                                                                       >  1:856589/1‑155 (MQ=255)                                                                                                                            gTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGc                                                                                                                                                                                             <  1:375095/181‑1 (MQ=255)                                                                                                                            gTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAAt                                                                                                                       >  1:191994/1‑251 (MQ=255)                                                                                                                              agaagaTAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCg                                                                                                                     >  1:16988/1‑251 (MQ=255)                                                                                                                                  gaTAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTgg                                                                                                                 >  1:622496/1‑251 (MQ=255)                                                                                                                                  gaTAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTgg                                                                                                                 >  1:331182/1‑251 (MQ=255)                                                                                                                                              tGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGc                                                                                                     <  1:101618/251‑1 (MQ=255)                                                                                                                                                tCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCaa                                                                                                   <  1:382299/251‑1 (MQ=255)                                                                                                                                                 tgCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGAttt                                                                                                                                                                                   <  1:430316/169‑1 (MQ=255)                                                                                                                                                 cGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCaaa                                                                                                  >  1:651539/1‑251 (MQ=255)                                                                                                                                                 cGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGGGCaaa                                                                                                  >  1:846731/1‑251 (MQ=255)                                                                                                                                                          cgcgCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTa                                                                                         <  1:330813/251‑1 (MQ=255)                                                                                                                                                          cgcgCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTa                                                                                         <  1:313568/251‑1 (MQ=255)                                                                                                                                                              cAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGC‑AAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCaa                                                                                                   <  1:859942/236‑1 (MQ=255)                                                                                                                                                                ggTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAg                                                                                                 <  1:867945/237‑1 (MQ=255)                                                                                                                                                                 gTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTc                                                                                     <  1:695202/248‑1 (MQ=255)                                                                                                                                                                                cATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTgcgc                                                                                                                                                                              >  1:561496/1‑144 (MQ=255)                                                                                                                                                                                    ttGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCAt                                                                                                                                                                                                                         >  1:103238/1‑97 (MQ=255)                                                                                                                                                                                      gCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTAtt                                                             >  1:260997/1‑251 (MQ=255)                                                                                                                                                                                        cATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCa                                                                                                                                                 >  1:63824/1‑165 (MQ=255)                                                                                                                                                                                         attattACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCa                                                          >  1:792384/1‑251 (MQ=255)                                                                                                                                                                                          ttattaCTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAAt                                                                                                                                                   <  1:366681/161‑1 (MQ=255)                                                                                                                                                                                          ttattaCTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGc                                                           >  1:297607/1‑249 (MQ=255)                                                                                                                                                                                             ttaCTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCAggg                                                       >  1:600726/1‑250 (MQ=255)                                                                                                                                                                                                     cTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTa                                                                                         >  1:883017/1‑208 (MQ=255)                                                                                                                                                                                                       gATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAg                                                                                                                                                                                              >  1:552535/1‑105 (MQ=255)                                                                                                                                                                                                       gATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCAGGGc                                                      >  1:247013/1‑241 (MQ=255)                                                                                                                                                                                                              aTCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCAGGGCATTCATCGTTTTAATTc                                     <  1:604424/251‑1 (MQ=255)                                                                                                                                                                                                              aTCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCAGGGCATTCATCGTTTTAATTc                                     <  1:556555/251‑1 (MQ=255)                                                                                                                                                                                                              aTCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCAGGGCATTCATCGTTTTAATTc                                     <  1:334995/251‑1 (MQ=255)                                                                                                                                                                                                                           gCATTGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCAGGGCATTCATCGTTTTAATTc                                     >  1:684166/1‑238 (MQ=255)                                                                                                                                                                                                                              ttGCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCAGGGCATTCATCGTTTTAATTCACCGCAGAGGCGCTGa                     >  1:185529/1‑251 (MQ=255)                                                                                                                                                                                                                                gCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCAGGGCATTCATCGTTTTAATTCACCGCAGAGGCGCTGa                     >  1:33250/1‑249 (MQ=255)                                                                                                                                                                                                                                gCCTGGTGGCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCAGGGCATTCATCGTTTTAATTCACCGCAGAGGCGCTGAc                    <  1:146379/250‑1 (MQ=255)                                                                                                                                                                                                                                     gtggCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGc                                                           <  1:768497/206‑1 (MQ=255)                                                                                                                                                                                                                                     gtggCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCAGGGCATTCATCGTTTTAATTCACCGCAGAGGCGCTGACTAaca               >  1:894020/1‑250 (MQ=255)                                                                                                                                                                                                                                      tggCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCAGGGCATTCATCGTTTTAATTCACCGCAGAGGCGCTGACTAacac              <  1:623623/250‑1 (MQ=255)                                                                                                                                                                                                                                      tggCAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCAGGGCATTCATCGTTTTAATTCACCGCAGAGGCGCTGACTAacac              <  1:719573/250‑1 (MQ=255)                                                                                                                                                                                                                                         cAAACGGGCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCAGGGCATTCATCGTTTTAATTCACCGCAGAGGCGCTGACTAACACGTCg          >  1:789493/1‑251 (MQ=255)                                                                                                                                                                                                                                               ggCTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAg                                                                                                                                                                                              >  1:701312/1‑65 (MQ=255)                                                                                                                                                                                                                                                 cTCGGTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCGGCGATTTTCCTGTTATTGCAGGGCATTCATCGTTTTAATTCACCGCAGAGGCGCTGACTAACACGTCGCTCGGCAt  <  1:782821/251‑1 (MQ=255)                                                                                                                                                                                                                                                    ggTCGGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTc                                                                                                                                      <  1:558108/116‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    TGAAAATTCCGCGCACCCATAAATTACGCAATATGAGCAATAAAAAAATCATCAGCCACTTGAAGAAATTCTTCGTAGATAAAGAAATGAAAAAGCTAACGCAGGCGGATTGGGCGTTGGGATGTAGAAGATAATGCCCGCTGTCGCTTATTCCGCGCAGGTGAAAAAAGCCGCGCATTTTGCCATTATTACTGCGCTGATTTTAATCTTAGCAAAAGGCATTGCCTGGTGGCAAACGGGCTCGGTCAGTATGCTTGCCTCCATAACGGATTCCATGCTGGACTTGCTGGCGTCCTTTATGAGCATGTTGATTTTGCGCTTTGCGTTAATGCCGGCGGATCATAATCATTCTTTCGGTCACGGCAAGGCGGAATCGTTGGCGTCGCTGGTGCAAAGCGCGTTTATTTCCGGTTCCGCGATTTTCCTGTTATTGCAGGGCATTCATCGTTTTAATTCACCGCAGAGGCGCTGACTAACACGTCGCTCGGCAT  >  1/242036‑242526

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