Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,615,823	C→G	T391S (ACT→AGT)	unknown →	Glycine betaine ABC transport system, permease protein OpuAB / Glycine betaine ABC transport system, glycine betaine‑binding protein OpuAC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,615,823	0	C	G	100.0%	214.3 / NA	70	T391S (ACT→AGT)	unknown	Glycine betaine ABC transport system, permease protein OpuAB / Glycine betaine ABC transport system, glycine betaine‑binding protein OpuAC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (31/39);  total (31/39)

CTCCTGAAACCACTTATCAGCTGAGTAAACAATTACTGGCGGAACAGTTCCGGCTTGAA‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTACTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAATGAAAGCCGATTTTGCAGCAAAAAATCCGCAGGTTCTGACCGCACTTAACCGCCTCGCCAATCGTATTTCCGAA  >  1/1615580‑1616067                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     ctcCTGAAACCACTTATCAGCTGAGTAAACAATTACTGGCGGAACAGTTCCGGCTTGAA‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATtagtttagtt                                                                                                                                                                                                                                               >  1:262430/1‑251 (MQ=255)  ccctgaacccccttaccagttgaaaaaaCAATTCCTGGCGAAACAGTTCCGGTTTAAT‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTc                                                                                                                                                                                                                                              <  1:524181/226‑1 (MQ=255)        aaCCATTTATCAGCTGAGAAAACAATTCCTGGCGGAACAGTTCCGGTTTAAT‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTg                                                                                                                                                                                                                                        <  1:490782/251‑1 (MQ=255)             cTTATCAGCTGAGAAAACAATTCCGGGCGAAACAGTTCCGGTTTAAA‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATtagtttagtt                                                                                                                                                                                                                                               <  1:3881/239‑1 (MQ=255)                     ttgaatAAACATTTACTGGCGGAACAGTTCCGGTTTGAT‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAAc                                                                                                                                                                                                                            <  1:547161/245‑1 (MQ=255)                     cTGAGTAAACAATTACTGGCGAAACATTCCCGGCTTGAA‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAAcc                                                                                                                                                                                                                           <  1:460958/251‑1 (MQ=255)                     cTGAGAAAACAATTACTGGCGGAACATTTCCGGTTTGAA‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAAc                                                                                                                                                                                                                            <  1:826530/250‑1 (MQ=255)                       gAGTAAACAATTACTGGCGGAACAGTTCCGGCTTGAA‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAAc                                                                                                                                                                                                                            >  1:46748/1‑248 (MQ=255)                          aaaaaaattattggggaaaCAGTTCCGGC‑TGAATTTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCtt                                                                                                                                                                                                                      <  1:635245/234‑1 (MQ=255)                              aaattagtggggaaaCAGTTCGGGTTTGAT‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCg                                                                                                                                                                                                                  <  1:532941/238‑1 (MQ=255)                               aaTTACGGGCGGAACAGTTCCGGTTTAAT‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGt                                                                                                                                                                                                                 <  1:180647/251‑1 (MQ=255)                                    cgggCGGAACAGTCCCGGCTTGAA‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATAc                                                                                                                                                                                                            <  1:68268/249‑1 (MQ=255)                                    cgggCGAAACAGTTCCGGTTTAAT‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATAc                                                                                                                                                                                                            <  1:175743/249‑1 (MQ=255)                                      ggCGGAACAGTTCCGGCTTGAA‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACgg                                                                                                                                                                                                          >  1:576540/1‑251 (MQ=255)                                          gAACAGTTCCGGTTTGAA‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGt                                                                                                                                                                                                       <  1:831791/250‑1 (MQ=255)                                               gTTCCGGCTTGAA‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTc                                                                                                                                                                                                  >  1:468864/1‑250 (MQ=255)                                                                                atCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGAtt                                                                                                                                                                                  >  1:667260/1‑234 (MQ=255)                                                                                 tCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAAcc                                                                                                                                                                                                                           <  1:901726/192‑1 (MQ=255)                                                                                 tCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTa                                                                                                                                                                >  1:216417/1‑251 (MQ=255)                                                                                  cAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATAtt                                                                                                                                                                                                                                         <  1:27956/177‑1 (MQ=255)                                                                                  cAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATAc                                                                                                                                                                                                            <  1:158850/206‑1 (MQ=255)                                                                                  cAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTAc                                                                                                                                                               <  1:73718/251‑1 (MQ=255)                                                                                   aaaaTACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTa                                                                                                                                                                >  1:307257/1‑249 (MQ=255)                                                                                          cTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAg                                                                                                                                                       >  1:341463/1‑251 (MQ=255)                                                                                            aTGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGttt                                                                                                                                                                                                                                                   <  1:240694/157‑1 (MQ=255)                                                                                            aTGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCAccc                                                                                                                                                                                                                                  <  1:207543/174‑1 (MQ=255)                                                                                            aTGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTc                                                                                                                                                                                                  <  1:799594/206‑1 (MQ=255)                                                                                               cATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATAc                                                                                                                                                                                                            >  1:348862/1‑193 (MQ=255)                                                                                                      aGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTg                                                                                                                                                                                                                               <  1:826319/167‑1 (MQ=255)                                                                                                         gAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATAc                                                                                                                                                                                                            <  1:123424/183‑1 (MQ=255)                                                                                                          aaaGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATc                                                                                                                                        <  1:416916/250‑1 (MQ=255)                                                                                                          aaaGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCa                                                                                                                                       >  1:435523/1‑251 (MQ=255)                                                                                                          aaaGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCa                                                                                                                                       >  1:420524/1‑251 (MQ=255)                                                                                                               agCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAaaca                                                                                                                                             >  1:125474/1‑240 (MQ=255)                                                                                                               agCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATc                                                                                                                                        <  1:690143/245‑1 (MQ=255)                                                                                                               agCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTc                                                                                                                                   >  1:857663/1‑250 (MQ=255)                                                                                                                       gAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATAc                                                                                                                                                                                                            <  1:850838/169‑1 (MQ=255)                                                                                                                            acacATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTc                                                                                                                                                                                                  <  1:470296/174‑1 (MQ=255)                                                                                                                            acacATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGAt                                                                                                                                                                                   >  1:77953/1‑189 (MQ=255)                                                                                                                            acacATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTc                                                                                                                                   <  1:327730/237‑1 (MQ=255)                                                                                                                            acacATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAga                                                                                                                     >  1:583701/1‑251 (MQ=255)                                                                                                                            acacATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAga                                                                                                                     >  1:725171/1‑251 (MQ=255)                                                                                                                             cacaTGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAgac                                                                                                                    <  1:827319/251‑1 (MQ=255)                                                                                                                                   ggCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACt                                                                                                                                                              <  1:506139/203‑1 (MQ=255)                                                                                                                                         caacaaTCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAAttt                                                                                                         <  1:24468/250‑1 (MQ=255)                                                                                                                                           acaaTCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTAtt                                                                                                      >  1:85715/1‑251 (MQ=255)                                                                                                                                            caaTCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTc                                                                                                     <  1:738203/251‑1 (MQ=255)                                                                                                                                                     aTTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACAtt                                                                                                              >  1:231225/1‑233 (MQ=255)                                                                                                                                                     aTTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTAt                                                                                            >  1:547176/1‑251 (MQ=255)                                                                                                                                                        tACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCaa                                                                                         <  1:525730/251‑1 (MQ=255)                                                                                                                                                          catcatGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAgg                                                                                       <  1:818342/251‑1 (MQ=255)                                                                                                                                                          catcatGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAgg                                                                                       <  1:70151/251‑1 (MQ=255)                                                                                                                                                                   gCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATCAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCAc                                                                               <  1:440554/250‑1 (MQ=255)                                                                                                                                                                     aGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTaa                                                                            >  1:736702/1‑251 (MQ=255)                                                                                                                                                                      ggTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAAc                                                                                                                                                                                                                            <  1:876138/106‑1 (MQ=255)                                                                                                                                                                      ggTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAAt                                                                           >  1:163211/1‑251 (MQ=255)                                                                                                                                                                        tGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAATg                                                                          <  1:472017/250‑1 (MQ=255)                                                                                                                                                                                      gCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAATGAAAGCCGATTTTgca                                                           >  1:561356/1‑251 (MQ=255)                                                                                                                                                                                           acacTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTAc                                                                                                                                                               <  1:388453/146‑1 (MQ=255)                                                                                                                                                                                            cacTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAATGAAAGCCGATTTTGCAGCaaaa                                                     >  1:418443/1‑251 (MQ=255)                                                                                                                                                                                                        tcgatcgaGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAATGAAAGCCGATTTTGCAGCAAAAAATCCGCAGGtt                                         <  1:696226/251‑1 (MQ=255)                                                                                                                                                                                                                      ggCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCGCCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACAt                                                                                                                                         <  1:644468/141‑1 (MQ=255)                                                                                                                                                                                                                       gCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAATGAAAGCCGATTTTGCAGCAAAAAATCCGCAGGTTCTGACCGCACTTaa                           >  1:207179/1‑250 (MQ=255)                                                                                                                                                                                                                       gCTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAATGAAAGCCGATTTTGCAGCAAAAAATCCGCAGGTTCTGACCGCACTTaa                           >  1:857517/1‑250 (MQ=255)                                                                                                                                                                                                                        cTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTcaac                                                                                                                                                                                               <  1:684418/85‑1 (MQ=255)                                                                                                                                                                                                                        cTATAAAGGAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAATGAAAGCCGATTTTGCAGCAAAAAATCCGCAGGTTCTGACCGCACTTAAc                          >  1:423928/1‑250 (MQ=255)                                                                                                                                                                                                                                gAATAAAATCACAAGGCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAATGAAAGCCGATTTTGCAGCAAAAAATCCGCAGGTTCTGACCGCACTTAACCGCCTCGc                  >  1:472300/1‑250 (MQ=255)                                                                                                                                                                                                                                              ggCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAATGAAAGCCGATTATGCAGCAAAAAATCCGCAGGTTCTGACCGCACTTAAc                          >  1:868065/1‑228 (MQ=255)                                                                                                                                                                                                                                               gCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAgg                                                                                       >  1:411598/1‑166 (MQ=255)                                                                                                                                                                                                                                               gCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAATGAAAGCCGATTTTGCAGCAAAAAATCCGCAGGTTCTGACCGCACTTAACCGCCTCGCCAATCGTATTTc      >  1:245905/1‑247 (MQ=255)                                                                                                                                                                                                                                               gCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAATGAAAGCCGATTTTGCAGCAAAAAATCCGCAGGTTCTGACCGCACTTAACCGCCTCGCCAATCGTATTTCCGaa  >  1:316814/1‑251 (MQ=255)                                                                                                                                                                                                                                               gCATTAGTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAATGAAAGCCGATTTTGCAGCAAAAAATCCGCAGGTTCTGACCGCACTTAACCGCCTCGCCAATCGTATTTCCGaa  >  1:727400/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     CTCCTGAAACCACTTATCAGCTGAGTAAACAATTACTGGCGGAACAGTTCCGGCTTGAA‑TTTTTACCGCCGATGGCATATCAAAATACCTATGCATTAGCAGTGAAAGAAGCCTATGGAGCAACACATGGGCTTGCAACAATCAGCGATTTACATCATGTGGCAGGTGACATTCGCGCCGGCTTCACACTGGAATTTCTCGATCGAGTCGATGGCTATAAAGGAATAAAATCACAAGGCATTACTTTAGTTCATATTGTCACCCTTGAACCTGCTTTGCGTTATACGGCGTTACTCAACAATAAAATTGATTTAATCGAAGCATTTTCTACTGACGCAGAACTCAAACAACATCAGCTCCGCTTACTGAAAGACGACATTAATTTATTCCCTGCCTATCAAGGTGCGCCACTAATGAAAGCCGATTTTGCAGCAAAAAATCCGCAGGTTCTGACCGCACTTAACCGCCTCGCCAATCGTATTTCCGAA  >  1/1615580‑1616067

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