Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,486,515	A→C	G148G (GGA→GGC)	unknown →	Flp pilus assembly protein TadD, contains TPR repeat

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,486,515	0	A	C	100.0%	187.8 / NA	65	G148G (GGA→GGC)	unknown	Flp pilus assembly protein TadD, contains TPR repeat
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (25/40);  total (25/40)

TTATAAATTATCTGAAATTTATTACAAGAAAGGGGACAGTAATTCATCTCTACTCTATTTACAGCCGTTGCTAACAAACAGTGGGCAGTTAAGTGAAAAGGCCAAAGTTCTACAAGCAAGAAATTTAAATCAACTCAAGCGCTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCCAAAAACGGGGAAGTTTATAATCTTCGTGGTGTGACTTATGCCTTAATGGGAAATCCTAATAAAGCCAACGAAGATATAAATAAAGCCCGCGAGTATTTTTTAAACGATGCTGTAGCGGTAAACAATCTTGCTATGCTAAGTATCGTCAACGGAGATTATCGTAATGCCGTTTCATTGTTATTACCTCAGTATTTAAATGGGGTTAGGGAACA  >  1/1486275‑1486676                                                                                                                                                                                                                                                 |                                                                                                                                                                  ttATAAATTATCTGAAATTTATTACAAGAAAGGGGACAGTAGTTCATCTCTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCaaata                                                                                                                                                         <  1:897787/251‑1 (MQ=255) ttATAAATTATCTGAAATTTATTACAAGAAAGGGGACAGTAGTTCATCTCTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCaaata                                                                                                                                                         <  1:713119/251‑1 (MQ=255)   ataAATTATCTGAAATTTATTACAAGAAAGGGGACAGTAGTTCATCTCTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCaaata                                                                                                                                                         >  1:650422/1‑249 (MQ=255)                     attaCAAGAAAGGGGACAGTAGTTCATCTCTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGatattaac                                                                                                                                     <  1:264888/251‑5 (MQ=255)                                     cAGTAGTTCATCTCTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCaaata                                                                                                                                                         >  1:330978/1‑215 (MQ=255)                                          gttCATCTCTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTa                                                                                                                >  1:627529/2‑251 (MQ=255)                                              aTCTCTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTaaatg                                                                                                            <  1:631311/251‑3 (MQ=255)                                               tctcTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGatattaac                                                                                                                                     <  1:641194/225‑5 (MQ=255)                                               tctcTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGatattaac                                                                                                                                     <  1:622692/225‑5 (MQ=255)                                               tctcTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGcacgagaat                                                                                                                        >  1:781689/1‑230 (MQ=255)                                               tctcTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTaaatgg                                                                                                           >  1:628643/1‑248 (MQ=255)                                               tctcTACTCTACTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTaaatga                                                                                                           >  1:89663/1‑248 (MQ=255)                                                  ctactctactTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGa                                                                                                                                               <  1:363830/212‑1 (MQ=255)                                                  ctactctactTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTaaatg                                                                                                            <  1:590627/247‑3 (MQ=255)                                                  ctactctactTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAatgatg                                                                                                         >  1:671526/1‑250 (MQ=255)                                                    actctaCTTACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGcag                                                                                                      <  1:173912/251‑3 (MQ=255)                                                         acttACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGa                                                                                                                                               >  1:521618/3‑205 (MQ=255)                                                         acttACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGcagtggca                                                                                                 >  1:652437/3‑244 (MQ=255)                                                          cttACAGCCGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGcagtggcag                                                                                                >  1:674347/2‑243 (MQ=255)                                                                 ccGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTaaat                                                                                           <  1:854547/249‑2 (MQ=255)                                                                 ccGTTGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGCAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTaaat                                                                                           <  1:784591/249‑2 (MQ=255)                                                                     tGCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatg                                                                                     >  1:23537/1‑248 (MQ=255)                                                                      gCTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTaaat                                                                                           <  1:352841/244‑2 (MQ=255)                                                                       cTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatggc                                                                                   <  1:364789/251‑6 (MQ=255)                                                                       cTAACAAACAGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatggc                                                                                   <  1:620879/251‑6 (MQ=255)                                                                         aacaaacaGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatggcta                                                                                 <  1:568789/251‑8 (MQ=255)                                                                         aacaaacaGCGGGCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatggc                                                                                   >  1:691556/1‑244 (MQ=255)                                                                                  cgggCAGTTAAGTGAAAGGGCCAAAAGTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATggctatgcta                                                                            <  1:157384/246‑1 (MQ=255)                                                                                    ggCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTAt                                                                       <  1:826444/250‑1 (MQ=255)                                                                                    ggCAGTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTAt                                                                       <  1:625293/250‑1 (MQ=255)                                                                                        gTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCaaataa                                                                                                                                                        <  1:303771/165‑1 (MQ=255)                                                                                        gTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatggct                                                                                  >  1:362741/1‑229 (MQ=255)                                                                                        gTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATcatta                                                                  >  1:535635/1‑247 (MQ=255)                                                                                        gTTAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATcatta                                                                  >  1:704594/1‑247 (MQ=255)                                                                                         ttAAGTGAAAGGGCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTCCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATcatta                                                                  <  1:18251/250‑5 (MQ=255)                                                                                                    ggCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTa                                                                        <  1:522092/233‑1 (MQ=255)                                                                                                    ggCCAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTAt                                                       <  1:620741/250‑1 (MQ=255)                                                                                                      ccAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGt                                                    <  1:117556/251‑1 (MQ=255)                                                                                                      ccAAAATTCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGt                                                    <  1:151741/251‑1 (MQ=255)                                                                                                         aaattCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAAt                                                 >  1:297629/4‑251 (MQ=255)                                                                                                           attCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGc                                               <  1:557515/250‑1 (MQ=255)                                                                                                           attCTACAAGCACGAAATTTAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGc                                               <  1:858571/250‑1 (MQ=255)                                                                                                                             ttAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCattattattgcc                             <  1:589520/251‑4 (MQ=255)                                                                                                                             ttAAATCAACTTAAGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCattattattgcc                             <  1:27851/251‑4 (MQ=255)                                                                                                                                     acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCaaat                                                                                                                                                          >  1:524525/5‑118 (MQ=255)                                                                                                                                     acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatatggcta                                                                                 <  1:270611/187‑8 (MQ=255)                                                                                                                                     acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATggctatgct                                                                             <  1:585531/191‑1 (MQ=255)                                                                                                                                     acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATtattat                                 <  1:566666/235‑1 (MQ=255)                                                                                                                                     acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTAtt                     >  1:266487/5‑251 (MQ=255)                                                                                                                                     acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTAtt                     >  1:572086/5‑251 (MQ=255)                                                                                                                                     acttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTAtt                     >  1:802420/5‑251 (MQ=255)                                                                                                                                      cttaaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTAttt                    <  1:829723/248‑1 (MQ=255)                                                                                                                                         aaGCGTTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTaaa                 <  1:844866/251‑1 (MQ=255)                                                                                                                                                aTCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGatattaac                                                                                                                                     <  1:74740/128‑5 (MQ=255)                                                                                                                                                aTCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGcagtgg                                                                                                   <  1:718475/162‑6 (MQ=255)                                                                                                                                                aTCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAaataatat                                                                                      >  1:122295/1‑173 (MQ=255)                                                                                                                                                aTCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCAtta                                     <  1:378780/224‑2 (MQ=255)                                                                                                                                                aTCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTAAATGGAGtt          >  1:289089/1‑251 (MQ=255)                                                                                                                                                aTCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCGGCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGATTGTTTCTTAAATGATGCAGTGGCAGTAaataatat                                                                                      <  1:796776/175‑3 (MQ=255)                                                                                                                                                  cAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGttt                                          >  1:784040/1‑217 (MQ=255)                                                                                                                                                        gCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCattattattgcc                             <  1:57228/224‑4 (MQ=255)                                                                                                                                                        gCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATCATTAATGGAGATTATCGTAATGCTGTTTCATTATTATTGCCTCAGTATTTAAATGGAGTTAGGGAACa  >  1:73577/1‑251 (MQ=255)                                                                                                                                                             ggAAGTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGcacgagaa                                                                                                                         <  1:422245/127‑8 (MQ=255)                                                                                                                                                                 gTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAGTGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGTATc                                                                      >  1:365875/1‑174 (MQ=255)                                                                                                                                                                 gTCGAAAACAGCTTATTAACAACATCTCCTAAAAACGGTGAAGTCTATAATCTTCGTGGTGTAACCTATGCTTTAATGGGCAATCCAAATAAAGCGAACGAAGATATTAACAAAGCACGAGAATATTTCTTAAATGATGCAGTGGCAGTAAATAATATGGCTATGCTAAGt                                                                         <  1:865588/171‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                  TTATAAATTATCTGAAATTTATTACAAGAAAGGGGACAGTAATTCATCTCTACTCTATTTACAGCCGTTGCTAACAAACAGTGGGCAGTTAAGTGAAAAGGCCAAAGTTCTACAAGCAAGAAATTTAAATCAACTCAAGCGCTATCAAGAAGCGTTGGAAGTCGAAAACAGCTTATTAACAACATCTCCCAAAAACGGGGAAGTTTATAATCTTCGTGGTGTGACTTATGCCTTAATGGGAAATCCTAATAAAGCCAACGAAGATATAAATAAAGCCCGCGAGTATTTTTTAAACGATGCTGTAGCGGTAAACAATCTTGCTATGCTAAGTATCGTCAACGGAGATTATCGTAATGCCGTTTCATTGTTATTACCTCAGTATTTAAATGGGGTTAGGGAACA  >  1/1486275‑1486676

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