Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,483,840	C→T	D263D (GAC→GAT)	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,840	0	C	T	100.0%	211.4 / NA	79	D263D (GAC→GAT)	unknown	Type II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/45);  total (34/45)

CGAAATGGCGAATTTCTTAATTATTGCTGCCAGATCCCGGGTAAATATAATTGTGTCAGGTGGTACAGGATCTGGTAAAACAACCTTACTTAATGCGCTATCTAACTATATTTCCCATACCGAACGTGTCATCACACTAGAAGATACTGCCGAGCTGCGTTTGGAACAGCCTCATGTTGTACGCCTGGAAACGCGTTTAGCCGGCGTTGAGCATACAGGCGAAGTTACTATGCAAGACTTAGTAATAAATGCATTGCGTATGCGTCCTGAGCGTATTATTGTGGGCGAATGCCGTGGTGGAGAAGCATTCCAAATGTTACAAGCAATGAATACCGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATAGCCCTCGCGATGCAACATCTCGCTTGGAAAGTATGGTTATGATGTCGAATGCATCCCTACCGTTAGAGGCAATTCGCCGTAATATCTCTTCAGCAGTAAACATA  >  1/1483603‑1484077                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              cgataTGGCAAATTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAAt                                                                                                                                                                                                                                   <  1:554324/246‑1 (MQ=255)          aaaTTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCactacgta                                                                                                                                                                                                                         <  1:96476/250‑8 (MQ=255)            aTTTTTTAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATg                                                                                                                                                                                                                       <  1:241731/251‑1 (MQ=255)              tttttAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCg                                                                                                                                                                                                                     <  1:389156/248‑1 (MQ=255)              tttttAATTATCGCTGCCAGATCTCGAGTAAATATAATTGTGTCAGGTGGTACAGGATCCGGTAAAACAACTCTCCTTAATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCg                                                                                                                                                                                                                     <  1:608606/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)                                           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)                                                               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)                                                                                   actctccttaATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGaacac                                                                                                                                                <  1:563851/245‑4 (MQ=255)                                                                                      ctccttaATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGaacactgg                                                                                                                                             >  1:103052/4‑245 (MQ=255)                                                                                        ccttaATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGc                                                                                                                                                                                                             <  1:679120/184‑1 (MQ=255)                                                                                        ccttaATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGaacactggt                                                                                                                                            >  1:102206/2‑243 (MQ=255)                                                                                          ttaATGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCa                                                                                                                                          >  1:519060/1‑250 (MQ=255)                                                                                             aTGCGCTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGaacac                                                                                                                                                >  1:475082/1‑238 (MQ=255)                                                                                                  cTATCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGtattat                                                                                                                                                                                                      <  1:422013/182‑1 (MQ=255)                                                                                                    aTCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGtattatcgtaggtg                                                                                                                                                                                              >  1:217916/1‑180 (MQ=255)                                                                                                    aTCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTc                                                                                                                                           <  1:471849/239‑1 (MQ=255)                                                                                                    aTCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTa                                                                                                                               >  1:188748/1‑251 (MQ=255)                                                                                                    aTCTAATTATATTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTa                                                                                                                               >  1:216418/1‑251 (MQ=255)                                                                                                       taattataTTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGt                                                                                                                            >  1:473526/5‑251 (MQ=255)                                                                                                            atatTTCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACg                                                                                                                       >  1:5915/1‑251 (MQ=255)                                                                                                                ttCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGtattatcgt                                                                                                                                                                                                   >  1:300534/1‑168 (MQ=255)                                                                                                                ttCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCa                                                                                                                                          <  1:109779/228‑1 (MQ=255)                                                                                                                ttCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTg                                                                                                                    >  1:89481/1‑250 (MQ=255)                                                                                                                ttCCCATACCGAACGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTg                                                                                                                    >  1:483751/1‑250 (MQ=255)                                                                                                                           aaCGTGTCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATac                                                                                                         <  1:363174/250‑2 (MQ=255)                                                                                                                               tgtCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCt                                                                                                     >  1:803416/1‑250 (MQ=255)                                                                                                                                gtCATCACATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGta                                                                                                                                                                                                          <  1:763185/148‑1 (MQ=255)                                                                                                                                  catcaCATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAg                                                                                                                                                          <  1:427258/194‑1 (MQ=255)                                                                                                                                  catcaCATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTcggg                                                                                                 <  1:298321/251‑3 (MQ=255)                                                                                                                                   atcaCATTAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTcggga                                                                                                <  1:847545/251‑4 (MQ=255)                                                                                                                                       cattAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTcggg                                                                                                 >  1:664584/4‑244 (MQ=255)                                                                                                                                       cattAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCa                                                                                            >  1:846900/4‑251 (MQ=255)                                                                                                                                         ttAGAAGATACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATATAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCa                                                                                            <  1:426233/248‑1 (MQ=255)                                                                                                                                              agaTACTGCTGAATTGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGt                                                                                                                            >  1:354185/1‑212 (MQ=255)                                                                                                                                                   ctgctgaattGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGaacact                                                                                                                                               <  1:718168/179‑5 (MQ=255)                                                                                                                                                   ctgctgaattGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGaacac                                                                                                                                                >  1:258161/10‑184 (MQ=255)                                                                                                                                                   ctgctgaattGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATaccc                                                                                                       <  1:113111/219‑4 (MQ=255)                                                                                                                                                   ctgctgaattGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTg                                                                                >  1:829310/10‑251 (MQ=255)                                                                                                                                                    tgctgaattGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTgg                                                                               <  1:149090/243‑1 (MQ=255)                                                                                                                                                     gctgaattGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGa                                                                              >  1:496729/8‑251 (MQ=255)                                                                                                                                                         aattGCGTTTGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGagag                                                                           <  1:526348/247‑4 (MQ=255)                                                                                                                                                                tttGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTa                                                                   <  1:383791/251‑1 (MQ=255)                                                                                                                                                                tttGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTa                                                                   <  1:452672/251‑1 (MQ=255)                                                                                                                                                                tttGGAACAACCTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTa                                                                   <  1:323908/251‑1 (MQ=255)                                                                                                                                                                         accTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTc                                                                                                                                 <  1:655922/179‑1 (MQ=255)                                                                                                                                                                         accTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACAtctc                                                                                     <  1:28837/223‑1 (MQ=255)                                                                                                                                                                         accTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTATGATGTCg                                                           <  1:475616/249‑1 (MQ=255)                                                                                                                                                                         accTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTATGATGTCGa                                                          >  1:782473/2‑251 (MQ=255)                                                                                                                                                                          ccTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGc                                                                                                                                                                                                             >  1:649693/1‑103 (MQ=255)                                                                                                                                                                          ccTCATGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGaac                                                                                                                                                  <  1:373503/162‑2 (MQ=255)                                                                                                                                                                              aTGTTGTACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGt                                                                     <  1:655968/235‑1 (MQ=255)                                                                                                                                                                                gttgtACGCCTTGAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTATGATGTCGAATGcgt                                                    >  1:407762/1‑248 (MQ=255)                                                                                                                                                                                        ccttgAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATacc                                                                                                        >  1:126127/5‑188 (MQ=255)                                                                                                                                                                                            gAAACACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAg                                                                                                                                                          >  1:163982/1‑136 (MQ=255)                                                                                                                                                                                             aaaCACGTTTAGCCGGCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTATGATGTCGAATGCGTCCCTACCGttggaa                                      >  1:41889/1‑247 (MQ=255)                                                                                                                                                                                                            gCGTTGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTATGATGt                                                             >  1:643934/1‑213 (MQ=255)                                                                                                                                                                                                                tGAGCATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattg                   <  1:117009/251‑3 (MQ=255)                                                                                                                                                                                                                    cATACAGGCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgca                 <  1:292012/249‑5 (MQ=255)                                                                                                                                                                                                                          ggCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTATGATGTCGAATGCGTCCCTACCGttgga                                       <  1:15671/221‑4 (MQ=255)                                                                                                                                                                                                                           gCGAGGTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcag          <  1:538650/249‑12 (MQ=255)                                                                                                                                                                                                                            cgaggTTACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACGCTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaa       <  1:28129/247‑15 (MQ=255)                                                                                                                                                                                                                                 ttACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaaatatc  <  1:689306/251‑20 (MQ=255)                                                                                                                                                                                                                                 ttACTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTATGATGTCGAATGCGTCCCTACCGTTGGAAGCAATTCGCCGTAatattgcatctgcagtaaatatc  <  1:215324/251‑20 (MQ=255)                                                                                                                                                                                                                                   aCTATGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGACGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAACATCTCGATTGGAGAGTATGGTTa                                                                   <  1:754580/184‑1 (MQ=255)                                                                                                                                                                                                                                      aTGCAGGATTTAGTAATAAATGCACTACGTATGCGTCCTGAGCGTATTATCGTAGGTGAATGCCGTGGTGGAGAAGCATTCCAGATGTTACAAGCGATGAACACTGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATACCCCTCGGGATGCAAc                                                                                          <  1:237478/158‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              CGAAATGGCGAATTTCTTAATTATTGCTGCCAGATCCCGGGTAAATATAATTGTGTCAGGTGGTACAGGATCTGGTAAAACAACCTTACTTAATGCGCTATCTAACTATATTTCCCATACCGAACGTGTCATCACACTAGAAGATACTGCCGAGCTGCGTTTGGAACAGCCTCATGTTGTACGCCTGGAAACGCGTTTAGCCGGCGTTGAGCATACAGGCGAAGTTACTATGCAAGACTTAGTAATAAATGCATTGCGTATGCGTCCTGAGCGTATTATTGTGGGCGAATGCCGTGGTGGAGAAGCATTCCAAATGTTACAAGCAATGAATACCGGTCACGATGGTTCTATGTCAACGTTGCACGCAAATAGCCCTCGCGATGCAACATCTCGCTTGGAAAGTATGGTTATGATGTCGAATGCATCCCTACCGTTAGAGGCAATTCGCCGTAATATCTCTTCAGCAGTAAACATA  >  1/1483603‑1484077

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