Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,134,500	G→T	L399L (CTC→CTA)	unknown ←	Type IV secretion system protein VirD4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,134,500	0	G	T	100.0%	189.6 / NA	70	L399L (CTC→CTA)	unknown	Type IV secretion system protein VirD4
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (41/29);  total (41/29)

GTACGATTGCCCAGCACTTTTGAATATTCTTCCGCGTCCTCCTGATTGTCGGGGGCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATGAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCCGGCGCAAGGCGCGTAAATCAAAATCACTGGCAGAGGTAGCGTGTTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCG  >  1/1134255‑1134746                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       gtgcGATTACCTAGCACTTTAGAATATTCTTCTGCGTCCTCCTGGTTGTCGGGTGCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATtagtag                                                                                                                                                                                                                                                   >  1:153221/4‑251 (MQ=255) gtgcGATTACCTAGCACTTTAGAATATTCTTCTGCGTCCTCCTGGTTGTCGGGTGCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATtagtag                                                                                                                                                                                                                                                   >  1:597058/4‑251 (MQ=255)        tacctaGCACTTTAGAATATTCTTCTGCGTCCTCCTGGTTGTCGGGTGCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGtaataa                                                                                                                                                                                                                                             >  1:264366/6‑250 (MQ=255)        tacctaGCACTTTAGAATATTCTTCTGCGTCCTCCTGGTTGTCGGGTGCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAAc                                                                                                                                                                                                                                            <  1:509854/246‑1 (MQ=255)         acctaGCACTTTAGAATATTCTTCTGCGTCCTCCTGGTTGTCGGGTGCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGtaataa                                                                                                                                                                                                                                             >  1:601565/5‑249 (MQ=255)         acctaGCACTTTAGAATATTCTTCTGCGTCCTCCTGGTTGTCGGGGGCCTATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGACGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGtaataa                                                                                                                                                                                                                                             >  1:491335/5‑249 (MQ=255)          cctaGCACTTTAGAATATTCTTCTGCGTCCTCCTGGTTGTCGGGTGCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTg                                                                                                                                                                                                                                          <  1:281098/248‑1 (MQ=255)                aCTTTAGAATATTCTTCTGCGTCCTCCTGGTTGTCGGGTGCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATAttt                                                                                                                                                                                                                                    <  1:324203/251‑1 (MQ=255)                     agAATATTCTTCTGCGTCCTCCTGGTTGTCGGGTGCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGtaat                                                                                                                                                                                                                                               >  1:522706/2‑235 (MQ=255)                      gAATATTCTTCTGCGTCCTCCTGGTTGTCGGGTGCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAgtagta                                                                                                                                                                                                                                                 >  1:502246/1‑232 (MQ=255)                                   cgTCCTCCTGGTTGTCGGGTGCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTc                                                                                                                                                                                                                  <  1:581203/250‑1 (MQ=255)                                       ctcctGGTTGTCGGGTGCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTa                                                                                                                                                                                                             >  1:74341/1‑251 (MQ=255)                                       ctcctGGTTGTCGGGTGCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTa                                                                                                                                                                                                             >  1:797333/1‑251 (MQ=255)                                                    ggtgCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAatat                                                                                                                                                                                                 <  1:401896/247‑1 (MQ=255)                                                              gAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATa                                                                                                                                                                                      >  1:290769/1‑251 (MQ=255)                                                                         gCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGc                                                                                                                                                                               <  1:501454/247‑1 (MQ=255)                                                                            atatTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAatat                                                                                                                                                                                                 >  1:411796/1‑226 (MQ=255)                                                                            atatTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATaaa                                                                                                                                                                        >  1:387344/1‑251 (MQ=255)                                                                              attatCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACa                                                                                                                                                                      >  1:605766/1‑251 (MQ=255)                                                                                 atCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGtaata                                                                                                                                                                                                                                              >  1:75807/1‑176 (MQ=255)                                                                                 atCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGatat                                                                                                                                                                                                                                      <  1:735325/184‑1 (MQ=255)                                                                                 atCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTc                                                                                                                                                                                                      <  1:275605/216‑1 (MQ=255)                                                                                 atCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGAt                                                                                                                                                                   >  1:336499/1‑251 (MQ=255)                                                                                        aTAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGa                                                                                                                                                             <  1:520947/250‑1 (MQ=255)                                                                                           gCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTg                                                                                                                                                         >  1:321486/1‑251 (MQ=255)                                                                                                       tCTTTACCGTAACCATTAGGCTTATTGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATATCGAGGGTAGTAAACTCATCCATTAGTAGTATTAACTGATAttt                                                                                                                                                                                                                                    >  1:278555/1‑164 (MQ=255)                                                                                                       tCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATtagtag                                                                                                                                                                                                                                                   >  1:399882/1‑149 (MQ=255)                                                                                                       tCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAAc                                                                                                                                                                       <  1:876365/225‑1 (MQ=255)                                                                                                       tCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGcc                                                                                                                                                  >  1:264212/1‑246 (MQ=255)                                                                                                       tCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGAt                                                                                                                                             >  1:751193/1‑251 (MQ=255)                                                                                                       tCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGAt                                                                                                                                             >  1:556168/1‑251 (MQ=255)                                                                                                       tCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGAt                                                                                                                                             >  1:21401/1‑251 (MQ=255)                                                                                                                   ccATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGggttggtt                                                                                                                                                                                                                    <  1:868593/168‑1 (MQ=255)                                                                                                                   ccATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTg                                                                                                                                  >  1:840182/1‑250 (MQ=255)                                                                                                                   ccATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACCGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATAttg                                                                                                                                                                                               >  1:561691/1‑189 (MQ=255)                                                                                                                        aGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATAttg                                                                                                                                                                                               >  1:809451/1‑184 (MQ=255)                                                                                                                         ggCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAg                                                                                                                                                                                                            <  1:245499/170‑1 (MQ=255)                                                                                                                                     cGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGaaataaat                                                                                                               <  1:399776/251‑1 (MQ=255)                                                                                                                                      gAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGatatctt                                                                                                                                                                                                                                   <  1:239604/134‑4 (MQ=255)                                                                                                                                       aaGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCg                                                                                                             >  1:833786/1‑251 (MQ=255)                                                                                                                                       aaGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCg                                                                                                             <  1:261/251‑1 (MQ=255)                                                                                                                                         gTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTg                                                                                                           >  1:22072/1‑251 (MQ=255)                                                                                                                                             cccTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATtagtag                                                                                                                                                                                                                                                   <  1:860830/111‑1 (MQ=255)                                                                                                                                              ccTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATACTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGAtttt                                                                                                      >  1:156221/1‑251 (MQ=255)                                                                                                                                                 tgtgATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTc                                                                                                   <  1:569374/251‑1 (MQ=255)                                                                                                                                                             ataatCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAAc                                                                                                                                                                       <  1:96932/171‑1 (MQ=255)                                                                                                                                                             ataatCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTccttgcgtaaagc                                                                                       <  1:534837/251‑12 (MQ=255)                                                                                                                                                              taatCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAAAATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCa                                                                                                                                                                                  >  1:193333/1‑159 (MQ=255)                                                                                                                                                                atCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTgg                                                                                                                       <  1:51803/216‑1 (MQ=255)                                                                                                                                                                  cAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGaaa                                                                                                                    <  1:653808/217‑1 (MQ=255)                                                                                                                                                                  cAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTa                                                                                  >  1:112942/1‑251 (MQ=255)                                                                                                                                                                       ggCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCCCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATc                                                                                                                                       <  1:66129/193‑1 (MQ=255)                                                                                                                                                                        gCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATgt                                                                                                                                                                                                                              <  1:175587/105‑1 (MQ=255)                                                                                                                                                                                    aCCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCag                                                                <  1:752925/251‑1 (MQ=255)                                                                                                                                                                                     ccGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCaga                                                               <  1:844166/251‑1 (MQ=255)                                                                                                                                                                                        gCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCagaagt                                                            >  1:48175/1‑248 (MQ=255)                                                                                                                                                                                                 ccTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCagaagtcgcatatt                                                    >  1:135605/1‑239 (MQ=255)                                                                                                                                                                                                         tttATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGt                                                                                                                                                                                                              >  1:338703/1‑88 (MQ=255)                                                                                                                                                                                                           tATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCa                                         <  1:786649/251‑1 (MQ=255)                                                                                                                                                                                                            aTTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGtaataa                                                                                                                                                                                                                                             >  1:377018/1‑54 (MQ=255)                                                                                                                                                                                                            aTTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTc                                                                                                   >  1:313429/1‑192 (MQ=255)                                                                                                                                                                                                              tGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAAt                                                                                                                                        >  1:691972/1‑153 (MQ=255)                                                                                                                                                                                                                aTAATGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTg                                    >  1:608420/1‑251 (MQ=255)                                                                                                                                                                                                                   aTGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCagaagtcgcatattcaac                                                >  1:598106/1‑221 (MQ=255)                                                                                                                                                                                                                   aTGCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGaaa                                 >  1:632859/1‑251 (MQ=255)                                                                                                                                                                                                                     gCTCAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGaaaaa                               >  1:509584/1‑251 (MQ=255)                                                                                                                                                                                                                        cAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATgg                            >  1:518562/1‑251 (MQ=255)                                                                                                                                                                                                                        cAACATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATgg                            <  1:285079/251‑1 (MQ=255)                                                                                                                                                                                                                          aCATACCGACGGCAGTAAACTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGtt                          >  1:570178/1‑251 (MQ=255)                                                                                                                                                                                                                                          aaaCTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATAttg                                                                                                                                                                                               >  1:429668/1‑70 (MQ=255)                                                                                                                                                                                                                                           aaCTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGt         <  1:195089/251‑1 (MQ=255)                                                                                                                                                                                                                                            aCTCATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGt         <  1:48118/250‑1 (MQ=255)                                                                                                                                                                                                                                               cATCCATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACgg                     >  1:465579/1‑235 (MQ=255)                                                                                                                                                                                                                                                  ccATTAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCTTGCGTAAAGCGCGTAAATCAAAATCACTGGCAGAAGTCGCATATTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCg  <  1:596482/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       GTACGATTGCCCAGCACTTTTGAATATTCTTCCGCGTCCTCCTGATTGTCGGGGGCATATAGAATTTTACAGGCGATATTATCCAGGATAGCTTTAGCACCATCTTTACCGTAACCATTAGGCTTATCGGCACGAAGTTGCCCTTGTGATTGATAAATAATCAATAGGCGTAAACCATAACCGGCGATAAACCCTACCGATTTATTGATAATGCTCAACATACCGACGGCAGTAAACTCATCCATGAGTAGTAATAACTGATATTTTAATGTTGGGTTGGTTTCCGGTAGTTCGTCAATATTGTTGGCAATAACCATGCTCATAAACAGATTAGCGACTTGCGGAGCCAGGATTAAATCATTTGCCGAAATGTGGAAATAAATCGTGATTTTTTCCCGGCGCAAGGCGCGTAAATCAAAATCACTGGCAGAGGTAGCGTGTTCAACAATTTCATTTTGAAAAATGGTTAACGGACTGATAAAGGTTGAATCG  >  1/1134255‑1134746

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