Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,807,056	G→A	Y62Y (TAC→TAT)	unknown ←	Bis‑ABC ATPase Uup

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,807,056	0	G	A	100.0%	217.8 / NA	77	Y62Y (TAC→TAT)	unknown	Bis‑ABC ATPase Uup
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/41);  total (36/41)

TATTTTCAAATTGCCAGCCATTGGCATGTTCTAATTTGGCTTGAATTTGCGCAAGCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCGTATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGTTAAGCTGATGAGTGCCACAAAAAAGTCCTATAAAAATAAAGAAAAAAATAACCGCACTTTTAAAGCTGCTGATACAC  >  1/1806830‑1807300                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     tattttCAAATTGCCAGCCATTGGCATGTTCTAATTTGGCTTGAATTTGCGCAAGCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATc                                                                                                                                                                                                                              <  1:682085/251‑1 (MQ=255)   ttttCAAATTGCCAGCCATTGGCATGTTCTAATTTGGCTTGAATTTGCGCAAGCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATcac                                                                                                                                                                                                                            <  1:543527/251‑1 (MQ=255)   ttttCAAATTGCCAGCCATTGGCATGTTCTAATTTGGCTTGAATTTGCGCAAGCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATcac                                                                                                                                                                                                                            <  1:771223/251‑1 (MQ=255)                  ccaTTGGCATGTTCTAATTTGGCTTGAATTTGCGCAAGCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCaat                                                                                                                                                                                                             >  1:388539/1‑251 (MQ=255)                  ccaTTGGCATGTTCTAATTTGGCTTGAATTTGCGCAAGCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCaat                                                                                                                                                                                                             >  1:666157/1‑251 (MQ=255)                            gTTCTAATTTGGCTTGAATTTGCGCAAGCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATcac                                                                                                                                                                                                                            >  1:300978/1‑226 (MQ=255)                            gTTCTAATTTGGCTTGAATTTGCGCAAGCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAg                                                                                                                                                                                                    >  1:370569/1‑250 (MQ=255)                            gTTCTAATTTGGCTTGAATTTGCGCAAGCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAg                                                                                                                                                                                                    >  1:252825/1‑250 (MQ=255)                            gTTCTAATTTGGCTTGAATTTGCGCAAGCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAg                                                                                                                                                                                                    >  1:597818/1‑250 (MQ=255)                            gTTCTAATTTGGCTTGAATTTGCGCAAGCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAg                                                                                                                                                                                                    >  1:118482/1‑250 (MQ=255)                                                    cAAGCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATACTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTa                                                                                                                                                                           <  1:485649/251‑1 (MQ=255)                                                      aGCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATcac                                                                                                                                                                                                                            >  1:816632/1‑200 (MQ=255)                                                       gCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTccc                                                                                                                                                                                     <  1:754639/238‑1 (MQ=255)                                                          gATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGccc                                                                                                                                                                      >  1:471152/1‑250 (MQ=255)                                                           aTTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCa                                                                                                                                                                    >  1:679866/1‑251 (MQ=255)                                                               aGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTa                                                                                                                                                                 <  1:720404/250‑1 (MQ=255)                                                                 cACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTaaac                                                                                                                                                              <  1:80442/251‑1 (MQ=255)                                                                 cACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTaaac                                                                                                                                                              <  1:496944/251‑1 (MQ=255)                                                                    cTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTaaacaaa                                                                                                                                                           >  1:109908/1‑251 (MQ=255)                                                                       cTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGTCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAAcg                                                                                                                                                         <  1:440667/250‑1 (MQ=255)                                                                       cTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAAcgc                                                                                                                                                        <  1:699106/251‑1 (MQ=255)                                                                           cTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGcc                                                                                                                                                                       <  1:115004/232‑1 (MQ=255)                                                                           cTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCAc                                                                                                                                                                   <  1:28533/236‑1 (MQ=255)                                                                           cTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGcgtt                                                                                                                                                     >  1:626012/1‑250 (MQ=255)                                                                               aaTTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCgttcgttc                                                                                                                                                <  1:537273/251‑1 (MQ=255)                                                                                         cAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAatat                                                                                                                                      <  1:689118/251‑1 (MQ=255)                                                                                         cAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAatat                                                                                                                                      <  1:496362/251‑1 (MQ=255)                                                                                            ttGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTa                                                                                                                                   <  1:165342/251‑1 (MQ=255)                                                                                              gCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTaaa                                                                                                                                 >  1:351757/1‑251 (MQ=255)                                                                                                   gAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGc                                                                                                                                                                                                                                         <  1:816378/142‑1 (MQ=255)                                                                                                   gAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATcac                                                                                                                                                                                                                            >  1:642465/1‑155 (MQ=255)                                                                                                   gAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGc                                                                                                                                                                        <  1:720912/207‑1 (MQ=255)                                                                                                   gAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTc                                                                                                                                           <  1:508178/236‑1 (MQ=255)                                                                                                   gAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAatat                                                                                                                                      >  1:239219/1‑241 (MQ=255)                                                                                                   gAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTaaa                                                                                                                                 <  1:32220/246‑1 (MQ=255)                                                                                                   gAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCg                                                                                                                             >  1:99501/1‑250 (MQ=255)                                                                                                   gAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCg                                                                                                                             >  1:698942/1‑250 (MQ=255)                                                                                                   gAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCg                                                                                                                             >  1:253428/1‑250 (MQ=255)                                                                                                   gAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCg                                                                                                                             >  1:403542/1‑250 (MQ=255)                                                                                                       tatGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGa                                                                                                                        >  1:14022/1‑251 (MQ=255)                                                                                                         tgatgaTATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATc                                                                                                                      <  1:173301/251‑1 (MQ=255)                                                                                                           atgatATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGAttt                                                                                                                                                                                        >  1:662456/1‑183 (MQ=255)                                                                                                           atgatATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAatat                                                                                                                                      >  1:888908/1‑233 (MQ=255)                                                                                                           atgatATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATcc                                                                                                                     <  1:815266/250‑1 (MQ=255)                                                                                                           atgatATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCa                                                                                                                    >  1:413622/1‑251 (MQ=255)                                                                                                            tgatATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATc                                                                                                                               >  1:716956/1‑239 (MQ=255)                                                                                                               tatTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATaa                                                                                                                <  1:388750/251‑1 (MQ=255)                                                                                                               tatTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATaa                                                                                                                <  1:549024/251‑1 (MQ=255)                                                                                                                  tCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAatat                                                                                                                                      <  1:3871/226‑1 (MQ=255)                                                                                                                       cAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATaa                                                                                                                <  1:68840/243‑1 (MQ=255)                                                                                                                       cAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCAt                                                                                                        >  1:614013/1‑251 (MQ=255)                                                                                                                             aTCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTga                                                                                                  >  1:416235/1‑251 (MQ=255)                                                                                                                               ccGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCgttcgttc                                                                                                                                                <  1:504433/203‑1 (MQ=255)                                                                                                                               ccGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTgaaa                                                                                                <  1:133518/251‑1 (MQ=255)                                                                                                                                 gCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATc                                                                                                                               <  1:348048/218‑1 (MQ=255)                                                                                                                                     gATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAaagaaaga                                                                                          >  1:193338/1‑251 (MQ=255)                                                                                                                                      aTGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGaaagaaag                                                                                           <  1:542561/249‑1 (MQ=255)                                                                                                                                         cTCAATGCCCTCCGCCACATAATCAAACACATTGCCCTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCcg                                                                                                                                                                                                                                       <  1:685667/106‑1 (MQ=255)                                                                                                                                             aTGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGtt                                                                                  >  1:362397/1‑251 (MQ=255)                                                                                                                                                 ccTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGtt                                                                              <  1:240462/251‑1 (MQ=255)                                                                                                                                                     cGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATACTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGTTAAGc                                                                          <  1:451101/251‑1 (MQ=255)                                                                                                                                                       ccACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTAc                                                                                                                                                                          >  1:782334/1‑153 (MQ=255)                                                                                                                                                       ccACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTgaaagaa                                                                                             <  1:747368/230‑1 (MQ=255)                                                                                                                                                              aTCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTaaac                                                                                                                                                              <  1:673585/158‑1 (MQ=255)                                                                                                                                                                    cacaTTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGTTAAGc                                                                          <  1:212094/236‑1 (MQ=255)                                                                                                                                                                      caTTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGTTAAGCTGATGAg                                                                   <  1:26126/241‑1 (MQ=255)                                                                                                                                                                          gCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCTCCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGTTAAGCTGATGAGTGCCACaaaaa                                                        >  1:396471/1‑248 (MQ=255)                                                                                                                                                                                                    gCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTaaac                                                                                                                                                              >  1:675537/1‑120 (MQ=255)                                                                                                                                                                                                    gCTCCAAACGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGTTAAGCTGATGAGTGCCACAAAAAAGTCCTATAAAAATAAAGAAAAAAATaa                            >  1:127604/1‑250 (MQ=255)                                                                                                                                                                                                         aaaCGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGTTAAGCt                                                                         >  1:781089/1‑200 (MQ=255)                                                                                                                                                                                                           aCGGGACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCA‑‑ATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGTTAAGCTGATGAGTGCCACAAAAAAGTCCTATAAAAATAAAGAAAAAAATAACCGCACtttt                  >  1:179661/1‑251 (MQ=255)                                                                                                                                                                                                              ggACACCACCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGTTAAGCTGATGAGTGCCACAAAAAAg                                                      >  1:250753/1‑214 (MQ=255)                                                                                                                                                                                                                 caccacCAAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATc                                                                                                                               <  1:364753/138‑1 (MQ=255)                                                                                                                                                                                                                  ccccccaAGTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCCCATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGTTAAGCTGATGAGTGCCACAAAAAAGTCCTATAAAAATAAAGAAAAAAATAACCGCACTTTTAAAgc             <  1:643986/247‑1 (MQ=255)                                                                                                                                                                                                                          gTCTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGTTAAGCTGATGAGTGCCACAAAAAAGTCCTATAAAAATAAAGAAAAAAATAACCGCACTTTTAAAGCTGCTGATa     >  1:380346/1‑251 (MQ=255)                                                                                                                                                                                                                            cTTTTTCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGTTAAGCTGATGAGTGCCACAAAAAAGTCCTATAAAAATAAAGAAAAAAATAACCGCACTTTTAAAGCTGCTGATaca   <  1:327934/251‑1 (MQ=255)                                                                                                                                                                                                                             tttttCATATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGTTAAGCTGATGAGTGCCACAAAAAAGTCCTATAAAAATAAAGAAAAAAATAACCGCACTTTTAAAGCTGCTGATacac  <  1:43874/251‑1 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     TATTTTCAAATTGCCAGCCATTGGCATGTTCTAATTTGGCTTGAATTTGCGCAAGCTGATTGAGCACCTGCTCGCTGTAATTTTGCTCCAATTGCACGGAAATATGATGATATTCCTTCAATAAATCCGCCAGATGCTCAATGCCCTCCGCCACATAATCAAACACATTGCCTTCCGCATGACGCGGCGGATCCTGCTCCAAACGGGACACCACCAAGTCTTTTTCGTATTGCACTTTGCCGTCGTCCATCACCACATCGCCGGCAATAATTTTGAGCAAGGTGGATTTCCCCGCCCCATTACGCCCCACTAAACAAACGCGTTCGTTCGGTTCAATATGTAAATCCGCATGATCCAATAAAGGCGCATCACTGAAAGAAAGATATCCGTTGGTTAAGCTGATGAGTGCCACAAAAAAGTCCTATAAAAATAAAGAAAAAAATAACCGCACTTTTAAAGCTGCTGATACAC  >  1/1806830‑1807300

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