Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,480,282	T→G	R89R (CGT→CGG)	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,282	0	T	G	100.0%	188.3 / NA	68	R89R (CGT→CGG)	unknown	Type II/IV secretion system secretin RcpA/CpaC, associated with Flp pilus assembly
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (29/39);  total (29/39)

TGCTGTGCAGTTCTCGGGGCAATGTTCTCATTAAATGCGTACGCACAAAATTTCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGTCCCTTAACATCCGATACTGTTAATGTTGATTCTGTTGTGACCAGTATTGCAGACACGAATAAACAACTTAAAAGCCGTTTCCCAAATACTAATCTATCTGTAAAAAAAGTAGGCAAGGCTTATGTAATTGAGGGAAAAGCCCGTTCTCAAGAAGAAAGCGATGAGGTTCGCAGAATTGTAGGAGAAGCTTTAGGAAGTGGTAGAAAGGTTACTGAA  >  1/1480058‑1480498                                                                                                                                                                                                                                 |                                                                                                                                                                                                                         tGCTGTGCAGTTCTCGGAGCAATGTTCTCATTAAATGCGTACGCACAAAATTTCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGt                                                                                                                                                                                                <  1:300522/251‑1 (MQ=255)  gCTGTGCAGTTCTCGGAGCAATGTTCTCATTAAATGCGTACGCACAAAATTTCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGtt                                                                                                                                                                                               >  1:98374/1‑251 (MQ=255)     gtgCAGTTCTCGGAGCAATGTTCTCATTAAATGCGTACGCACAAAATTTCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGAt                                                                                                                                                                                            >  1:12346/1‑251 (MQ=255)       gCAGTTCTCGGAGCAATGTTCTCATTAAATGCGTACGCACAAAATTTCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGAtt                                                                                                                                                                                           >  1:551242/1‑250 (MQ=255)                 gagCAATGTTCTCATTAAATGCGTACGCACAAAATTTCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCtgttgtcac                                                                                                                                                                                 >  1:398083/3‑247 (MQ=255)                      aTGTTCTCATTAAATGCGTACGCACAAAATTTCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAgcat                                                                                                                                                                           >  1:861191/1‑248 (MQ=255)                           ctcATTAAATGCGTACGCACAAAATTTCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAg                                                                                                                                                                      <  1:297018/251‑1 (MQ=255)                              aTTAAATGCGTACGCACAAAATTTCTCTCTGGATAAAGGCGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGata                                                                                                                                                                   <  1:232526/251‑3 (MQ=255)                                   aTGCGTACGCACAAAATTTCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTc                                                                                                                                                                                          <  1:261636/223‑1 (MQ=255)                                        tACGCACAAAATTTCTCTCTGGATAAAGGTGCAGCCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAaac                                                                                                                                                          <  1:669579/250‑1 (MQ=255)                                        tACGCACAAAATTTCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAaac                                                                                                                                                          <  1:331450/250‑1 (MQ=255)                                           gCACAAAATTTCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAaacaa                                                                                                                                                        >  1:585934/1‑249 (MQ=255)                                               aaaaTTTCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACtgaag                                                                                                                                                  <  1:312855/251‑5 (MQ=255)                                                 aaTTTCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGACTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACtgaagag                                                                                                                                                <  1:809690/251‑7 (MQ=255)                                                    ttCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGcca                                                                                                                                             <  1:442610/251‑2 (MQ=255)                                                       tctctGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCAttt                                                                                                                                          >  1:91018/1‑251 (MQ=255)                                                          ctGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCAttttc                                                                                                                                        >  1:252177/1‑248 (MQ=255)                                                                 aaGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACt                                                                                                                                <  1:809596/251‑1 (MQ=255)                                                                    gTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAat                                                                                                                             >  1:726530/1‑251 (MQ=255)                                                                    gTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAat                                                                                                                             >  1:419808/1‑251 (MQ=255)                                                                    gTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAat                                                                                                                             >  1:379363/1‑251 (MQ=255)                                                                     tGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAatc                                                                                                                            <  1:237979/251‑1 (MQ=255)                                                                     tGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAatc                                                                                                                            <  1:235684/251‑1 (MQ=255)                                                                          cccAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAatctatct                                                                                                                       <  1:398006/251‑1 (MQ=255)                                                                            cAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAatc                                                                                                                            <  1:748234/244‑1 (MQ=255)                                                                                  gTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCGCCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGAt                                                                                                                                                                                            <  1:899541/174‑1 (MQ=255)                                                                                  gTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGtgaagaa                                                                                                               >  1:587261/1‑245 (MQ=255)                                                                                     caaacaaaAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGt                                                                                                            <  1:230276/251‑1 (MQ=255)                                                                                          aaaagaaaagATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAggta                                                                                                       <  1:391479/251‑3 (MQ=255)                                                                                           aaagaaaagaTCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATAc                                                                                                                                 <  1:790791/224‑1 (MQ=255)                                                                                               aaaagaTCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGcc                                                                                                  <  1:226588/251‑2 (MQ=255)                                                                                               aaaagaTCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGc                                                                                                   <  1:321199/250‑1 (MQ=255)                                                                                                    aTCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGt                                                                                             <  1:368532/251‑1 (MQ=255)                                                                                                    aTCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGt                                                                                             <  1:641929/251‑1 (MQ=255)                                                                                                    aTCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGt                                                                                             <  1:580437/251‑1 (MQ=255)                                                                                                    aTCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGt                                                                                             <  1:737393/251‑1 (MQ=255)                                                                                                         cacTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTg                                                                                        <  1:639839/251‑1 (MQ=255)                                                                                                         cacTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTg                                                                                        <  1:554590/251‑1 (MQ=255)                                                                                                         cacTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTg                                                                                        <  1:873240/251‑1 (MQ=255)                                                                                                         cacTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTg                                                                                        <  1:682816/251‑1 (MQ=255)                                                                                                             aTTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAatc                                                                                                                            >  1:301991/1‑211 (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)                                                                                                                                          cTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACtgaag                                                                                                                                                  <  1:551428/160‑5 (MQ=255)                                                                                                                                                     tAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGa                                            <  1:144644/251‑1 (MQ=255)                                                                                                                                                               cacCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTc                                                 <  1:244589/236‑1 (MQ=255)                                                                                                                                                               cacCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGgag                                  >  1:605228/1‑251 (MQ=255)                                                                                                                                                               cacCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGgag                                  >  1:597552/1‑251 (MQ=255)                                                                                                                                                                 ccTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGAAGGTCGGCCTTTGACGCCTGATACCGTTAATGTTGATTCTGTTGTCATCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAatc                                                                                                                            >  1:643571/1‑159 (MQ=255)                                                                                                                                                                  cTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAatc                                                                                                                            >  1:273481/1‑158 (MQ=255)                                                                                                                                                                   tttATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTaa                                                                                           >  1:337503/1‑190 (MQ=255)                                                                                                                                                                   tttATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAATCTATCTGTGAAGAAAGTAGGTAAGGCCTATGTAATTGAAGGTAAAGCTCGTTCACAAGAAGAAAGCGATGAAGTTCGTAGAATTGTAGGAGAAg                               >  1:703990/1‑250 (MQ=255)                                                                                                                                                                                cAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATACTAatc                                                                                                                            >  1:182938/1‑144 (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)                                                                                                                                                                                               gAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGGCCTTTGACGTCTGATACCGTTAATGTTGATTCTGTTGTCACCAGCATTGCAGATACGAATAAACAACTGAAGAGCCATTTTCCAAATAc                                                                                                                                 >  1:322637/1‑124 (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)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                         TGCTGTGCAGTTCTCGGGGCAATGTTCTCATTAAATGCGTACGCACAAAATTTCTCTCTGGATAAAGGTGCAACCCAGTTGGTTCAAACAAAAGAAAAGATCGACACTATTTTCGTCTCTTCACCAAATATTGCCGACTATGAAATTTTAGACGACAACACCTTTATTATCTACGCAAAGGAGGAAGGTAGAACAGAAGTTACTGCTTTCGGAGCAGATGGTCGTCCCTTAACATCCGATACTGTTAATGTTGATTCTGTTGTGACCAGTATTGCAGACACGAATAAACAACTTAAAAGCCGTTTCCCAAATACTAATCTATCTGTAAAAAAAGTAGGCAAGGCTTATGTAATTGAGGGAAAAGCCCGTTCTCAAGAAGAAAGCGATGAGGTTCGCAGAATTGTAGGAGAAGCTTTAGGAAGTGGTAGAAAGGTTACTGAA  >  1/1480058‑1480498

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