Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	889,780	A→C	D295A (GAC→GCC)	unknown →	Transcription‑repair coupling factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	889,780	0	A	C	100.0%	216.1 / NA	72	D295A (GAC→GCC)	unknown	Transcription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (39/33);  total (39/33)

AATTAATTTATTGCCGGCACATGAATTTCCGACGGATAGCACGGGTATTGAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTACTACCGCCGAGCAAAATCAAGACCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGTGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCATTCAATCTCAACAAAAAGAACC  >  1/889535‑890020                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 aatcaattTATTGCCGGCACATGAATTTCCGACGGATAGCACGGGTATTGAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAg                                                                                                                                                                                                                                             >  1:788486/5‑251 (MQ=255) aatcaattTATTGCCGGCACATGAATTTCCGACGGATAGCACGGGTATTGAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAg                                                                                                                                                                                                                                             >  1:554248/5‑251 (MQ=255)  atcaattTATTGCCGGCACATGAATTTCCGACGGATAGCACGGGTATTGAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAgg                                                                                                                                                                                                                                            <  1:245689/248‑1 (MQ=255)      attTATTGCCGGCACATGAATTTCCGACGGATAGCACGGGTATTGAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAg                                                                                                                                                                                                                                        >  1:679804/1‑251 (MQ=255)                  cacaTGAATTTCCGACGGATAGCACGGGTATTGAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTacac                                                                                                                                                                                                                             <  1:391504/250‑1 (MQ=255)                  cacaTGAATTTCCGACGGATAGCACGGGTATTGAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACAcc                                                                                                                                                                                                                            <  1:566027/251‑1 (MQ=255)                            tCCGACGGATAGCACGGGTATTGAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCacaac                                                                                                                                                                                                                  <  1:457468/251‑1 (MQ=255)                                     tAGCACGGGTATTGAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTacac                                                                                                                                                                                                                             >  1:791186/1‑231 (MQ=255)                                        cACGGGTATTGAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGaaa                                                                                                                                                                                                      <  1:780242/251‑1 (MQ=255)                                            ggTATTGAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTacac                                                                                                                                                                                                                             <  1:179449/224‑1 (MQ=255)                                                  gAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCacaac                                                                                                                                                                                                                  <  1:880400/229‑1 (MQ=255)                                                  gAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAg                                                                                                                                                                                                                >  1:691001/1‑231 (MQ=255)                                                  gAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAg                                                                                                                                                                                                                <  1:213304/231‑1 (MQ=255)                                                              aGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATg                                                                                                                                                                                  <  1:26194/249‑1 (MQ=255)                                                                ccAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGcc                                                                                                                                                                              <  1:267099/251‑1 (MQ=255)                                                                   aTTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCgt                                                                                                                                                                           >  1:412985/1‑251 (MQ=255)                                                                    tttCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCgtt                                                                                                                                                                          >  1:232758/1‑251 (MQ=255)                                                                                                 ggATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCa                                                                                                                                             >  1:183042/1‑251 (MQ=255)                                                                                                   aTCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCAtt                                                                                                                                           >  1:370321/1‑251 (MQ=255)                                                                                                             caTTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGa                                                                                                                                 >  1:706786/1‑251 (MQ=255)                                                                                                               tttATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAg                                                                                                                                                                                                                <  1:240075/170‑1 (MQ=255)                                                                                                               tttATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCgttgtt                                                                                                                                                                       >  1:622465/1‑211 (MQ=255)                                                                                                               tttATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAAt                                                                                                                               >  1:317913/1‑251 (MQ=255)                                                                                                               tttATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAAt                                                                                                                               >  1:546969/1‑251 (MQ=255)                                                                                                                         ggTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACtt                                                                                                                                                       <  1:894355/217‑1 (MQ=255)                                                                                                                            cAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGaaa                                                                                                                   >  1:236026/1‑250 (MQ=255)                                                                                                                            cAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGaaa                                                                                                                   >  1:714379/1‑250 (MQ=255)                                                                                                                               cAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAAc                                                                                                                  <  1:719508/248‑1 (MQ=255)                                                                                                                               cAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAg                                                                                                                <  1:212377/250‑1 (MQ=255)                                                                                                                                 aaGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTAt                                                                                                             <  1:497313/251‑1 (MQ=255)                                                                                                                                    gCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATcc                                                                                                           <  1:37576/250‑1 (MQ=255)                                                                                                                                           ggctgCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATtg                                                                                                   >  1:725690/1‑251 (MQ=255)                                                                                                                                            gctgCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATtgatga                                                                                                                                      >  1:866171/1‑215 (MQ=255)                                                                                                                                            gctgCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATtg                                                                                                   <  1:750674/250‑1 (MQ=255)                                                                                                                                                  ggCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGccc                                                                                                                                                                             >  1:459443/1‑170 (MQ=255)                                                                                                                                                   gCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGt                                                                                                               >  1:901495/1‑231 (MQ=255)                                                                                                                                                   gCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTaaaa                                                                                           >  1:594994/1‑251 (MQ=255)                                                                                                                                                        gAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGa                                                                                      >  1:632586/1‑251 (MQ=255)                                                                                                                                                        gAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACC‑TGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTacac                                                                                                                                                                                                                             <  1:459771/115‑1 (MQ=255)                                                                                                                                                          atatTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGaaa                                                                                    >  1:234191/1‑251 (MQ=255)                                                                                                                                                               ggCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAg                                                                                                                                    <  1:233758/198‑1 (MQ=255)                                                                                                                                                                     cGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGgtcg                                                                         <  1:522322/251‑1 (MQ=255)                                                                                                                                                                         gTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCg                                                                                                                                                                   >  1:684342/1‑157 (MQ=255)                                                                                                                                                                         gTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGgtcgt                                                                        >  1:505390/1‑248 (MQ=255)                                                                                                                                                                          ttttttGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTgc                                                                    <  1:526903/251‑1 (MQ=255)                                                                                                                                                                          ttttttGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTgc                                                                    <  1:61613/251‑1 (MQ=255)                                                                                                                                                                                     aaTGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGc                                                                                                        >  1:8395/1‑204 (MQ=255)                                                                                                                                                                                       tGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAAt                                                       <  1:649680/251‑1 (MQ=255)                                                                                                                                                                                                      attattTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAAc                                                                                                                  >  1:604405/1‑177 (MQ=255)                                                                                                                                                                                                      attattTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATtg                                                                                                   <  1:104828/192‑1 (MQ=255)                                                                                                                                                                                                      attattTTCCGGTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCg                                        >  1:183182/1‑251 (MQ=255)                                                                                                                                                                                                                 gTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAg                                                                                                                                    >  1:146177/1‑148 (MQ=255)                                                                                                                                                                                                                 gTAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGtt                             >  1:153542/1‑251 (MQ=255)                                                                                                                                                                                                                  tAAATACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGtt                             <  1:436074/250‑1 (MQ=255)                                                                                                                                                                                                                    aaTACGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGAc                          <  1:198742/251‑1 (MQ=255)                                                                                                                                                                                                                       aCGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAg                                                                                          <  1:704890/184‑1 (MQ=255)                                                                                                                                                                                                                       aCGCTGTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCAt                       >  1:871373/1‑251 (MQ=255)                                                                                                                                                                                                                            gTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATtga                                                                                                                                         >  1:185522/1‑132 (MQ=255)                                                                                                                                                                                                                            gTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGgt                                                                           >  1:285974/1‑194 (MQ=255)                                                                                                                                                                                                                            gTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCATTc                     <  1:358604/248‑1 (MQ=255)                                                                                                                                                                                                                            gTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCATTc                     >  1:272325/1‑248 (MQ=255)                                                                                                                                                                                                                            gTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCATTCAAt                  >  1:765487/1‑251 (MQ=255)                                                                                                                                                                                                                            gTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCATTCAAt                  >  1:830780/1‑251 (MQ=255)                                                                                                                                                                                                                            gTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCATTCAAt                  >  1:8026/1‑251 (MQ=255)                                                                                                                                                                                                                            gTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCATTCAAt                  >  1:545798/1‑251 (MQ=255)                                                                                                                                                                                                                            gTTTATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCATTCAAt                  >  1:586051/1‑251 (MQ=255)                                                                                                                                                                                                                             tttATTACCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCATTCAAt                  <  1:137290/250‑1 (MQ=255)                                                                                                                                                                                                                               tattacCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCATTCAAtctc               <  1:807908/248‑1 (MQ=255)                                                                                                                                                                                                                                attacCGCCGAGCAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCATTCAATCTca              >  1:414345/3‑251 (MQ=255)                                                                                                                                                                                                                                            cAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCATTCAATCTCAACAAAAAGAAcc  <  1:618440/251‑1 (MQ=255)                                                                                                                                                                                                                                            cAAAATCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGCGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCATTCAATCTCAACAAAAAGAAcc  <  1:856284/251‑1 (MQ=255)                                                                                                                                                                                                                                              aaaTCAAGCCCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAAc                                                                                                                  <  1:447898/137‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 AATTAATTTATTGCCGGCACATGAATTTCCGACGGATAGCACGGGTATTGAATTTTTCCGTAGCCAATTTCGGGCGGAATTCGGGGAAATTCGACGGGATCCGGAACACATTTATCAACAGGTCAGCAAAGGCACGCTGGCTGCGGGCATTGAATATTGGCAGCCGCTGTTTTTTGAGCAAATGGCAACCTTGTTCAATTATTTTCCGGTAAATACGCTGTTTACTACCGCCGAGCAAAATCAAGACCAAGGCGAGCGTTTTTACACCGATGCACAACAGCGTTTTGAAAGTCGTCGCGTGGATCCTATGCGCCCGTTGTTACCGCCCGATAAACTTTGGTTGCGCATTGATGAAGTGAATCGTCACTTGAAACAGTATCCCGCTATTGTGCTAAAAGCGGAAAAAGTGCGGTCGTCTGCGCGGCAGAAAAATTTAGTGATTGCGTCGCTGCCACCGTTGACCATTCAATCTCAACAAAAAGAACC  >  1/889535‑890020

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