Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,481,276	A→G	I421V (ATT→GTT)	unknown →	Type II/IV secretion system secretin RcpA/CpaC, associated with Flp pilus assembly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,481,276	0	A	G	100.0%	212.9 / NA	74	I421V (ATT→GTT)	unknown	Type II/IV secretion system secretin RcpA/CpaC, associated with Flp pilus assembly
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/50);  total (24/50)

GAATATTCCGTTCTTTAATACCCGTAAAGCAAAATCTTTATTTGAGGTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAGTGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTGGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTATTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAAACATTAACATCTAACTTCTTAAAAAATAGTGGGTTCATACAATGAGAAAATTAGTTATCACGGCCTCCATTGCATTGGCCTCAATTACCGCCCCGGTAACAACATTTGCAATCAGTGAATTTGTTAACGATACGACTGTTGCACAAC  >  1/1481032‑1481500                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 gAATATTCCGTTCTTTAATACCCGTAAAGCAAAGTCTTTATTTGAAGTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTAc                                                                                                                                                                                                                             >  1:402589/1‑250 (MQ=255)  aaTATTCCGTTCTTTAATACCCGTAAAGCAAAGTCTTTATTTGAAGTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTAccc                                                                                                                                                                                                                           <  1:730686/251‑1 (MQ=255)             cTTTAATACCCGTAAAGCAAAGTCTTTATTTGAAGTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTAc                                                                                                                                                                                                                             <  1:87451/238‑1 (MQ=255)             cTTTAATACCCGTAAAGCAAAGTCTTTATTTGAAGTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGaaa                                                                                                                                                                                                                <  1:40447/251‑1 (MQ=255)                  aTACCCGTAAAGCAAAGTCTTTATTTGAAGTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAAcc                                                                                                                                                                                                           <  1:73225/251‑1 (MQ=255)                               aaagtCTTTATTTGAAGTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGaaa                                                                                                                                                                                              <  1:578680/247‑1 (MQ=255)                                  gtCTTTATTTGAAGTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGAt                                                                                                                                                                                           >  1:698302/2‑251 (MQ=255)                                    cTTTATTTGAAGTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGAttc                                                                                                                                                                                         <  1:266667/251‑1 (MQ=255)                                       tATTTGAAGTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGAttcttc                                                                                                                                                                                      <  1:523771/251‑1 (MQ=255)                                        aTTTGAAGTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTcc                                                                                                                                                                                     <  1:118996/251‑1 (MQ=255)                                          ttGAAGTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCAt                                                                                                                                                                                   <  1:367296/251‑1 (MQ=255)                                             aagTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATAc                                                                                                                                                                                 <  1:56696/248‑1 (MQ=255)                                             aagTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATAc                                                                                                                                                                                 <  1:726876/248‑1 (MQ=255)                                             aagTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATAc                                                                                                                                                                                 <  1:822528/248‑1 (MQ=255)                                                tAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGAc                                                                                                                                                                              <  1:872205/250‑1 (MQ=255)                                                tAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGAc                                                                                                                                                                              <  1:61851/250‑1 (MQ=255)                                                   cTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTc                                                                                                                                                                            <  1:301665/249‑1 (MQ=255)                                                   cTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTc                                                                                                                                                                            <  1:498393/249‑1 (MQ=255)                                                          agaAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAg                                                                                                                                                                    <  1:621188/250‑1 (MQ=255)                                                          agaAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAg                                                                                                                                                                    <  1:575743/250‑1 (MQ=255)                                                          agaAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAg                                                                                                                                                                    <  1:56193/250‑1 (MQ=255)                                                            aaaGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGaa                                                                                                                                                                  <  1:443364/250‑1 (MQ=255)                                                            aaaGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGaa                                                                                                                                                                  <  1:629169/250‑1 (MQ=255)                                                                ttttATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGt                                                                                                                                                               <  1:283383/249‑1 (MQ=255)                                                                   tattatCGGTAGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATATTGTTTACCCTGATTTTGAAAAATCCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTttatt                                                                                                                                                          <  1:890265/251‑1 (MQ=255)                                                                     ttatCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGCTGTTTACCCTGATTTTGAAAAAACTGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTAc                                                                                                                                                        <  1:148992/251‑1 (MQ=255)                                                                     ttatCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGACACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTAc                                                                                                                                                        <  1:612029/251‑1 (MQ=255)                                                                           gTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACaagac                                                                                                                                                   <  1:547811/250‑4 (MQ=255)                                                                             ggCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGaa                                                                                                                                                                  <  1:281560/233‑1 (MQ=255)                                                                             ggCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACAtt                                                                                                                                                >  1:348098/1‑251 (MQ=255)                                                                              gCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTa                                                                                                                                               >  1:418738/1‑251 (MQ=255)                                                                                ttATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAAc                                                                                                                                             <  1:19889/251‑1 (MQ=255)                                                                                 tATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACa                                                                                                                                            <  1:713673/251‑1 (MQ=255)                                                                                  aTTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACAt                                                                                                                                           >  1:343454/1‑251 (MQ=255)                                                                                  aTTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACAt                                                                                                                                           <  1:570967/251‑1 (MQ=255)                                                                                  aTTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACAt                                                                                                                                           <  1:734551/251‑1 (MQ=255)                                                                                     cAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACAt                                                                                                                                           >  1:627371/1‑248 (MQ=255)                                                                                        cTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAAc                                                                                                                                      <  1:156225/250‑1 (MQ=255)                                                                                          cAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAActt                                                                                                                                    >  1:283624/1‑250 (MQ=255)                                                                                               gACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATgg                                                                                                                                                                                                 <  1:92958/184‑1 (MQ=255)                                                                                                aCCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGa                                                                                                                                                                                                                        <  1:450934/160‑1 (MQ=255)                                                                                                 ccTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACaagac                                                                                                                                                   >  1:402423/1‑225 (MQ=255)                                                                                                   ttGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGAttcttc                                                                                                                                                                                      >  1:692825/1‑191 (MQ=255)                                                                                                                cAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATAc                                                                                                             <  1:633376/251‑1 (MQ=255)                                                                                                                cAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATAc                                                                                                             <  1:397528/251‑1 (MQ=255)                                                                                                                        ctctTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAATGAGaa                                                                                                     >  1:584154/1‑251 (MQ=255)                                                                                                                          ctTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGaa                                                                                                                                                                  <  1:254106/188‑1 (MQ=255)                                                                                                                          ctTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACAt                                                                                                                                           >  1:19574/1‑211 (MQ=255)                                                                                                                          ctTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAAc                                                                                                                                      >  1:191150/1‑216 (MQ=255)                                                                                                                          ctTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAAc                                                                                                                                      >  1:374242/1‑216 (MQ=255)                                                                                                                                     cATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTc                                                                                                                                                                            >  1:432591/1‑167 (MQ=255)                                                                                                                                     cATCCCAATTTTAGGTTCTTTCTTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATAc                                                                                                             <  1:350451/230‑1 (MQ=255)                                                                                                                                                    ttctttctTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAATGAGAAAATTAGTTATTACGGCCTCcattgcat                                                                          <  1:299605/250‑1 (MQ=255)                                                                                                                                                    ttctttctTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAATGAGAAAATTAGTTATTACGGCCTCcattgcat                                                                          <  1:390986/250‑1 (MQ=255)                                                                                                                                                    ttctttctTCCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAATGAGAAAATTAGTTATTACGGCCTCcattgcat                                                                          <  1:677226/250‑1 (MQ=255)                                                                                                                                                     tctttcttCCGTAATG‑AACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGtt                                                                                                                  >  1:893751/1‑208 (MQ=255)                                                                                                                                                        ttcttcCGTAATGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAATGAGAAAATTAGTTATTACGGCCTCCATTGCATTGGc                                                                      >  1:524493/1‑250 (MQ=255)                                                                                                                                                                  atGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACa                                                                                                                                            >  1:673363/2‑170 (MQ=255)                                                                                                                                                                  atGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAATGAGAAAATTAGTTATTACGGCCTCCATTGCATTGGCCTCAATTatcg                                                           >  1:25320/2‑248 (MQ=255)                                                                                                                                                                       acaacTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAATGAGAAAATTAGTTATTACgg                                                                                      <  1:477827/219‑1 (MQ=255)                                                                                                                                                                        caacTGAGCGCGATAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTa                                                                                                                                                                       >  1:460959/1‑137 (MQ=255)                                                                                                                                                                                    aTAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAg                                                                                                                                                                    <  1:475864/128‑1 (MQ=255)                                                                                                                                                                                    aTAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAATGAGAAAATTAGTTATTACgg                                                                                      <  1:47135/206‑1 (MQ=255)                                                                                                                                                                                    aTAAAAAAGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAATGAGAAAATTAGTTATTACGGCCTCCATTGCATTGGCCTCAATTATCGTTCCGGCAAAAACATTTg                                         <  1:277222/251‑1 (MQ=255)                                                                                                                                                                                      aaaaaaGAGTTAGTTATCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAATGAGAAAATTAGTTATTACGGCCTCCATTGCATTGGCCTCAATTATCGTTCCGGCAAAAACATTTGCa                                       >  1:743436/1‑251 (MQ=255)                                                                                                                                                                                              gttagttaTCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTaa                                                                                                                               >  1:521202/1‑155 (MQ=255)                                                                                                                                                                                                 agttaTCGTAGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACaagac                                                                                                                                                   >  1:759512/1‑129 (MQ=255)                                                                                                                                                                                                         tagCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTAc                                                                                                                                                        <  1:169512/117‑1 (MQ=255)                                                                                                                                                                                                              aCAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACAt                                                                                                                                           <  1:328686/127‑1 (MQ=255)                                                                                                                                                                                                              aCAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAATGAGAAAATTAGTTATTACGGCCTCCATTGCATTGGCCTCAATTatcg                                                           >  1:47207/1‑204 (MQ=255)                                                                                                                                                                                                               cAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAGTGAGAAAATTAGTTATTACGGCCTCCATTGCATTGGCCTCAATTATCGTTCCGGCAAAAACATTTGCAATCAATGAATTTGTTGACGATACGa              >  1:26189/1‑251 (MQ=255)                                                                                                                                                                                                                        tagtaAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAATGAGAAAATTAGTTATTACGGCCTCCATTGCATTGGCCTCAATTATCGTTCCGGCAAAAACATTTGCAATCAATGAATTTGTTGACGATACGACTGTTGcac     <  1:315531/251‑1 (MQ=255)                                                                                                                                                                                                                           taAAACCGGTAAATGAAAAAGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAATGAGAAAATTAGTTATTACGGCCTCCATTGCATTGGCCTCAATTATCGTTCCGGCAAAAACATTTGCAATCAATGAATTTGTTGACGATACGACTGTTGCACAAc  <  1:569516/251‑1 (MQ=255)                                                                                                                                                                                                                                          aaaaaGATGTTGTTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAGACATTAACATCTAACTTCTTAAAGAATAGTGGGTTCATACAATGAGAAAATTAg                                                                                               <  1:160012/143‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 GAATATTCCGTTCTTTAATACCCGTAAAGCAAAATCTTTATTTGAGGTAGCTAACGGAGAAAGTTTTATTATCGGTGGCTTATTCAGCTCAAATGACCTTGAAGGTATAAACAAAGTTCCTCTTTTAGGAGACATCCCAATTTTAGGTTCTTTCTTCCGTAGTGCAACAACTGAGCGCGATAAAAAAGAGTTAGTTATCGTGGCAACAGTAAATGTAGTAAAACCGGTAAATGAAAAAGATGTTATTTACCCTGATTTTGAAAAAACCGGAACAATGGAAAGATTCTTCCATACGACTCCGTTAAAGAATGTTTATTACAAAACATTAACATCTAACTTCTTAAAAAATAGTGGGTTCATACAATGAGAAAATTAGTTATCACGGCCTCCATTGCATTGGCCTCAATTACCGCCCCGGTAACAACATTTGCAATCAGTGAATTTGTTAACGATACGACTGTTGCACAAC  >  1/1481032‑1481500

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