Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,688,038	A→G	A108A (GCA→GCG)	unknown →	Fe2+ ABC transporter, permease protein 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,688,038	0	A	G	100.0%	223.0 / NA	77	A108A (GCA→GCG)	unknown	Fe2+ ABC transporter, permease protein 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (43/34);  total (43/34)

ATTTCCGCCTTGCTTGCGGTGTCCATCGACATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCACAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTGCCAAACAACATCAGTGGAAAATAGG  >  1/1687805‑1688274                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             gtcttgcgccttgccgtttggtCCATCGACATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGc                                                                                                                                                                                                                                   >  1:162367/21‑245 (MQ=255) gtcttgcgccttgccgtttggtCCATCGACATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGc                                                                                                                                                                                                                                   <  1:280194/225‑1 (MQ=255)   ttCCGCCTTGCTTGCGGTGTCCATCGACATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGAc                                                                                                                                                                                                                           <  1:473823/251‑1 (MQ=255)           tgcttgcGGTGTCCATCGACATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGc                                                                                                                                                                                                                                   <  1:173626/235‑1 (MQ=255)            gcttgcGGTGTCCATCGACATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCa                                                                                                                                                                                                                  >  1:672520/1‑251 (MQ=255)                     gtCCATCGACATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAg                                                                                                                                                                                                         >  1:111632/1‑251 (MQ=255)                     gtCCATCGACATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAg                                                                                                                                                                                                         >  1:635770/1‑251 (MQ=255)                     gtCCATCGACATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACA‑TTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGa                                                                                                                                                                                                        >  1:805291/1‑251 (MQ=255)                      tCCATCGACATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCTATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAg                                                                                                                                                                                                         <  1:435479/250‑1 (MQ=255)                       ccATCGACATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGa                                                                                                                                                                                                        >  1:746946/1‑250 (MQ=255)                       ccATCGACATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGAc                                                                                                                                                                                                       <  1:247511/251‑1 (MQ=255)                       ccATCGACATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGAc                                                                                                                                                                                                       <  1:653168/251‑1 (MQ=255)                             aCATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGAc                                                                                                                                                                                                 <  1:356868/251‑1 (MQ=255)                                 cGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATtc                                                                                                                                                                                             <  1:876739/251‑1 (MQ=255)                                        aaaaTGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCa                                                                                                                                                                                       >  1:335933/1‑250 (MQ=255)                                              gCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGttttt                                                                                                                                                                                >  1:880776/1‑251 (MQ=255)                                                 aaaGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGa                                                                                                                                                                             <  1:807431/251‑1 (MQ=255)                                                       cTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGAtcatca                                                                                                                                                                       >  1:452809/1‑251 (MQ=255)                                                            aTCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGAc                                                                                                                                                                                                       >  1:614127/1‑214 (MQ=255)                                                            aTCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCg                                                                                                                                                                  >  1:678913/1‑251 (MQ=255)                                                             tCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCg                                                                                                                                                                  >  1:885292/1‑250 (MQ=255)                                                             tCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGc                                                                                                                                                                 <  1:777821/251‑1 (MQ=255)                                                                       ccAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGAt                                                                                                                                                       >  1:286365/1‑251 (MQ=255)                                                                       ccAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGAt                                                                                                                                                       <  1:263040/251‑1 (MQ=255)                                                                          atatTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGaa                                                                                                                                                    >  1:37033/1‑251 (MQ=255)                                                                            atTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGa                                                                                                                                                  <  1:534515/251‑1 (MQ=255)                                                                                     gAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTAcc                                                                                                                                          <  1:900905/250‑1 (MQ=255)                                                                                        cccGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCacac                                                                                                                                      <  1:728098/251‑1 (MQ=255)                                                                                        cccGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCacac                                                                                                                                      <  1:833907/251‑1 (MQ=255)                                                                                        cccGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCacac                                                                                                                                      <  1:342298/251‑1 (MQ=255)                                                                                          cGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGAc                                                                                                                                                                                                       >  1:226833/1‑184 (MQ=255)                                                                                                  gcgcCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCaaaaa                                                                                                                            >  1:499982/1‑251 (MQ=255)                                                                                                     ccATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCgg                                                                                                                                                                                         <  1:145886/187‑1 (MQ=255)                                                                                                     ccATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTAcc                                                                                                                                          >  1:695372/1‑234 (MQ=255)                                                                                                     ccATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCacac                                                                                                                                      <  1:203480/238‑1 (MQ=255)                                                                                                                gatgaAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTAc                                                                                                                                           >  1:34148/1‑222 (MQ=255)                                                                                                                   gaAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGc                                                                                                                                                                                                                                   >  1:568264/1‑131 (MQ=255)                                                                                                                   gaAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATtc                                                                                                                                                                                             >  1:535553/1‑169 (MQ=255)                                                                                                                        gAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAgg                                                                                                      >  1:209145/1‑251 (MQ=255)                                                                                                                           ccGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTgg                                                                                                   <  1:586462/251‑1 (MQ=255)                                                                                                                             gTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCacac                                                                                                                                      >  1:631725/1‑214 (MQ=255)                                                                                                                             gTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCacac                                                                                                                                      >  1:486399/1‑214 (MQ=255)                                                                                                                              taataaTTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCacac                                                                                                                                      >  1:450860/1‑213 (MQ=255)                                                                                                                                ataatTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCacac                                                                                                                                      <  1:363803/211‑1 (MQ=255)                                                                                                                                       tctcGCCACGCAAGATTGTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTGCGCGCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGACAAGGCGCAAGACATCGGCATTCTCGGCGCGTTTTTTGATCATCAAATCGCCGTGCGGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATggg                                                                                       <  1:688282/251‑1 (MQ=255)                                                                                                                                              aCGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGa                                                                                <  1:865044/251‑1 (MQ=255)                                                                                                                                              aCGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGa                                                                                <  1:764558/251‑1 (MQ=255)                                                                                                                                              aCGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGCTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGa                                                                                <  1:453718/251‑1 (MQ=255)                                                                                                                                               cGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGAt                                                                               <  1:319/251‑1 (MQ=255)                                                                                                                                                    gATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGAc                                                                                                                                                                                                       <  1:344113/126‑1 (MQ=255)                                                                                                                                                                      gAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACa                                                        >  1:805610/1‑251 (MQ=255)                                                                                                                                                                           gCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGtt                                                    >  1:253693/1‑250 (MQ=255)                                                                                                                                                                           gCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGtt                                                    >  1:383277/1‑250 (MQ=255)                                                                                                                                                                            ccGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTcc                                                  <  1:11916/251‑1 (MQ=255)                                                                                                                                                                                  aTTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGtt                                            >  1:567593/1‑251 (MQ=255)                                                                                                                                                                                      aaGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTc                                        <  1:707446/251‑1 (MQ=255)                                                                                                                                                                                            ttttttGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGca                                  >  1:279570/1‑251 (MQ=255)                                                                                                                                                                                                      tAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTc                                                                >  1:383333/1‑211 (MQ=255)                                                                                                                                                                                                      tAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTc                                        <  1:470082/235‑1 (MQ=255)                                                                                                                                                                                                      tAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTTCCa                        >  1:112388/1‑251 (MQ=255)                                                                                                                                                                                                        aCAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCg                                                            >  1:247625/1‑213 (MQ=255)                                                                                                                                                                                                        aCAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTGCCaaa                      >  1:473743/1‑251 (MQ=255)                                                                                                                                                                                                        aCAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTGCCaaa                      >  1:632640/1‑251 (MQ=255)                                                                                                                                                                                                         cAATATTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCa                                                                                                        >  1:478176/1‑168 (MQ=255)                                                                                                                                                                                                           atatTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTGCCAaacaa                   >  1:606454/1‑251 (MQ=255)                                                                                                                                                                                                           atatTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTGCCAaacaa                   >  1:170987/1‑251 (MQ=255)                                                                                                                                                                                                           atatTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTGCCAaacaa                   >  1:169136/1‑251 (MQ=255)                                                                                                                                                                                                             atTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACa                                                                                                                                                                                                <  1:37716/76‑1 (MQ=255)                                                                                                                                                                                                             atTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGc                                               <  1:533765/221‑1 (MQ=255)                                                                                                                                                                                                             atTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGt                                             <  1:733159/223‑1 (MQ=255)                                                                                                                                                                                                             atTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTGCCAAacaaca                 >  1:62844/1‑251 (MQ=255)                                                                                                                                                                                                             atTGTAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTGCCAAacaaca                 >  1:380525/1‑251 (MQ=255)                                                                                                                                                                                                               tGTAGGCTTCGCCCCGTTTCCCAGCGCGGAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTGCCAAACAACATc               <  1:322189/251‑1 (MQ=255)                                                                                                                                                                                                                 tAGGCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTGCCAAACAACATCAg             <  1:605861/251‑1 (MQ=255)                                                                                                                                                                                                                   ggCTTCGCCCCGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTGCCAAACAACATCAGTg           >  1:92968/1‑251 (MQ=255)                                                                                                                                                                                                                          ccccGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTGCCAAACAACATCAGTGGAAAATa    >  1:266304/1‑251 (MQ=255)                                                                                                                                                                                                                            ccGTTACTCAGCGCGCAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAGCGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTGCCAAACAACATCAGTGGAAAATAgg  >  1:441359/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             ATTTCCGCCTTGCTTGCGGTGTCCATCGACATCGGCGACAAAATGGCGAAAGAGCTTAAATCCTACGGCGCCAATATTCTGGTGGAGCCCGCCAGCAGCGCCATTTTGCCTGATGAAGTGAGCCGTAATAATTCTCTCGCCACGCAAGATTTTTTGGACGAAAAAGAACTGCCGAACATTAAAGACATTTTTTGGCGTAACAATATTGTAGGCTTCGCCCCGTTACTCAGCGCACAAGTCAAAGCCGATGGACCAAACGGCAAGGCGCAAGACATCGACATTCTCGGCACGTTTTTTGATCATCAAATCGCCGTGCCGGATGAAGACGATTACCACACCGGGCAAAAAATCATCAACCCTTATTGGCAGGTGGAAGGCGAATGGGTGAACGATGCCACCGATGATTTCAGCGAACAGGTTCCTGCGTTACTCGGCGCACAACTTGCCAAACAACATCAGTGGAAAATAGG  >  1/1687805‑1688274

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