Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,349,659	G→A	A32T (GCC→ACC)	unknown →	GTP‑binding protein EngA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,349,659	0	G	A	100.0%	185.4 / NA	68	A32T (GCC→ACC)	unknown	GTP‑binding protein EngA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/35);  total (33/35)

AAATAAATTATCACGCCAAAGCCTTGTCATTTCTCATCATGGCTAAAATCCGCTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCGCCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACGTCGGCGGACATTGGCATCGCCAATTATTTACGCCAACGCCAAAAC  >  1/1349442‑1349901                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   aaataaatTATCACGCCAAAGCCTTGTCATTTCTCATCATGGCTAAAATCCGCTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAat                                                                                                                                                                                                                   >  1:651162/1‑251 (MQ=255)              cGCCAAAGCCTTGTCATTTCTCATCATGGCTAAAATCCGCTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATgg                                                                                                                                                                                                               >  1:711519/1‑242 (MQ=255)                 cAAAGCCTTGTCATTTCTCATCATGGCTAAAATCCGCTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCt                                                                                                                                                                                                   <  1:246535/251‑1 (MQ=255)                     gCCTTGTCATTTCTCATCATGGCTAAAATCCGCTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTc                                                                                                                                                                                                             <  1:148282/237‑1 (MQ=255)                         tGTCATTTCTCATCATGGCTAAAATCCGCTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAata                                                                                                                                                                                                                  <  1:595848/228‑1 (MQ=255)                         tGTCATTTCTCATCATGGCTAAAATCCGCTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGcc                                                                                                                                                                                           <  1:756732/251‑1 (MQ=255)                                   catcatGGCTAAAATCCGCTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTAt                                                                                                                                                                                 <  1:209316/251‑1 (MQ=255)                                          gCTAAAATCCGCTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCAtt                                                                                                                                                                          >  1:886551/1‑251 (MQ=255)                                          gCTAAAATCCGCTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCAtt                                                                                                                                                                          >  1:749356/1‑251 (MQ=255)                                              aaaTCCGCTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATa                                                                                                                                                                      <  1:347716/251‑1 (MQ=255)                                              aaaTCCGCTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATa                                                                                                                                                                      <  1:454175/251‑1 (MQ=255)                                                     cTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTc                                                                                                                                                                                                             <  1:420611/205‑1 (MQ=255)                                                     cTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCgg                                                                                                                                                               >  1:231153/1‑251 (MQ=255)                                                     cTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCg                                                                                                                                                                <  1:328078/250‑1 (MQ=255)                                                     cTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCACGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCgg                                                                                                                                                               >  1:655266/1‑251 (MQ=255)                                                           tgcctgccGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGa                                                                                                                                                         >  1:248298/1‑251 (MQ=255)                                                                    tttttCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGga                                                                                                                                                >  1:849650/1‑251 (MQ=255)                                                                      tttCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGgaag                                                                                                                                              <  1:425296/251‑1 (MQ=255)                                                                            gTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGgtg                                                                                                                                        >  1:488501/1‑251 (MQ=255)                                                                             tGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGgtgt                                                                                                                                       <  1:410311/251‑1 (MQ=255)                                                                             tGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGgt                                                                                                                                         <  1:617025/249‑1 (MQ=255)                                                                             tGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGgt                                                                                                                                         <  1:231717/249‑1 (MQ=255)                                                                                  atatCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTagaag                                                                                                                                  >  1:219671/1‑251 (MQ=255)                                                                                      cAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGaaaa                                                                                                                              >  1:874423/1‑251 (MQ=255)                                                                                                     acaTTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGc                                                                                                                                                                 >  1:647527/1‑201 (MQ=255)                                                                                                       atttcatttATGAGAGGTTTTTATAACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTTAATCGACTGACCCGTACCCGCGACGCCCACGTCACCGGTTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTAGCAACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATcac                                                                                                             <  1:454632/251‑3 (MQ=255)                                                                                                        tttcatttATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCg                                                                                                              >  1:571456/1‑249 (MQ=255)                                                                                                          tcatttATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCgctg                                                                                                           >  1:703919/1‑250 (MQ=255)                                                                                                          tcatttATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCgctg                                                                                                           >  1:615799/1‑250 (MQ=255)                                                                                                           catttATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGgtgt                                                                                                                                       <  1:646821/221‑1 (MQ=255)                                                                                                               tATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGt                                                                                                                                                                                                              >  1:118747/1‑146 (MQ=255)                                                                                                                      ggTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAag                                                                                              >  1:558079/1‑251 (MQ=255)                                                                                                                       gTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCg                                                                                                                                                                               >  1:612056/1‑169 (MQ=255)                                                                                                                       gTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCAc                                                                                                                                                   <  1:477247/197‑1 (MQ=255)                                                                                                                       gTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAag                                                                                              <  1:386854/250‑1 (MQ=255)                                                                                                                       gTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAag                                                                                              <  1:528090/250‑1 (MQ=255)                                                                                                                              tGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAag                                                                                              <  1:683577/243‑1 (MQ=255)                                                                                                                                aCAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGatat                                                                                    <  1:365025/251‑1 (MQ=255)                                                                                                                                 cAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTAt                                                                                                                                                            <  1:436051/178‑1 (MQ=255)                                                                                                                                 cAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAgg                                                                                                                                          <  1:360552/196‑1 (MQ=255)                                                                                                                                    cTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGgt                                                                                                                                         <  1:269171/194‑1 (MQ=255)                                                                                                                                      cccGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCg                                                                                         >  1:134975/1‑240 (MQ=255)                                                                                                                                               gCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTtggg                                                                     >  1:662048/1‑250 (MQ=255)                                                                                                                                                ccTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCGCCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAAGGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTtggt                                                                     <  1:664396/250‑1 (MQ=255)                                                                                                                                                ccTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTtggtg                                                                    <  1:117280/251‑1 (MQ=255)                                                                                                                                                  ttGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGgaag                                                                                                                                              >  1:319224/1‑175 (MQ=255)                                                                                                                                                                 aaTGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGACCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATTCCCGCGCAGGAtt                                                   >  1:766478/1‑251 (MQ=255)                                                                                                                                                                     tCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACg                                               <  1:613963/251‑1 (MQ=255)                                                                                                                                                                     tCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACg                                               <  1:82044/251‑1 (MQ=255)                                                                                                                                                                     tCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACg                                               <  1:276249/251‑1 (MQ=255)                                                                                                                                                                     tCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCCTTTTGGTGGATGCCCGCGCAGGATTAACg                                               <  1:801422/251‑1 (MQ=255)                                                                                                                                                                        gCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACGTc                                             <  1:297623/250‑1 (MQ=255)                                                                                                                                                                            aTCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGACGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACGTcggcgg                                        >  1:444534/1‑251 (MQ=255)                                                                                                                                                                                   cTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACGTCGGCGGACATTgg                                 >  1:496950/1‑251 (MQ=255)                                                                                                                                                                                   cTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGAAATTGGTCTCTTTTTTGTGGGTGCCCGCGCAAGGTTAACGACGGCGGACAtttg                                 >  1:862767/1‑249 (MQ=255)                                                                                                                                                                                     aTTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTtggtg                                                                    <  1:115266/214‑1 (MQ=255)                                                                                                                                                                                     aTTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACGTCGGCGGACATTGGCa                               <  1:720191/251‑1 (MQ=255)                                                                                                                                                                                     aTTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACGTCGGCGGACATTGGCa                               >  1:199260/1‑251 (MQ=255)                                                                                                                                                                                     aTTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACGGCGGCGGACATTGGCa                               >  1:439550/1‑251 (MQ=255)                                                                                                                                                                                     aTTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGGTCTCTTTTTTGTGGATGCCCGCGCAGGAATAACGTCCGCGGACACTGGCa                               >  1:327437/1‑251 (MQ=255)                                                                                                                                                                                            cGTCTGACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACGTCGGCGGACATTGGCATCGCCaa                        >  1:546642/1‑251 (MQ=255)                                                                                                                                                                                                 gACCCGTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACGTCGGCGGACATTGGCATCGCCAattatt                   >  1:796424/1‑251 (MQ=255)                                                                                                                                                                                                      gTACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGGTCTCTTTTTTGTGGATGCCCGCGCAGGATTAACGGCGGCGGACATTGGCATCGCCAATTATTTacg               >  1:826177/1‑250 (MQ=255)                                                                                                                                                                                                       tACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACGTCGGCGGACATTGGCATCGCCAATTATTTacgcc             <  1:23879/251‑1 (MQ=255)                                                                                                                                                                                                       tACCCGCGACGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACGTCGGCGGACATTGGCATCGCCAATTATTTacgcc             <  1:664082/251‑1 (MQ=255)                                                                                                                                                                                                               aCGCCCTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGc                                                                                                                                                     >  1:620477/1‑107 (MQ=255)                                                                                                                                                                                                                  cccTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCgc                                                                                                             >  1:12780/1‑144 (MQ=255)                                                                                                                                                                                                                  cccTCGTCACCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACGTCGGCGGACATTGGCATCGCCAATTATTTACGCCAACGCcaaaac  <  1:526915/251‑1 (MQ=255)                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   AAATAAATTATCACGCCAAAGCCTTGTCATTTCTCATCATGGCTAAAATCCGCTATAATGCCTGCCGTTTTTCAAGTGCCGATATCAATCTTCCCGTGGCACATTTCATTTATGAGAGGTTTTTATGACAACTCCCGTAGTTGCCTTGGTGGGTCGCCCGAATGTCGGCAAATCCACGCTATTCAATCGTCTGACCCGTACCCGCGACGCCCTCGTCGCCGATTTTCCCGGTTTAACCCGCGATCGCAAATATGGTCATGCCCACCTTGCCGGCCATGATTTTATCGTCATTGATACGGGCGGTATCGACGGCACGGAAGAAGGTGTAGAAGAAAAAATGGCGGAGCAATCGCTGCTCGCCATTGAAGAGGCGGATATTGTTCTCTTTTTGGTGGATGCCCGCGCAGGATTAACGTCGGCGGACATTGGCATCGCCAATTATTTACGCCAACGCCAAAAC  >  1/1349442‑1349901

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