Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,483,690	A→C	L213L (TTA→CTC)	unknown →	Type II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,483,690	0	A	C	100.0%	218.3 / NA	76	L213L (TTA→CTC)	unknown	Type II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (34/42);  total (34/42)

GTAGTCCACTAGTGGACTCCCGACTACCTGATGGCAGCCGTTTAAATGTTGTTATTGCACCGATTGCATTAGATGGCACATCCATTTCTATTCGTAAATTCAGTAAGAATAAAAAAACATTGCAAGAGCTTGTAAATTTTGGTTCTATGACTCGCGAAATGGCGAATTTCTTAATTATTGCTGCCAGATCCCGGGTAAATATAATTGTGTCAGGTGGTACAGGATCTGGTAAAACAACCTTACTTAATGCGCTATCTAACTATATTTCCCATACCGAACGTGTCATCACACTAGAAGATACTGCCGAGCTGCGTTTGGAACAGCCTCATGTTGTACGCCTGGAAACGCGTTTAGCCGGCGTTGAGCATACAGGCGAAGTTACTATGCAAGACTTAGTAATAAATGCATTGCGTATGCGTCCTGAGCGTATTATTGTGGGCGAATGCCGTGGTGGAGAAGCATTCCAAATGTTACAAG  >  1/1483449‑1483925                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            gTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATgc                                                                                                                                                                                                                                     <  1:635676/250‑1 (MQ=255)    gTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTa                                                                                                                                                                                                                                 >  1:40588/1‑251 (MQ=255)    gTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTa                                                                                                                                                                                                                                 >  1:601913/1‑251 (MQ=255)            gTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTaatta                                                                                                                                                                                                                         >  1:408684/1‑248 (MQ=255)             tGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTTGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTaatta                                                                                                                                                                                                                         <  1:448071/250‑4 (MQ=255)             tGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTaatta                                                                                                                                                                                                                         <  1:697782/250‑4 (MQ=255)             tGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTaatta                                                                                                                                                                                                                         <  1:422294/250‑4 (MQ=255)               gACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATtata                                                                                                                                                                                                                       <  1:2897/250‑1 (MQ=255)                      gACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCAt                                                                                                                                                                                                               >  1:249121/1‑251 (MQ=255)                        cTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATAc                                                                                                                                                                                                             <  1:235047/251‑1 (MQ=255)                          acccgATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTaatta                                                                                                                                                                                                                         >  1:640523/5‑234 (MQ=255)                           cccgATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGa                                                                                                                                                                                                          >  1:761072/4‑251 (MQ=255)                             cgATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAAc                                                                                                                                                                                                        <  1:268674/250‑1 (MQ=255)                              gATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCt                                                                                                                                                                                                                              >  1:50244/1‑228 (MQ=255)                                       cGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATcac                                                                                                                                                                                              <  1:418083/251‑1 (MQ=255)                                         tttAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATcaca                                                                                                                                                                                             <  1:516466/250‑1 (MQ=255)                                               tgttgtAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTagaag                                                                                                                                                                                      <  1:489297/251‑1 (MQ=255)                                                 ttgtaATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTaatta                                                                                                                                                                                                                         <  1:571432/214‑4 (MQ=255)                                                 ttgtaATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGAt                                                                                                                                                                                    <  1:544612/251‑1 (MQ=255)                                                   gtaaTTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATgc                                                                                                                                                                                                                                     <  1:619802/197‑1 (MQ=255)                                                   gtaaTTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCt                                                                                                                                                                                                                              >  1:731078/4‑207 (MQ=255)                                                    taaTTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATAc                                                                                                                                                                                  <  1:521795/248‑1 (MQ=255)                                                    taaTTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATAc                                                                                                                                                                                  <  1:401642/248‑1 (MQ=255)                                                                  gcgttAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTagaag                                                                                                                                                                                      >  1:595869/4‑232 (MQ=255)                                                                  gcgttAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGttt                                                                                                                                                                   >  1:437599/4‑251 (MQ=255)                                                                    gttAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAGCTCTCCTTAATGCGCTATCTaatta                                                                                                                                                                                                                         >  1:429363/2‑192 (MQ=255)                                                                       aGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGaa                                                                                                                                                               >  1:883240/1‑250 (MQ=255)                                                                         aTGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAAc                                                                                                                                                              <  1:321550/249‑1 (MQ=255)                                                                         aTGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACaa                                                                                                                                                            >  1:144450/1‑250 (MQ=255)                                                                          tggtacaTCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAAc                                                                                                                                                              >  1:652590/5‑248 (MQ=255)                                                                          tggtacaTCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTAAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGa                                                                                                                                                                >  1:333926/5‑246 (MQ=255)                                                                            gtacaTCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACaac                                                                                                                                                           >  1:277005/3‑247 (MQ=255)                                                                                aTCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATAttt                                                                                                                                                                                                                    >  1:377932/1‑188 (MQ=255)                                                                                aTCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAAc                                                                                                                                                                                                        >  1:49684/1‑200 (MQ=255)                                                                                 tccatttccattCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCAtgt                                                                                                                                                    <  1:843206/251‑1 (MQ=255)                                                                                 tccatttccattCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCAtgt                                                                                                                                                    <  1:207061/251‑1 (MQ=255)                                                                                 tccatttccattCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCAtgt                                                                                                                                                    <  1:450213/251‑1 (MQ=255)                                                                                 tccatttccattCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCAtgt                                                                                                                                                    <  1:290863/251‑1 (MQ=255)                                                                                        ccaTTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGc                                                                                                                                             <  1:508393/249‑1 (MQ=255)                                                                                                  gttCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATgc                                                                                                                                                                                                                                     <  1:340412/152‑1 (MQ=255)                                                                                                  gttCAGTAAGAGTAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAAca                                                                                                                                    >  1:384333/2‑249 (MQ=255)                                                                                                           gagtAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGcc                                                                                                                           >  1:628033/4‑250 (MQ=255)                                                                                                             gtAAAAAAAGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCgg                                                                                                                         >  1:46410/2‑250 (MQ=255)                                                                                                               aaaaaaaGCTTGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAAc                                                                                                                                                                                                        >  1:161846/1‑169 (MQ=255)                                                                                                                aaaaaagcttGCAAGAGCTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGaaa                                                                                                                                      <  1:626718/226‑1 (MQ=255)                                                                                                                                 cTTGTAAATTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGaggt                                                                                                    >  1:579715/1‑248 (MQ=255)                                                                                                                                        aTTTCGGTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATg                                                                                             <  1:815454/251‑1 (MQ=255)                                                                                                                                         tttcggTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGc                                                                                            <  1:354902/247‑1 (MQ=255)                                                                                                                                              gTTCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCAtgt                                                                                                                                                    >  1:560235/1‑190 (MQ=255)                                                                                                                                                tCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATAc                                                                                                                                                                                                             >  1:765229/1‑131 (MQ=255)                                                                                                                                                tCTATGACCCGCGATATGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCAtgttg                                                                                                                                                  <  1:656595/190‑1 (MQ=255)                                                                                                                                                           cgataTGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAAt                                                                           <  1:554324/246‑1 (MQ=255)                                                                                                                                                                    aaaTTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCactacgta                                                                 <  1:96476/250‑8 (MQ=255)                                                                                                                                                                      aTTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATg                                                               <  1:241731/251‑1 (MQ=255)                                                                                                                                                                        tttttAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCg                                                             <  1:608606/248‑1 (MQ=255)                                                                                                                                                                        tttttAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCg                                                             <  1:389156/248‑1 (MQ=255)                                                                                                                                                                           ttAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTcc                                                          <  1:480766/251‑1 (MQ=255)                                                                                                                                                                                       gcCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATg                                                               >  1:311426/1‑234 (MQ=255)                                                                                                                                                                                       gcCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCg                                                             <  1:818038/236‑1 (MQ=255)                                                                                                                                                                                       gcCAGATCTCGAGTAAATATAATTGTGTCAGGTGGCACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGtattat                                              >  1:63340/1‑251 (MQ=255)                                                                                                                                                                                                cgagTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTAc                                                                                                                                               <  1:409918/142‑1 (MQ=255)                                                                                                                                                                                                  agTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGc                                                                                                                                             >  1:722098/2‑145 (MQ=255)                                                                                                                                                                                                     aaaTATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAAc                                                                                                                                     <  1:100438/150‑1 (MQ=255)                                                                                                                                                                                                     aaaTATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGcc                                <  1:493036/251‑1 (MQ=255)                                                                                                                                                                                                     aaaTATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGc                                 >  1:889650/1‑250 (MQ=255)                                                                                                                                                                                                       atatAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCgt                              >  1:644994/1‑251 (MQ=255)                                                                                                                                                                                                         ataATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCgtgg                            <  1:144208/251‑1 (MQ=255)                                                                                                                                                                                                                      gtggtACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCAtt                >  1:829725/1‑250 (MQ=255)                                                                                                                                                                                                                         gtACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGc                                                                                            <  1:558827/171‑1 (MQ=255)                                                                                                                                                                                                                         gtACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATg                                                               >  1:794791/1‑200 (MQ=255)                                                                                                                                                                                                                         gtACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCg                                                             <  1:146158/202‑1 (MQ=255)                                                                                                                                                                                                                         gtACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGtattatcgtaggt                                       >  1:361783/1‑217 (MQ=255)                                                                                                                                                                                                                         gtACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCag            >  1:334492/1‑250 (MQ=255)                                                                                                                                                                                                                          tACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCaga           <  1:494249/251‑3 (MQ=255)                                                                                                                                                                                                                                atccggTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTAc     <  1:245260/247‑1 (MQ=255)                                                                                                                                                                                                                                   cggTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAg  <  1:205371/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            GTAGTCCACTAGTGGACTCCCGACTACCTGATGGCAGCCGTTTAAATGTTGTTATTGCACCGATTGCATTAGATGGCACATCCATTTCTATTCGTAAATTCAGTAAGAATAAAAAAACATTGCAAGAGCTTGTAAATTTTGGTTCTATGACTCGCGAAATGGCGAATTTCTTAATTATTGCTGCCAGATCCCGGGTAAATATAATTGTGTCAGGTGGTACAGGATCTGGTAAAACAACCTTACTTAATGCGCTATCTAACTATATTTCCCATACCGAACGTGTCATCACACTAGAAGATACTGCCGAGCTGCGTTTGGAACAGCCTCATGTTGTACGCCTGGAAACGCGTTTAGCCGGCGTTGAGCATACAGGCGAAGTTACTATGCAAGACTTAGTAATAAATGCATTGCGTATGCGTCCTGAGCGTATTATTGTGGGCGAATGCCGTGGTGGAGAAGCATTCCAAATGTTACAAG  >  1/1483449‑1483925

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