Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,630,210	T→G	L169R (CTC→CGC)	unknown →	Oxygen‑independent coproporphyrinogen‑III oxidase‑like protein YggW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,630,210	0	T	G	100.0%	178.1 / NA	62	L169R (CTC→CGC)	unknown	Oxygen‑independent coproporphyrinogen‑III oxidase‑like protein YggW
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (27/35);  total (27/35)

CATCGCGTATTTGTTAGCGCGGATAGAACAACACATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACTCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAACTGTGGGATATTTTCGAGCAAGGGCATCAATTACTCACTGCCGCCGGTTATCAACAATACGAAACCTC  >  1/1629965‑1630453                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    cATCGCGTATTTGTTAGCGCGGATAGAACAACACATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGtt                                                                                                                                                                                                                                                <  1:637223/251‑1 (MQ=255)      cgTATTTGTTAGCGCGGATAGAACAACACATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAAt                                                                                                                                                                                                                                           <  1:238924/251‑1 (MQ=255)        tATTTGTTAGCGCGGATAGAACAACACATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCACTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATc                                                                                                                                                                                                                                          <  1:401544/250‑1 (MQ=255)            tGTTAGCGCGGATAGAACAACACATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGAt                                                                                                                                                                                                                                     >  1:601924/1‑251 (MQ=255)              ttAGCGCGGATAGAACAACACATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGAtt                                                                                                                                                                                                                                    <  1:836542/250‑1 (MQ=255)                     ggATAGAACAACACATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCa                                                                                                                                                                                                                             >  1:535338/1‑250 (MQ=255)                      gATAGAACAACACATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACg                                                                                                                                                                                                                           >  1:637069/1‑251 (MQ=255)                        tAGAACAACACATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGt                                                                                                                                                                                                                         >  1:850879/1‑251 (MQ=255)                          gAACAACACATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATACCGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGttt                                                                                                                                                                                                                       >  1:424481/1‑251 (MQ=255)                          gAACAACACATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCaa                                                                                                                                                                                                                                            >  1:653818/1‑230 (MQ=255)                          gAACAACACATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATg                                                                                                                                                                                                                               >  1:189372/1‑243 (MQ=255)                          gAACAACACATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGCTTCGCATCTACGCAGATTCAATCTGGATTTAATGCACGGttt                                                                                                                                                                                                                       >  1:555270/1‑251 (MQ=255)                               acacATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACgg                                                                                                                                                                                                                          >  1:718238/1‑243 (MQ=255)                                      ccTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAAc                                                                                                                                                                                                           <  1:450718/251‑1 (MQ=255)                                      ccTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGATATGCGCAGGCGGGCGAGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCCAAGAAGCCAAAAGAGCGGTGGGTTTTGCTGGCGGTTCGCATCTAAGCAGTTTCAATCTGGATTTAATg                                                                                                                                                                                                                               >  1:188116/1‑231 (MQ=255)                                             cgcCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAAc                                                                                                                                                                                                    <  1:384847/251‑1 (MQ=255)                                             cgcCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAAc                                                                                                                                                                                                    <  1:896184/251‑1 (MQ=255)                                                     atCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCt                                                                                                                                                                                            <  1:20862/251‑1 (MQ=255)                                                               ccTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTa                                                                                                                                                                                  >  1:864643/1‑251 (MQ=255)                                                               ccTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTa                                                                                                                                                                                  >  1:270098/1‑251 (MQ=255)                                                                    gAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGcc                                                                                                                                                                              <  1:896843/250‑1 (MQ=255)                                                                    gAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGcc                                                                                                                                                                              <  1:351356/250‑1 (MQ=255)                                                                       gCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAg                                                                                                                                                                           >  1:779277/1‑250 (MQ=255)                                                                        ccAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGc                                                                                                                                                                          <  1:47968/250‑1 (MQ=255)                                                                           aTCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCAtt                                                                                                                                                                      >  1:512277/1‑251 (MQ=255)                                                                             cccGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTg                                                                                                                                                                     <  1:285254/250‑1 (MQ=255)                                                                               cGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAACCGAGGCGCTGGACGATTTACGCCAAGCCATTGCAc                                                                                                                                                                  <  1:86952/251‑1 (MQ=255)                                                                                    cggcggAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGgcg                                                                                                                                                             <  1:830300/251‑1 (MQ=255)                                                                                    cggcggAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGgcg                                                                                                                                                             >  1:582589/1‑251 (MQ=255)                                                                                     ggcggAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGAc                                                                                                                                                                                        <  1:52713/223‑1 (MQ=255)                                                                                               gcgcGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCAtt                                                                                                                                                   <  1:625348/250‑1 (MQ=255)                                                                                                    ttttCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCCTTTGTCtt                                                                                                                                             >  1:694253/1‑251 (MQ=255)                                                                                                     tttCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATc                                                                                                                                                                                                                                          >  1:170378/1‑157 (MQ=255)                                                                                                          cgGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCa                                                                                                                                                                                                              >  1:651918/1‑180 (MQ=255)                                                                                                          aaactcagcggaGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATc                                                                                                                                       <  1:422656/240‑1 (MQ=255)                                                                                                            gTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGcc                                                                                                                                                                              >  1:618896/1‑210 (MQ=255)                                                                                                            gTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGcc                                                                                                                                                                              >  1:682243/1‑210 (MQ=255)                                                                                                             ttATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATCCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGc                                                                                                                                    <  1:437585/251‑1 (MQ=255)                                                                                                              tATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCt                                                                                                                                   <  1:31419/251‑1 (MQ=255)                                                                                                                    cAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATC‑ACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCAt                                                                                                                             <  1:631841/250‑1 (MQ=255)                                                                                                                          ggcgTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATcgaac                                                                                                                        <  1:610146/250‑1 (MQ=255)                                                                                                                                aaTCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAacac                                                                                                                 <  1:898330/251‑1 (MQ=255)                                                                                                                                 aTCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACAcc                                                                                                                <  1:400916/251‑1 (MQ=255)                                                                                                                                         tCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACg                                                                                                                                                                                                   >  1:649729/1‑160 (MQ=255)                                                                                                                                                   ttCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCAc                                                                                                                                                                  >  1:822259/1‑183 (MQ=255)                                                                                                                                                   ttCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAacac                                                                                                                 <  1:272065/232‑1 (MQ=255)                                                                                                                                                                     aaaaaTTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCcgacga                                                                             <  1:365937/250‑1 (MQ=255)                                                                                                                                                                             cAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAGCGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCGCCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAAc                                                                     <  1:372010/250‑1 (MQ=255)                                                                                                                                                                                 aGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATcgaac                                                                                                                        <  1:622414/195‑1 (MQ=255)                                                                                                                                                                                           cATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAACTGTGGGATATTTTCg                                                      >  1:862680/1‑251 (MQ=255)                                                                                                                                                                                             ttCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCtt                                                                                                      >  1:518728/1‑201 (MQ=255)                                                                                                                                                                                               cATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAACTGTGGGATATTTTCGAGCa                                                  >  1:130691/1‑251 (MQ=255)                                                                                                                                                                                               cATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAACTGTGGGATATTTTCGAGCa                                                  <  1:250278/251‑1 (MQ=255)                                                                                                                                                                                                              gCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAACTGTGGGATATTTTCGAGCAAggg                                              <  1:665584/240‑1 (MQ=255)                                                                                                                                                                                                                              gTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAACTGTGGGATATTTTCGAGCAAGGGCATCAATTACTCACTGCCGCCGGTTa                    <  1:78067/250‑1 (MQ=255)                                                                                                                                                                                                                              gTTTTGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAACTGTGGGATATTTTCGAGCAAGGGCATCAATTACTCACTGCCGCCGGTTAt                   >  1:494891/1‑251 (MQ=255)                                                                                                                                                                                                                                  tGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGATGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAACTGTGGGATATTTTCGAGCAAGGGCATCAATTACTCACTGCCGCCGGTTATcaac               <  1:423376/251‑1 (MQ=255)                                                                                                                                                                                                                                  tGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAACTGTGGGATATTTTTGAGCAAGGGCATCAATTACTCACTGCCGCCGGTTATcaac               <  1:507305/251‑1 (MQ=255)                                                                                                                                                                                                                                  tGCCGGCGTTTCGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCCCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAACTGTGGGATATTTTCGAGCAAGGGCATCAATTACTCACTGCCGCCGGTTATcaac               <  1:62540/251‑1 (MQ=255)                                                                                                                                                                                                                                             cGCATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAACTGTGGGATATTTTCGAGCAAGGGCATCAAGTACTCGCTGCCGCCGGTTATCAACAATACGAAAcc    <  1:507838/251‑1 (MQ=255)                                                                                                                                                                                                                                               cATCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAACCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAACTGTGGGATATTTTCGAGCAAGGGCATCAATTACTCACTGCCGCCGGTTATCAACAATACGAAACctc  <  1:448419/251‑1 (MQ=255)                                                                                                                                                                                                                                                 tCTACGCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAACTGTGGGATATTTTCGAGCAAGGGCATCAATTACTCACTGCCGCCGGTTATCAACAATACGAAACct   >  1:589726/1‑248 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    CATCGCGTATTTGTTAGCGCGGATAGAACAACACATTCCTTTTGCGCCGAATATCGAAATCACCTTGGAAGCCAATCCCGGCACGGCGGAAGCGGCGCGTTTTCGCGGTTATGCGCAGGCGGGCGTGAATCGAATTTCCATGGGCATTCAAAGTTTTTCCGATGAAAAATTGCAAAAGCTGGGGCGCATTCATAACGCCGAAGAAGCCAAAAGTGCGGTGGGTTTTGCCGGCGTTTCGCATCTACTCAGTTTCAATCTGGATTTAATGCACGGTTTGCCGAATCAAACCTTAAACGAGGCGCTGGACGATTTACGCCAAGCCATTGCACTGGCGCCGCCGCATTTGTCTTGGTATCAGCTCACCATCGAACCGAACACCATGTTTGCTTATCGCCCGCCGACGCTGCCCGACGATGACGAACTGTGGGATATTTTCGAGCAAGGGCATCAATTACTCACTGCCGCCGGTTATCAACAATACGAAACCTC  >  1/1629965‑1630453

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