Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	203,596	A→C	E295A (GAA→GCA)	unknown →	Multidrug resistance transporter => Bicyclomycin resistance protein Bcr‑1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	203,596	0	A	C	100.0%	174.6 / NA	63	E295A (GAA→GCA)	unknown	Multidrug resistance transporter => Bicyclomycin resistance protein Bcr‑1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (35/28);  total (35/28)

GTGCTGGGGTATATGTTCGCCGCTTCTTTTGCTTTCGGGGGCTTATTTTCTTTCATCACGGCGGGGTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGAAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAATTGCAATTGGCATTGCGATTTTTGTGGGACAAAACTCGCTGATTTCTTCCAACGCCATGGCGTCCATTTTAGAAAAATTCCCGAATGCCGCCTGCTCGGCAAATTCTATGGTCGGTAGTGTACGTTTCGGCACCGGCGCCTGCGTGGGACACTGGTTGCCCTTATGAATATGAGCAGCGCAACGCCA  >  1/203361‑203830                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           gTGCTGGGGTATATGTTCGCCGCTTCTTTTGCTTTCGGGGGCTTATTTTCTTTCATCACGGCGGGGTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTc                                                                                                                                                                                                                              <  1:479680/250‑1 (MQ=255)   gCTGGGGTATATGTTCGCCGCTTCTTTTGCTTTCGGGGGCTTATTTTCTTTCATCACGGCGGGGTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTgg                                                                                                                                                                                                                           >  1:344846/1‑251 (MQ=255)   gCTGGGGTATATGTTCGCCGCTTCTTTTGCTTTCGGGGGCTTATTTTCTTTCATCACGGCGGGGTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGGTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTgg                                                                                                                                                                                                                           >  1:265878/1‑251 (MQ=255)      ggggTATATGTTCGCCGCTTCTTTTGCTTTCGGGGGCTTATTTTCTTTCATCACGGCGGGGTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTg                                                                                                                                                                                                                         >  1:467261/1‑250 (MQ=255)       gggTATATGTTCGCCGCTTCTTTTGCTTTCGGGGGCTTATTTTCTTTCATCACGGCGGGGTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTg                                                                                                                                                                                                                         >  1:519211/1‑249 (MQ=255)       gggTATATGTTCGCCGCTTCTTTTGCTTTCGGGGGCTTATTTTCTTTCATCACGGCGGGGTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTg                                                                                                                                                                                                                         >  1:572976/1‑249 (MQ=255)        ggTATATGTTCGCCGCTTCTTTTGCTTTCGGGGGCTTATTTTCTTTCATCACGGCGGGGTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTg                                                                                                                                                                                                                         >  1:112726/1‑248 (MQ=255)         gTATATGTTCGCCGCTTCTTTTGCTTTCGGGGGCTTATTTTCTTTCATCACGGCGGGGTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATTGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTg                                                                                                                                                                                                                         >  1:594132/1‑247 (MQ=255)                          cttTTGCTTTCGGGGGCTTATTTTCTTTCATCACGGCGGGGTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTTGGCttt                                                                                                                                                                                                    >  1:399684/1‑251 (MQ=255)                                          cTTATTTTCTTTCATCACGGCGGGGTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacac                                                                                                                                                                                    <  1:535167‑M1/251‑11 (MQ=255)                                                               ggggTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaag                                                                                                                                                                 <  1:856183‑M1/249‑30 (MQ=255)                                                                 ggTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGc                                                                                                                                                                                                       <  1:218690/209‑1 (MQ=255)                                                                 ggTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaat                                                                                                                                                             >  1:697935‑M1/1‑218 (MQ=255)                                                                 ggTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaat                                                                                                                                                             >  1:561024‑M1/1‑218 (MQ=255)                                                                    cTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTg                                                                                                                                                                                                                         <  1:718642/188‑1 (MQ=255)                                                                      aTTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATttggtt                                                                                                                                                                                                                                   <  1:631549/176‑1 (MQ=255)                                                                      aTTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTg                                                                                                                                                                                                                         <  1:166105/186‑1 (MQ=255)                                                                      aTTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGaa                                                                                                                                                                                                              <  1:187513/197‑1 (MQ=255)                                                                      aTTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaa                                                                                                                                                          >  1:634167‑M1/1‑213 (MQ=255)                                                                      aTTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacac                                                                                                                                                                                    >  1:469450‑M1/1‑213 (MQ=255)                                                                      aTTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATAGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTCTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacac                                                                                                                                                                                    >  1:282712‑M1/1‑213 (MQ=255)                                                                         gtgtATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGaa                                                                                                                                                                                                              <  1:889312/194‑1 (MQ=255)                                                                         gtgtATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaa                                                                                                                                                              <  1:79665‑M1/242‑33 (MQ=255)                                                                          tgtATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactaggg                                                                                                                                                    <  1:274667‑M1/251‑43 (MQ=255)                                                                           gtATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactaggg                                                                                                                                                    <  1:170403‑M1/250‑43 (MQ=255)                                                                             atatCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaa                                                                                                                                                                                         >  1:305983‑M1/1‑206 (MQ=255)                                                                             atatCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTGCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggaga                                                                                                                                                 <  1:453982‑M1/251‑46 (MQ=255)                                                                                 ctgcaTTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgt                                                                                                                                             <  1:789438‑M1/249‑50 (MQ=255)                                                                                                 aaGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaat                                                                                                                             <  1:488163‑M1/251‑66 (MQ=255)                                                                                                    ccGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataag                                                                                                                          <  1:61844‑M1/251‑69 (MQ=255)                                                                                                               cGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaac                                                                                                               <  1:652583‑M1/251‑80 (MQ=255)                                                                                                                ggTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaac                                                                                                               >  1:624182‑M1/1‑171 (MQ=255)                                                                                                                       tCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaa                                                                                                       >  1:897654‑M1/1‑164 (MQ=255)                                                                                                                       tCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaa                                                                                                       >  1:220939‑M1/1‑164 (MQ=255)                                                                                                                       tCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaa                                                                                                                  <  1:690272‑M1/240‑77 (MQ=255)                                                                                                                            aTGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaag                                                                                                  >  1:240038‑M1/1‑159 (MQ=255)                                                                                                                               gatgtgtataaggggcaggcctttggttccctcctgaatgGGCGGTTGGGGCATAAAGGCGGCGCGGAAACTATGTTGGGGGGTGGGTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaaggg                                                                                                <  1:196400‑M1/231‑95 (MQ=255)                                                                                                                                     atCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaa                                                                                                                                                                   >  1:20038‑M1/1‑150 (MQ=255)                                                                                                                                      tctcGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaa                                                                                        >  1:677987‑M1/1‑149 (MQ=255)                                                                                                                                      tctcGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataag                                                                                                                          >  1:741566‑M1/1‑149 (MQ=255)                                                                                                                                           gTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTggg                                                                                                                                                                                                >  1:368420/1‑142 (MQ=255)                                                                                                                                            tgatgaCCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaac                                                                                   <  1:725008‑M1/250‑108 (MQ=255)                                                                                                                                                 aCCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgt                                                                              >  1:167400‑M1/1‑138 (MQ=255)                                                                                                                                                  ccTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgtc                                                                             >  1:832489‑M1/1‑137 (MQ=255)                                                                                                                                                  ccTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgtc                                                                             >  1:212152‑M1/1‑137 (MQ=255)                                                                                                                                                  ccTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagt                                                                                                                         <  1:484441‑M1/206‑70 (MQ=255)                                                                                                                                                  ccTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaag                                                                                                                                    <  1:875583‑M1/195‑59 (MQ=255)                                                                                                                                                  ccTTTGGTTCCTTCCTGAATGGGCGATTGGGGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgtc                                                                             >  1:335550‑M1/1‑137 (MQ=255)                                                                                                                                                   cTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgtc                                                                             >  1:147901‑M1/1‑136 (MQ=255)                                                                                                                                                             ttcctGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgtccttttttattta                                                                 >  1:765031‑M1/1‑126 (MQ=255)                                                                                                                                                                       ggCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgtccttttttatttaagaatagac                                                        <  1:441102‑M1/250‑135 (MQ=255)                                                                                                                                                                               gTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgtccttttttatttaagaatagacaagaataaa                                               >  1:601919‑M1/1‑108 (MQ=255)                                                                                                                                                                               gTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaat                                                                                                                                        >  1:735345‑M1/1‑108 (MQ=255)                                                                                                                                                                                 gCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactaggg                                                                                                                                                    >  1:285213‑M1/1‑106 (MQ=255)                                                                                                                                                                                  cATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgtccttttttatttaagaatagacaagaataaaaaa                                            <  1:677010‑M1/251‑147 (MQ=255)                                                                                                                                                                                  cATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgtccttttttatttaagaatagacaagaataaaaaa                                            <  1:320532‑M1/251‑147 (MQ=255)                                                                                                                                                                                  cATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgtccttttttatttaagaatagacaagaataaaaaa                                            <  1:19564‑M1/251‑147 (MQ=255)                                                                                                                                                                                       ggTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtagaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgtccttttttatttaagaatagacaagaat                                                  >  1:451625‑M1/1‑100 (MQ=255)                                                                                                                                                                                                 gAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaaggg                                                                                                >  1:159427‑M1/1‑90 (MQ=255)                                                                                                                                                                                                   aaCTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacgg                                                                                 >  1:297207‑M1/1‑88 (MQ=255)                                                                                                                                                                                                    aCTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTCGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgtccttttttatttaagaatagacaagaataaaaaatgcgaaaagttttaaaaa                          <  1:444571‑M1/251‑165 (MQ=255)                                                                                                                                                                                                                         aTGGTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaataaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgtccttttttatttaagaatagacaagaataaaaaatgcgaaaagttttaaaaatcgcaaaaaaaaagcgaaac      <  1:446769‑M1/250‑185 (MQ=255)                                                                                                                                                                                                                            gTTCAGTTTTTGGCGGCAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAtgtaaaacaccgaacaaaatgagaaaaaggaataaactagggagaatgttgaatcaagcgggaataagtggaacaaaaccaatgaaaacaagggcgcaagaaaaaacggcgtccttttttatttaagaatagacaagaataaaaaatgcgaaaagttttaaaaatcgcaaaaaaaaagcgaaacaaaa  >  1:788025‑M1/1‑63 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           GTGCTGGGGTATATGTTCGCCGCTTCTTTTGCTTTCGGGGGCTTATTTTCTTTCATCACGGCGGGGTCTATTGTGTATATCGGCATTTACGGCATTAAGCCGGAGAATTTCGGTTATTTCTTTATGCTGAATATCTCGGTGATGACCTTTGGTTCCTTCCTGAATGGGCGATTGGTGCATAAGGTCGGCGCGGAAACTATGTTGCGCGTTGGTTTGATGGTTCAGTTTTTGGCGGAAATTTGGTTGGTTCTGGTGCTGTTGCTGAATTTGGGCTTTTGGGCAATTGCAATTGGCATTGCGATTTTTGTGGGACAAAACTCGCTGATTTCTTCCAACGCCATGGCGTCCATTTTAGAAAAATTCCCGAATGCCGCCTGCTCGGCAAATTCTATGGTCGGTAGTGTACGTTTCGGCACCGGCGCCTGCGTGGGACACTGGTTGCCCTTATGAATATGAGCAGCGCAACGCCA  >  1/203361‑203830

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