Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	787,923	T→C	G20G (GGA→GGG)	unknown ←	Toxin HigB / Protein kinase domain of HipA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	787,923	0	T	C	100.0%	189.8 / NA	69	G20G (GGA→GGG)	unknown	Toxin HigB / Protein kinase domain of HipA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (42/27);  total (42/27)

CTTCTCAGTTGGCACAGGTCCACCAAATATTGAGTTATACTCAACAAAAAGCAGAACAGTGTATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGTCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTACACTTCTTATCATTTAACGGCTCCGCATAAATTTGCCGAACAGAGCGACTTT  >  1/787700‑788156                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          cTTCTCAGTTGGCACAGGTCCACCAAATATTGAGTTATACTCAACAAAAAGCAGAACAGTGTATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTaaaa                                                                                                                                                                                                                >  1:61015/1‑251 (MQ=255)      cAGTTGGCACAGGTCCACCAAATATTGAGTTATACTCAACAAAAAGCAGAACAGTGTATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATg                                                                                                                                                                                                           >  1:664088/1‑251 (MQ=255)                      accaAATATTGAGTTATACTCAACAAAAAGCAGAACAGTGTATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAAt                                                                                                                                                                                           >  1:144853/1‑251 (MQ=255)                           atatTGAGTTATACTCAACAAAAAGCAGAACAGTGTATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTatg                                                                                                                                                                                      >  1:294567/1‑249 (MQ=255)                           atatTGAGTTATACTCAACAAAAAGCAGAACAGTGTATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAgg                                                                                                                                                                                      >  1:179260/1‑251 (MQ=255)                                         tCAACAAAAAGCAGAACAGTGTATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTg                                                                                                                                                                        <  1:282664/251‑1 (MQ=255)                                            acaaAAAGCAGAACAGTGTATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAGGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAgggg                                                                                                                                                                                    >  1:430461/1‑236 (MQ=255)                                            acaaAAAGCAGAACAGTGTATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAgggg                                                                                                                                                                                    >  1:589618/1‑236 (MQ=255)                                            acaaAAAGCAGAACAGTGTATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAggg                                                                                                                                                                                     <  1:795611/235‑1 (MQ=255)                                             caaAAAGCAGAACAGTGTATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGAttt                                                                                                                                                                    <  1:23143/251‑1 (MQ=255)                                                 aaGCAGAACAGTGTATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCgggg                                                                                                                                                                                              >  1:414770/1‑221 (MQ=255)                                                 aaGCAGAACAGTGTATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGAt                                                                                                                                                                <  1:661914/251‑1 (MQ=255)                                                           gtgtATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATAc                                                                                                                                                                                                                                 <  1:145493/176‑1 (MQ=255)                                                           gtgtATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAgggg                                                                                                                                                                                    >  1:767390/1‑221 (MQ=255)                                                             gtATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGtagagtagag                                                                                                                                                    >  1:160354/1‑251 (MQ=255)                                                             gtATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGtagagtagag                                                                                                                                                    >  1:624420/1‑251 (MQ=255)                                                                             tAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCtgtgt                                                                                                                                    >  1:763007/1‑251 (MQ=255)                                                                                  ggAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGa                                                                                                                               >  1:870520/1‑251 (MQ=255)                                                                                      gTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACa                                                                                                                             <  1:527071/249‑1 (MQ=255)                                                                                        tAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCg                                                                                                                         <  1:338710/251‑1 (MQ=255)                                                                                          ggATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAg                                                                                                                                                                                                                           >  1:127252/1‑151 (MQ=255)                                                                                          ggATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTcc                                                                                                                                         >  1:155903/1‑233 (MQ=255)                                                                                          ggATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGaa                                                                                                                       >  1:69572/1‑251 (MQ=255)                                                                                                    gggTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGt                                                                                                                                                                                                                         <  1:402307/143‑1 (MQ=255)                                                                                                          ttttAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCgg                                                                                                       <  1:469011/251‑1 (MQ=255)                                                                                                                   gTGCAGGACTGATTGAGGTTATAGATATACTTATTCACTGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTTAAAATATGGTTGGTAGCGGGGAATGTAGGGGACATTACCATTGACTTTGATAGTAGCGCAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCa                                                                                              >  1:620402/1‑251 (MQ=255)                                                                                                                   gTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTaaaaa                                                                                                                                                                                                               >  1:4577/1‑138 (MQ=255)                                                                                                                   gTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGt                                                                                                                                                                                         <  1:736110/160‑1 (MQ=255)                                                                                                                   gTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATagtagagtag                                                                                                                                                      >  1:302502/1‑195 (MQ=255)                                                                                                                   gTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTcc                                                                                                                                         >  1:527379/1‑208 (MQ=255)                                                                                                                   gTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTcc                                                                                                                                         <  1:486825/208‑1 (MQ=255)                                                                                                                   gTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGAc                                                                                                                              <  1:404773/219‑1 (MQ=255)                                                                                                                   gTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCa                                                                                              >  1:138385/1‑251 (MQ=255)                                                                                                                   gTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCa                                                                                              >  1:254215/1‑251 (MQ=255)                                                                                                                             gATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACgg                                                                                    >  1:493693/1‑251 (MQ=255)                                                                                                                                tGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGaaa                                                                                 <  1:436684/251‑1 (MQ=255)                                                                                                                                   ggTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATAc                                                                                                                                                                                                                                 >  1:543230/1‑104 (MQ=255)                                                                                                                                   ggTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGAtagtag                                                                                                                                                           <  1:5620/174‑1 (MQ=255)                                                                                                                                   ggTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATagtagagtag                                                                                                                                                      >  1:548858/1‑179 (MQ=255)                                                                                                                                       ataGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTcc                                                                                                                                         >  1:681244/1‑188 (MQ=255)                                                                                                                                       ataGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGcc                                                                                          <  1:207480/235‑1 (MQ=255)                                                                                                                                        taGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGt                                                                         >  1:335517/1‑251 (MQ=255)                                                                                                                                          gATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAg                                                                       >  1:83741/1‑251 (MQ=255)                                                                                                                                           atatACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAg                                                                       <  1:534117/250‑1 (MQ=255)                                                                                                                                              taCTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAg                                                                               >  1:81540/1‑239 (MQ=255)                                                                                                                                                   aTTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAg                                                                                                                                                                                       <  1:601714/130‑1 (MQ=255)                                                                                                                                                   aTTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCaa                                                               >  1:514116/1‑250 (MQ=255)                                                                                                                                                   aTTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCaa                                                               >  1:45473/1‑250 (MQ=255)                                                                                                                                                   aTTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCaa                                                               >  1:646902/1‑250 (MQ=255)                                                                                                                                                             tGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATg                                                                                                                    <  1:26430/187‑1 (MQ=255)                                                                                                                                                             tGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATg                                                                                                                    <  1:150421/187‑1 (MQ=255)                                                                                                                                                             tGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTaca                                                    <  1:478291/251‑1 (MQ=255)                                                                                                                                                              gTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTaca                                                    >  1:476932/1‑250 (MQ=255)                                                                                                                                                                 ttAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTACActtc                                                <  1:400346/251‑1 (MQ=255)                                                                                                                                                                   aGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTACActtctt                                              >  1:578053/1‑251 (MQ=255)                                                                                                                                                                    gCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATagtagagt                                                                                                                                                        <  1:353890/144‑1 (MQ=255)                                                                                                                                                                        aaCTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTACACTTCTTATCAt                                         <  1:97579/251‑1 (MQ=255)                                                                                                                                                                              aGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGAt                                                                                                                                                                      >  1:228482/1‑120 (MQ=255)                                                                                                                                                                                gCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTACACTTCTTATCATTTAAc                                    <  1:433094/248‑1 (MQ=255)                                                                                                                                                                                        aTTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTACACTTCTTATCATTTAACGGCTCCGCATa                         >  1:475696/1‑251 (MQ=255)                                                                                                                                                                                          tCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTACACTTCTTATCATTTAACGGCTCCGCATaaa                       <  1:552250/251‑1 (MQ=255)                                                                                                                                                                                                 cGCAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTACACTTCTTATCATTTAACGGCTCCGCATAAATTTGCCg                <  1:597276/251‑1 (MQ=255)                                                                                                                                                                                                   cAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTACACTTCTTATCATTTAAc                                    <  1:480829/229‑1 (MQ=255)                                                                                                                                                                                                   cAATTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTACACTTCTTATCATTTAACGGCTCCGCATAAATTTGCCGaa              >  1:585328/1‑251 (MQ=255)                                                                                                                                                                                                     aTTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGAtagt                                                                                                                                                             >  1:266527/1‑106 (MQ=255)                                                                                                                                                                                                     aTTATGGGAAAGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTACACTTCTTATCATTTAACGGCTCCGCATAAATTTGCCGAACa            >  1:603237/1‑251 (MQ=255)                                                                                                                                                                                                               agactAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTACACTTCTTATCATTTAACGGCTCCGCATAAATTTGCCGAACAGAGCGACttt  >  1:841452/5‑251 (MQ=255)                                                                                                                                                                                                               aGATTAAGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGt                                                                                                                                                                                         >  1:859281/1‑68 (MQ=255)                                                                                                                                                                                                                    aaGTTGCTCCTGCCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAg                                                                       >  1:558886/1‑177 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          CTTCTCAGTTGGCACAGGTCCACCAAATATTGAGTTATACTCAACAAAAAGCAGAACAGTGTATTAATGAAAAACGTAATTGGAGGTTAGGATTAATCGGGGTTATTTTAGCCGGTGCAGGACTGATTGAGGTTATAGATATACTTATTCAATGGTTGTCTTAGCCGAACTTCAGGCAAAGCCATTCGTTCTCGCAATTATGGGAAAGATTAAGTTGCTCCTGTCCGACGATACCGATAGGTAATTTAAAAATATGGCTGGTAGCGGGGAATGTAGGGGACATTACCATTGATTTTGATAGTAGAGTAGAGCCGTTTTTTCCTGTGTGCCGACAAGCGAAATGAGAAAATCATCGGTTTCTTCCATGCCTAACGGAAAAGTGCGGTAGTTTTTCAACACGTTTTACACTTCTTATCATTTAACGGCTCCGCATAAATTTGCCGAACAGAGCGACTTT  >  1/787700‑788156

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