Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,480,396	C→T	G127G (GGC→GGT)	unknown →	Type II/IV secretion system secretin RcpA/CpaC, associated with Flp pilus assembly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,480,396	0	C	T	100.0%	116.5 / NA	43	G127G (GGC→GGT)	unknown	Type II/IV secretion system secretin RcpA/CpaC, associated with Flp pilus assembly
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/24);  total (19/24)

AAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGTCCCTTAACATCCGATACTGTTAATGTTGATTCTGTTGTGACCAGTATTGCAGACACGAATAAACAACTTAAAAGCCGTTTCCCAAATACTAATCTATCTGTAAAAAAAGTAGGCAAGGCTTATGTAATTGAGGGAAAAGCCCGTTCTCAAGAAGAAAGCGATGAGGTTCGCAGAATTGTAGGAGAAGCTTTAGGAAGTGGTAGAAAGGTTACTGAAACAAAATTAGGAGAAGACTCTTTACCATTCTTAGACAAATATCATTACGATGGTGTCGTTGATAACGCAAATATTGCCGATACAACCCAAATTAATGTGAAGCTTTCTGTTGTAGAGGTAAATAAAAAACTTTCTGAGGCAAT  >  1/1480152‑1480641                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      aaaagaTCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGcc                                                                                                                                                                                                                                                 <  1:226588/251‑2 (MQ=255) aaaagaTCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGc                                                                                                                                                                                                                                                  <  1:321199/250‑1 (MQ=255)      aTCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGt                                                                                                                                                                                                                                            <  1:737393/251‑1 (MQ=255)      aTCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGt                                                                                                                                                                                                                                            <  1:641929/251‑1 (MQ=255)      aTCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGt                                                                                                                                                                                                                                            <  1:580437/251‑1 (MQ=255)      aTCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGt                                                                                                                                                                                                                                            <  1:368532/251‑1 (MQ=255)           cacTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTg                                                                                                                                                                                                                                       <  1:873240/251‑1 (MQ=255)           cacTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTg                                                                                                                                                                                                                                       <  1:682816/251‑1 (MQ=255)           cacTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTg                                                                                                                                                                                                                                       <  1:639839/251‑1 (MQ=255)           cacTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTg                                                                                                                                                                                                                                       <  1:554590/251‑1 (MQ=255)               aTTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGt                                                                                                                                                                                                                                            <  1:673619/242‑1 (MQ=255)                            cACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGcc                                                                                                                                                                                                                                                 >  1:259137/1‑223 (MQ=255)                            cACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTtca                                                                                                                                                                                                                      >  1:75060/1‑250 (MQ=255)                                aaaTATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGcc                                                                                                                                                                                                                                                 >  1:266301/1‑219 (MQ=255)                                  atatTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGcc                                                                                                                                                                                                                                                 <  1:296150/218‑2 (MQ=255)                                                       tAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGa                                                                                                                                                                                           <  1:144644/251‑1 (MQ=255)                                                                 cacCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTc                                                                                                                                                                                                <  1:244589/236‑1 (MQ=255)                                                                 cacCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGgag                                                                                                                                                                                 >  1:597552/1‑251 (MQ=255)                                                                 cacCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGgag                                                                                                                                                                                 >  1:605228/1‑251 (MQ=255)                                                                     tttATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTaa                                                                                                                                                                                                                                          >  1:337503/1‑190 (MQ=255)                                                                     tttATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAg                                                                                                                                                                              >  1:703990/1‑250 (MQ=255)                                                                                  cAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCt                                                                                                                                                                            >  1:676765/1‑239 (MQ=255)                                                                                       gaggaAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGaaggt                                                                                                                                                                                                                                  <  1:224996/180‑5 (MQ=255)                                                                                              gTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTg                                                                                                                                                                                                                                       >  1:519009/1‑168 (MQ=255)                                                                                              gTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGAAACGGTCGAAAGGTTa                                                                                                                                                      >  1:726356/1‑249 (MQ=255)                                                                                                aGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGcc                                                                                                                                                                                                                                                 <  1:658551/156‑2 (MQ=255)                                                                                                 gAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGc                                                                                                                                                                                                                                                  >  1:112598/1‑154 (MQ=255)                                                                                                 gAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGAAACGGTCGAAAGGTTACTGaa                                                                                                                                                 >  1:818469/1‑251 (MQ=255)                                                                                                  aaCAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGAAACGGTCGAAAGGTTACTGa                                                                                                                                                  <  1:80305/249‑1 (MQ=255)                                                                                                   aCAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGa                                                                                                                                                                      <  1:395681/228‑1 (MQ=255)                                                                                                    cAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTg                                                                                                                                                                                                                                       <  1:455650/162‑1 (MQ=255)                                                                                                      gAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGa                                                                                                                                                                      <  1:556071/225‑1 (MQ=255)                                                                                                                                                                                                                  tttccAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGAAACGGTCGAAAGGTTACTGAAACAAAGTTAGGGGAAGATGCTTTGCCGTTCTTAGACAAATACCATTACGATGGTGTTGTTGATAACGCAAATATTGCAGACACAACCCAAATTAATGTGAAGCTTTCtgttgt                                >  1:567193/4‑251 (MQ=255)                                                                                                                                                                                                                    tccAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGAAACGGTCGAAAGGTTACTGAAACAAAGTTAGGGGAAGATGCTTTGCCGTTCTTAGACAAATACCATTACGATGGTGTTGTTGATAACGCAAATATTGCAGACACAACCCAAATTAATGTGAAGCTTTCTGTTGTag                              <  1:591569/250‑1 (MQ=255)                                                                                                                                                                                                                    tccAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGAAACGGTCGAAAGGTTACTGAAACAAAGTTAGGGGAAGATGCTTTGCCGTTCTTAGACAAATACCATTACGATGGTGTTGTTGATAACGCAAATATTGCAGACACAACCCAAATTAATGTGAAGCTTTCTGTTGTag                              <  1:735964/250‑1 (MQ=255)                                                                                                                                                                                                                    tccAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGAAACGGTCGAAAGGTTACTGAAACAAAGTTAGGGGAAGATGCTTTGCCGTTCTTAGACAAATACCATTACGATGGTGTTGTTGATAACGCAAATATTGCAGACACAACCCAAATTAATGTGAAGCTTTCTGTTGTag                              <  1:742506/250‑1 (MQ=255)                                                                                                                                                                                                                      cAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGAAACGGTCGAAAGGTTACTGAAACAAAGTTAGGGGAAGATGCTTTGCCGTTCTTAGACAAATACCATTACGATGGTGTTGTTGATAACGCAAATATTGCAGACACAACCCAAATTAATGTGAAGCTTTCTGTTGTagag                            <  1:251446/251‑1 (MQ=255)                                                                                                                                                                                                                         aTACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGAAACGGTCGAAAGGTTACTGAAACAAAGTTAGGGGAAGATGCTTTGCCGTTCTTAGACAAATACCATTACGATGGTGTTGTTGATAACGCAAATATTGCAGACACAACCCAAATTAATGTGAAGCTTTCTGTTGTAGAGGta                         >  1:795982/1‑251 (MQ=255)                                                                                                                                                                                                                         aTACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGAAACGGTCGAAAGGTTACTGAAACAAAGTTAGGGGAAGATGCTTTGCCGTTCTTAGACAAATACCATTACGATGGTGTTGTTGATAACGCAAATATTGCAGACACAACCCAAATTAATGTGAAGCTTTCTGTTGTAGAGGta                         >  1:862888/1‑251 (MQ=255)                                                                                                                                                                                                                            cTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGAAACGGTCGAAAGGTTACTGAAACAAAGTTAGGGGAAGATGCTTTGCCGTTCTTAGACAAATACCATTACGATGGTGTTGTTGATAACGCAAATATTGCAGACACAACCCAAATTAATGTGAAGCTTTCTGTTGTAGAGGtaaaa                      >  1:837636/1‑250 (MQ=255)                                                                                                                                                                                                                                   atctGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGAAACGGTCGAAAGGTTACTGAAACAAAGTTAGGGGAAGATGCTTTGCCGTTCTTAGACAAATACCATTACGATGGTGTTGTTGATAACGCAAATATTGCAGACACAACCCAAATTAATGTGAAGCTTTCTGTTGTAGAGGTAAATaagaaa                >  1:176042/1‑250 (MQ=255)                                                                                                                                                                                                                                       gtgaagaaaGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGAAACGGTCGAAAGGTTACTGAAACAAAGTTAGGGGAAGATGCTTTGCCGTTCTTAGACAAATACCATTACGATGGTGTTGTTGATAACGCAAATATTGCAGACACAACCCAAATTAATGTGAAGCTTTCTGTTGTAGAGGTAAATAAGAAACttt            >  1:216054/7‑250 (MQ=255)                                                                                                                                                                                                                                                gTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAGCTTTAGGAAACGGTCGAAAGGTTACTGAAACAAAGTTAGGGGAAGATGCTTTGCCGTTCTTAGACAAATACCATTACGATGGTGTTGTTGATAACGCAAATATTGCAGACACAACCCAAATTAATGTGAAGCTTTCTGTTGTAGAGGTAAATAAGAAACTTTCTGAGGcgat  >  1:244491/1‑248 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      AAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGTCCCTTAACATCCGATACTGTTAATGTTGATTCTGTTGTGACCAGTATTGCAGACACGAATAAACAACTTAAAAGCCGTTTCCCAAATACTAATCTATCTGTAAAAAAAGTAGGCAAGGCTTATGTAATTGAGGGAAAAGCCCGTTCTCAAGAAGAAAGCGATGAGGTTCGCAGAATTGTAGGAGAAGCTTTAGGAAGTGGTAGAAAGGTTACTGAAACAAAATTAGGAGAAGACTCTTTACCATTCTTAGACAAATATCATTACGATGGTGTCGTTGATAACGCAAATATTGCCGATACAACCCAAATTAATGTGAAGCTTTCTGTTGTAGAGGTAAATAAAAAACTTTCTGAGGCAAT  >  1/1480152‑1480641

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