Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	659,761	C→T	G111R (GGG→AGG)	unknown ←	High‑affinity Fe2+/Pb2+ permease precursor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	659,761	0	C	T	100.0%	214.8 / NA	77	G111R (GGG→AGG)	unknown	High‑affinity Fe2+/Pb2+ permease precursor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (44/33);  total (44/33)

TCCATTTCGGAATTTTTGTACCCCTCATACTGCGCTTGTTGCACCAGTTTTTTCGCCGTGGCGTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCCGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGCCGATTTTCAATATCCTGCACCCATG  >  1/659520‑660003                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   gccatttcGGAATTTTTGTACCCCTCATACTGCGCTTGTTGCACCAGTTTTTTCGCCGTGGCGTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGAtttt                                                                                                                                                                                                                                           >  1:383041/2‑251 (MQ=255)    atttcGGAATTTTTGTACCCCTCATACTGCGCTTGTTGCACCAGTTTTTTCGCCGTGGCGTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGc                                                                                                                                                                                                                                        <  1:606148/251‑1 (MQ=255)           aaTTTTTGTACCCCTCATACTGCGCTTGTTGCACCAGTTTTTTCGCCGTGGCGTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTAcc                                                                                                                                                                                                                                 <  1:246074/251‑1 (MQ=255)              ttttGTACCCCTCATACTGCGCTTGTTGCACCAGTTTTTTCGCCGTGGCGTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGa                                                                                                                                                                                                                              >  1:153487/1‑251 (MQ=255)                     ccccTCATACTGCGCTTGTTGCACCAGTTTTTTCGCCGTGGCGTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCg                                                                                                                                                                                                                       >  1:208183/1‑251 (MQ=255)                      cccTCATACTGCGCTTGTTGCACCAGTTTTTTCGCCGTGGCGTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCgg                                                                                                                                                                                                                      >  1:746918/1‑251 (MQ=255)                      cccTCATACTGCGCTTGTTGCACCAGTTTTTTCGCCGTGGCGTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCgg                                                                                                                                                                                                                      >  1:650718/1‑251 (MQ=255)                      cccTCATACTGCGCTTGTTGCACCAGTTTTTTCGCCGTGGCGTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCgg                                                                                                                                                                                                                      >  1:640521/1‑251 (MQ=255)                       ccTCATACTGCGCTTGTTGCACCAGTTTTTTCGCCGTGGCGTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCgg                                                                                                                                                                                                                      >  1:230220/1‑250 (MQ=255)                        cTCATACTGCGCTTGTTGCACCAGTTTTTTCGCCGTGGCGTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGc                                                                                                                                                                                                                     <  1:653866/250‑1 (MQ=255)                                gcgcTTGTTGCACCAGTTTTTTCGCCGTGGCGTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCa                                                                                                                                                                                                            >  1:232387/1‑251 (MQ=255)                                gcgcTTGTTGCACCAGTTTTTTCGCCGTGGCGTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCa                                                                                                                                                                                                            >  1:324947/1‑251 (MQ=255)                                                               gTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTgg                                                                                                                                                                             >  1:207542/1‑251 (MQ=255)                                                                        ggCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCAttt                                                                                                                                                                    >  1:656019/1‑251 (MQ=255)                                                                         gCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTg                                                                                                                                                                   >  1:793073/1‑251 (MQ=255)                                                                                  gCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCtttt                                                                                                                                                          >  1:156655/1‑251 (MQ=255)                                                                                   cGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCttttt                                                                                                                                                         >  1:109842/1‑251 (MQ=255)                                                                                           gCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGga                                                                                                                                                 >  1:811208/1‑251 (MQ=255)                                                                                              tGGGCTAATAG‑TCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAg                                                                                                                                             <  1:883424/251‑1 (MQ=255)                                                                                                ggCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCGGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGt                                                                                                                                            >  1:842221/1‑251 (MQ=255)                                                                                                ggCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCt                                                                                                                                                                                                                                       <  1:845211/160‑1 (MQ=255)                                                                                                ggCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGc                                                                                                                                                                            >  1:696798/1‑219 (MQ=255)                                                                                                ggCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGgaga                                                                                                                                               <  1:217232/248‑1 (MQ=255)                                                                                                ggCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAg                                                                                                                                             <  1:777607/250‑1 (MQ=255)                                                                                                         aTCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGaa                                                                                                                                   >  1:652057/1‑251 (MQ=255)                                                                                                                   ggtttAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCtt                                                                                                                          >  1:575782/3‑250 (MQ=255)                                                                                                                                           gcggcgATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAAt                                                                                                  <  1:184852/250‑1 (MQ=255)                                                                                                                                              gcgATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACg                                                                                               >  1:352723/1‑250 (MQ=255)                                                                                                                                                       tggtttTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCAc                                                                                     <  1:676098/251‑1 (MQ=255)                                                                                                                                                        ggtttTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCa                                                                                      >  1:91514/1‑249 (MQ=255)                                                                                                                                                            ttCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCg                                                                                >  1:17351/1‑251 (MQ=255)                                                                                                                                                               aTACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGttt                                                                             <  1:335456/251‑1 (MQ=255)                                                                                                                                                                    cGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCaaa                                                                                                    <  1:333939/223‑1 (MQ=255)                                                                                                                                                                        aTGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAAt                                                                                                  >  1:244460/1‑221 (MQ=255)                                                                                                                                                                        aTGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGcatcat                                                                    >  1:562143/1‑251 (MQ=255)                                                                                                                                                                        aTGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTCGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATg                                                                                                                                       <  1:879975/184‑1 (MQ=255)                                                                                                                                                                             gACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGtt                                                               >  1:787246/1‑251 (MQ=255)                                                                                                                                                                              aCCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGttt                                                              <  1:37214/251‑1 (MQ=255)                                                                                                                                                                              aCCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGttt                                                              <  1:367279/251‑1 (MQ=255)                                                                                                                                                                                  gAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTa                                                                                                                                          >  1:181976/1‑171 (MQ=255)                                                                                                                                                                                    ggCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTcc                                                                                                                        <  1:267916/187‑1 (MQ=255)                                                                                                                                                                                          gagttgaTGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCa                                                                                      <  1:247256/215‑1 (MQ=255)                                                                                                                                                                                            gttgaTGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCAGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTaa                                                >  1:48345/1‑251 (MQ=255)                                                                                                                                                                                                   aCCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGAt                                                                                                                    <  1:215972/176‑1 (MQ=255)                                                                                                                                                                                                      ttcttTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAgg                                      <  1:23120/251‑1 (MQ=255)                                                                                                                                                                                                        cttTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATc                                                                                                                                                                                                      <  1:173215/89‑1 (MQ=255)                                                                                                                                                                                                        cttTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTcc                                                                                                                        <  1:480799/167‑1 (MQ=255)                                                                                                                                                                                                        cttTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAAt                                                                                                  >  1:782836/1‑189 (MQ=255)                                                                                                                                                                                                        cttTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGc                                                                         >  1:187638/1‑214 (MQ=255)                                                                                                                                                                                                        cttTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGtt                                    >  1:209071/1‑251 (MQ=255)                                                                                                                                                                                                        cttTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGtt                                    >  1:677691/1‑251 (MQ=255)                                                                                                                                                                                                             aaaGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTAt                                >  1:54642/1‑250 (MQ=255)                                                                                                                                                                                                             aaaGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTAt                                >  1:26007/1‑250 (MQ=255)                                                                                                                                                                                                             aaaGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTAt                                >  1:898977/1‑250 (MQ=255)                                                                                                                                                                                                             aaaGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATc                               <  1:773400/251‑1 (MQ=255)                                                                                                                                                                                                              aaGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATcc                              <  1:398817/251‑1 (MQ=255)                                                                                                                                                                                                              aaGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATcc                              <  1:681025/251‑1 (MQ=255)                                                                                                                                                                                                                gagGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAAt                                                                                                  >  1:780516/1‑181 (MQ=255)                                                                                                                                                                                                                  gggAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGc                                                                                                                                      >  1:525822/1‑143 (MQ=255)                                                                                                                                                                                                                  gggAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAg                                              >  1:282424/1‑231 (MQ=255)                                                                                                                                                                                                                  gggAAGCACTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGc                           >  1:501558/1‑250 (MQ=255)                                                                                                                                                                                                                          cTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGa                                                                                                                     >  1:831003/1‑152 (MQ=255)                                                                                                                                                                                                                          cTTCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGCCGATTTTCa                  >  1:702721/1‑251 (MQ=255)                                                                                                                                                                                                                            tCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGCCGATTTTCAat                <  1:674518/251‑1 (MQ=255)                                                                                                                                                                                                                            tCCTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGCCGATTTTCAat                >  1:253612/1‑251 (MQ=255)                                                                                                                                                                                                                             ccTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAg                                              >  1:577113/1‑220 (MQ=255)                                                                                                                                                                                                                              cTGAAGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGCCGATTTTCAatat              <  1:58179/251‑1 (MQ=255)                                                                                                                                                                                                                                 aaGCTCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGCCGATTTTCAATATcc            <  1:266207/250‑1 (MQ=255)                                                                                                                                                                                                                                     tCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGCCGATTTTCAATATCCTGCAc       <  1:13091/251‑1 (MQ=255)                                                                                                                                                                                                                                     tCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGCCGATTTTCAATATCCTGCAc       <  1:725930/251‑1 (MQ=255)                                                                                                                                                                                                                                     tCGGCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGCCGATTTTCAATATCCTGCAc       <  1:626215/251‑1 (MQ=255)                                                                                                                                                                                                                                        gCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCAc                                                                                     <  1:529678/170‑1 (MQ=255)                                                                                                                                                                                                                                        gCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGCCGATTTTCAATATCCTGCACCCa    >  1:77233/1‑251 (MQ=255)                                                                                                                                                                                                                                        gCTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGATAATGCg                                                                                                                                     >  1:836443/1‑122 (MQ=255)                                                                                                                                                                                                                                         cTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGCCGATTTTCAATATCCTGCACCCa    <  1:531530/250‑1 (MQ=255)                                                                                                                                                                                                                                         cTTTCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGCCGATTTTCAATATCCTGCACCCa    <  1:198429/250‑1 (MQ=255)                                                                                                                                                                                                                                          tttCAGCCTGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGCCGATTTTCAATATCCTGCACCCATg  <  1:779865/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   TCCATTTCGGAATTTTTGTACCCCTCATACTGCGCTTGTTGCACCAGTTTTTTCGCCGTGGCGTAATCTTGGGCTTCATAAGCGGTAATTGCTTGGGCTAATAGATCGTCAATGGTTTTAAAACTTTTTTGCCAGTGGGCGGCGATGGTTTGGTTTTCATACACGCCATGTTGACCGGAGGCGTTGAGTTGATGACCTTCTTTTAAAGAGGGAAGCACTTCCTGAAGCTCGGCTTTCAGCCCGTCGATTTTGGCTTTTACCTGATCTTGCGGCAATCCGTCACGAATCATTTTGCGAATTTCACCGAAAGTGGCTTCCATTTGATAGCTTTTTTGTGCGGAGAAGTTAATGCGAATCGGACCTTCCAGATTTTCAAACACCTCAAAATACGCCATTTGCACTTCCGTTTGAGCATCATCGGTTTTGTTCTCGGCGTAAAGTGCGGAGGTTTTATCCAGCCGATTTTCAATATCCTGCACCCATG  >  1/659520‑660003

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