Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,482,132	G→A	E73E (GAG→GAA)	unknown →	Type II/IV secretion system ATPase TadZ/CpaE, associated with Flp pilus assembly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,482,132	0	G	A	100.0%	182.4 / NA	68	E73E (GAG→GAA)	unknown	Type II/IV secretion system ATPase TadZ/CpaE, associated with Flp pilus assembly
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/31);  total (37/31)

GTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCTCGCACTATTACTATAGTGTCAAGCCGTGATGATATCCAAGGTGAAGTTGCACAAATTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTCACCTCATCTGATGAAATTTCATTTTCAGCAGAAGACACTGTAGGGGTTATTATTGATATCACACATGAGACTAACATAAAGACTATTGTCGAACGAGTTTTTAGTATAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCCATCTCTTTATCACAAAAACTACTAGATGAGGGCATCCTTTACTTTAATTCGCATACTCAGTTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGCTGATATCCCTCGCTTAAGAGATACTGTTAAAATTGCTGTATTAAGTTGTAAAGGTGGTATTGGCG  >  1/1481892‑1482376                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     gTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATaa                                                                                                                                                                                                                                            >  1:625355/1‑251 (MQ=255) gTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATaa                                                                                                                                                                                                                                            >  1:324068/1‑251 (MQ=255) gTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATaa                                                                                                                                                                                                                                            >  1:408840/1‑251 (MQ=255)                            attaCTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCTTTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTtagtg                                                                                                                                                                                                                 >  1:228663/1‑250 (MQ=255)                             ttaCTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTtagtgt                                                                                                                                                                                                                >  1:282716/1‑249 (MQ=255)                                     tCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACcac                                                                                                                                                                                                        <  1:78445/251‑1 (MQ=255)                                                             tgcacgtaCTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGgtgttgtgt                                                                                                                                                                                     <  1:780066/239‑1 (MQ=255)                                                              gcacgtaCTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACaa                                                                                                                                                                                                      <  1:278232/222‑1 (MQ=255)                                                                   taCTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGTGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACt                                                                                                                                                                                                                                       <  1:782614/189‑1 (MQ=255)                                                                   taCTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGAc                                                                                                                                                                                                                                        >  1:876917/2‑189 (MQ=255)                                                                   taCTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGa                                                                                                                                                                                                                           >  1:691729/2‑202 (MQ=255)                                                                   taCTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGAt                                                                                                                                                                           >  1:560341/2‑250 (MQ=255)                                                                   taCTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATa                                                                                                                                                                          >  1:791165/2‑251 (MQ=255)                                                                     cTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTtag                                                                                                                                                                                                                   >  1:220775/1‑208 (MQ=255)                                                                     cTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTAt                                                                                                                                                                                  >  1:728572/1‑241 (MQ=255)                                                                     cTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGAt                                                                                                                                                                           <  1:820219/248‑1 (MQ=255)                                                                           ctgtaGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGt                                                                                                                                                                                                                                  >  1:511697/4‑187 (MQ=255)                                                                           ctgtaGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTAc                                                                                                                                                                                                           <  1:566742/207‑1 (MQ=255)                                                                                gtatCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAg                                                                                                                                                                         >  1:303041/4‑239 (MQ=255)                                                                                     aaGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACaa                                                                                                                                                                                                      >  1:498813/1‑205 (MQ=255)                                                                                       gCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTaatat                                                                                                                                                                                                                                              >  1:509928/1‑160 (MQ=255)                                                                                              tgatATTCAAGGTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAac                                                                                                                                               <  1:101743/251‑1 (MQ=255)                                                                                                         gTGAAGTTGCACAAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATg                                                                                                                                      <  1:58988/249‑1 (MQ=255)                                                                                                                  cacaAACTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGCGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCctt                                                                                                                           >  1:619879/1‑251 (MQ=255)                                                                                                                      aactCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATg                                                                                                                                      <  1:594277/233‑1 (MQ=255)                                                                                                                       actCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCctttactt                                                                                                                      <  1:42275/249‑1 (MQ=255)                                                                                                                        ctCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACt                                                                                                                                                                    >  1:551800/2‑204 (MQ=255)                                                                                                                          cTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGAt                                                                                                                                                                           <  1:105799/195‑1 (MQ=255)                                                                                                                          cTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTaa                                                                                                                   >  1:547299/1‑251 (MQ=255)                                                                                                                                aCTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCa                                                                                                             >  1:834043/1‑251 (MQ=255)                                                                                                                                         cTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATAcccaatt                                                                                                    >  1:305685/1‑245 (MQ=255)                                                                                                                                             aaaaTATTGAGATTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGt                                                                                                <  1:498609/251‑1 (MQ=255)                                                                                                                                                       gaTTGTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTg                                                                                      >  1:81503/1‑251 (MQ=255)                                                                                                                                                           gTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATgg                                                                                                                                                                                             >  1:406740/1‑144 (MQ=255)                                                                                                                                                           gTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCttt                                                                                                                                                          >  1:534043/1‑179 (MQ=255)                                                                                                                                                           gTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCtt                                                                                                                                                           >  1:133029/1‑178 (MQ=255)                                                                                                                                                           gTTAAAAAGGATTTTTTTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTgaga                                                                                   <  1:331546/250‑1 (MQ=255)                                                                                                                                                                        ttttTACCTCGTCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATcac                                                                                                                                                     >  1:711300/1‑171 (MQ=255)                                                                                                                                                                                 cgtCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGatata                                                            >  1:294837/3‑250 (MQ=255)                                                                                                                                                                                   tCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGatatac                                                           <  1:642527/250‑3 (MQ=255)                                                                                                                                                                                   tCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGatatac                                                           <  1:894702/250‑3 (MQ=255)                                                                                                                                                                                   tCTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGatatac                                                           <  1:223233/250‑3 (MQ=255)                                                                                                                                                                                    cTGATGAGATTTCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGAGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGATATACCt                                                         <  1:232166/251‑1 (MQ=255)                                                                                                                                                                                              ttCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGATATACCTCGTTTAagag                                               <  1:809055/251‑1 (MQ=255)                                                                                                                                                                                              ttCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGATATACCTCGTTTAagag                                               <  1:497566/251‑1 (MQ=255)                                                                                                                                                                                              ttCATTTTCAGCAGAAGATACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGATATACCTCGTTTAagag                                               <  1:368818/251‑1 (MQ=255)                                                                                                                                                                                                           gaagaTACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACaa                                                                                                                                                                                                      >  1:811510/1‑87 (MQ=255)                                                                                                                                                                                                           gaagaTACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGATATACCTCGTTTAAGAGATACTGTTAAAAt                                  >  1:293829/1‑251 (MQ=255)                                                                                                                                                                                                           gaagaTACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGACGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGAGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGATATACCTCGTTTAAGAGATACTGTTAAAAt                                  >  1:272846/1‑251 (MQ=255)                                                                                                                                                                                                            aagaTACTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGATATACCTCGTTTAAGAGATACTGTTaaaa                                   <  1:494435/249‑1 (MQ=255)                                                                                                                                                                                                             agatacTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGAt                                                                                                                                       >  1:370620/5‑148 (MQ=255)                                                                                                                                                                                                             agatacTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGATATACCTCGTTTAAGAGATACTGTTAAAAt                                  <  1:569257/245‑1 (MQ=255)                                                                                                                                                                                                               atacTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGATATACCTCGTTTAag                                                 >  1:484173/3‑232 (MQ=255)                                                                                                                                                                                                                 acTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGatatac                                                           >  1:654308/1‑218 (MQ=255)                                                                                                                                                                                                                  cTGTAGGGATTATTATTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCAt                                                                            <  1:648091/202‑1 (MQ=255)                                                                                                                                                                                                                       gggattattaTTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCctttactt                                                                                                                      >  1:310216/5‑155 (MQ=255)                                                                                                                                                                                                                         gattattaTTGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGa                                                                                                                                        <  1:327238/133‑1 (MQ=255)                                                                                                                                                                                                                             attattGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTATCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGATATACCTCGTTTAAGAGATACTGTTAAAAt                                  <  1:456535/233‑1 (MQ=255)                                                                                                                                                                                                                             attattGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATgg                                                                                                                                                                                             <  1:51721/78‑1 (MQ=255)                                                                                                                                                                                                                                attGATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATcc                                                                                                                             <  1:659086/139‑1 (MQ=255)                                                                                                                                                                                                                                   gATATCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGATATACCTCGTTTAagag                                               >  1:586589/1‑214 (MQ=255)                                                                                                                                                                                                                                    atatCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATAcccaatta                                                                                                   <  1:271883/161‑8 (MQ=255)                                                                                                                                                                                                                                    atatCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATAc                                                                                                          <  1:878003/154‑1 (MQ=255)                                                                                                                                                                                                                                    atatCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGATATACCTCGTTTAagag                                               >  1:103839/1‑213 (MQ=255)                                                                                                                                                                                                                                     tatCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCAtttt                                                                         >  1:295976/1‑186 (MQ=255)                                                                                                                                                                                                                                       tCACACATGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCAC‑AAAACTACTAGATGAGGGTATCctttactt                                                                                                                      <  1:773383/138‑1 (MQ=255)                                                                                                                                                                                                                                            caTGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGATATACCt                                                         <  1:807484/195‑1 (MQ=255)                                                                                                                                                                                                                                            caTGAAACTAATATAAAGACTATTGTCGAACGAGTTTTTAGTGTAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCAATCTCTTTATCACAAAAACTACTAGATGAGGGTATCCTTTACTTTAATTCGCATACCCAATTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGTTGATATACCTCGTTTAAGAGATACTGTTAAAATTGCTGTATTAAGTTGTAAGGGCGGTATTGgcg  >  1:307409/1‑250 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     GTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCTCGCACTATTACTATAGTGTCAAGCCGTGATGATATCCAAGGTGAAGTTGCACAAATTCTACGTACTCGCGGTCTTGAAAATATTGAGATTGTTAAAAAGGATTTTTTCACCTCATCTGATGAAATTTCATTTTCAGCAGAAGACACTGTAGGGGTTATTATTGATATCACACATGAGACTAACATAAAGACTATTGTCGAACGAGTTTTTAGTATAGTACCACAAAATGTATGGTGTTGTGTTATCGGTGATAGTGACTCCATCTCTTTATCACAAAAACTACTAGATGAGGGCATCCTTTACTTTAATTCGCATACTCAGTTAAGTCAGATGGTTGAGAAAATCATTTTAGGTGCTGATATCCCTCGCTTAAGAGATACTGTTAAAATTGCTGTATTAAGTTGTAAAGGTGGTATTGGCG  >  1/1481892‑1482376

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: ‑